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1. Introduction. Reynolds’ lubrication theory gives rise to a differential equation 
for the pressure distribution in a lubricating film. This film is usually restricted to lie 
between a stationary surface and a moving surface. In any practical problem, the two 
surfaces are, of course, of finite dimensions. However, all early solutions of Reynolds’ 
equation involved the assumption of one infinite dimension since this was equivalent 
to not permitting fluid flow in this direction, thereby reducing the problem to the solution 
of an ordinary differential equation. The results have severe limitations, however. In 
the case of journal bearings, the reduction of the axial length from infinity to a finite 
value can be expected to decrease the circumferential extent of the lubricating film, 
decrease the pressure developed throughout the film, and cause the point of closest 
approach of the journal and bearing surfaces to move closer to the plane of the applied 
load. These changes will in turn affect the performance characteristics of the journal 
bearing. For example, Shaw and Macks [1]** have compared properties of the full journal 
bearing predicted by the one-dimensional theory with measured values found by Brad- 


/ 
9-0, 


Fig. 1. Cross-section of a partial fitted journal bearing. 


*Received January 14, 1957; revised manuscript received April 8, 1957. 
**Numbers in square brackets refer to the bibliography at the end of the paper. 
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ford [2] and have concluded that the simple theory can lead to optimistic values of 
load capacity that err by more than 100 per cent. 

It is the purpose of this paper to determine an analytical expression for the pressure 
distribution in the lubricating film of a finite length partial fitted journal bearing when 
the lubricant is introduced under constant pressure through an axial bearing groove. 
A partial be: 
and the term “fitted”? implies that the radii of the journal and bearing are equal. A 


ring is one in which the bearing does not completely surround the journal, 


cross-section of such a journal bearing is shown in Fig. 1 where O and C are the centers 
of the journal and bearing, respectively, e is the eccentricity, U is the linea 
of the jour! il, and @is an angle measured from the directed line of centers Of The angle 


velocity 
6, is called the inlet edge of the bearing surface and @, is the outlet edge. 
2. Reynolds’ equation and its solution. or an incompressible and constant viscosity 


lubricant, Reynolds’ equation for a finite length journal bearing may be written as [3] 


1 a f Op i a (, op) Gull dh 





1 | : : ‘ 

06 \ 06 OY OY r dé 
where DP, M, id / re the pressure, viscosity, and radial thickness of the lubricant, ¥ 
is the axial variable measured from the central axial section of the journal bearing, 


and r is the common radius of the journal and bearing. If h is replaced by e cos 6 and if 


the dimensionless riable w yr is introduced, then (2.1) reduces to 
) OD \" . ‘a] ra] ) / cos 6 
(cost 0 22) 4 2 (cot 9 22) — tem? , 2.2) 
Og OG/ Ow ow ag 
where 8 Bul 
In order to remove the inhomogeneous term in (2.2) we set 
PO, W P(6, w) + f(6). ae 
Substitutir o this expr ion into (2.2) and eparating terms in @ alone from those in 4 
and w gives rise to two differential equations. The function f() satisfies 
d ( aft \ d Cos @ 
cos @ | re} 2.4 


dé ag dé 
and P(@, w) satisfies 


0 oP fa} : oP . 
( cos d } 4 . ( cos §@ ——] o. 9.5 


aa 
14 og OW Ow 
Equation (2.4) is simply the differential equation for the pressure distribution in the 
corresponding infinite length bearing and its solution is readily found to be 


es , N ae LK, [ dé _ K, 


Jy COS cos @ 2.6) 


Bg(0) + K,k(@) + K,, 


where 


g( 6) tan 6, 2.4) 

tan 6 | 1 + sin 6 

MG) « a——— + — beg i: (2.8) 
2 cos 8 } \l — sin 6 
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and K, and K, are constants of integration. The solution to (2.5) represents a term in 
the pressure solution which corrects for finiteness of the bearing or lubricant flow through 
the end axial sections. It may be determined by seeking a solution having the form 


P(6, w) = Q(0)W(w). (2.9) 


If (2.9) is substituted into (2.5) and the variables are separated, it is found that W and 


@ satisly the equations 


Ww’ —-NW =0 (2.10) 

Q” — 3tan 6Q’ + »°Q = 0, (2.11) 

where was chosen to be the separation constant. Equation (2.11) together with the 
mditions to be applied here is a problem of the Sturm-Liouville type and 


boundary 


es rise to a sequence of solutions Q; corresponding to a sequence of real, non-negative 


eigenvalues \? . We therefore write the solution of (2.10) as 


Ww. A; cosh (A,w) + B; sinh (,w) (2.12) 

f ~ (0) (X. + 0) and 
Wo A, + Bow (2.13) 
for 0 (A » 0). In order to solve (2.11) we make the substitution z = sin’@ 


which results in the equation 


1 5 ¥: 
~) a - — oe , NS o_ 

1 — 2)Q/’ + (; 9° Qi + 4 Q, = 0. (2.14) 

[his is the hypergeometric differential equation and has as its solution, for 7 ¥ 0, 
O.(f C.F(a, , b, : 3: sin? 6) + D, sin 6F(a; + 4, b; + 4; 3; sin’ 8) . 
' , 7 is idle (2.15) 

C.F5"’'(6) + D, sin 6F;°’(8), 

a; = 113 + (9 + 4)'”] (2.16) 
b, 113 — (9 + 4A%)'”?]. (2.17) 
Che above result is valid for z = sin*@ < 1 and this is consistent with the physical 
requirement that @ must lie in the open interval — 7/2, 2/2 for there to be a non-zero 


t all points between the journal and the bearing. When ¢ = 0 (A; = Ao = 0) 


Im thickness at 


orresponding solution Q, is given by 


6 dé 
; Sa" 
O c. i a D, 
(2.18) 


tan @ 1 1 + sin *) | ’ 
1 — log |———__] 1+ D,. 
cd 5 cos@ «4 log G — sin 6 iets 
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The general solution of (2.5) may therefore be written as 


P(6, w) = QWo + Dd Wiw)Q,(0). (2.19) 


t=1 
3. Application of boundary conditions. The simplest example of a journal bearing 

having an axial oil inlet is that where the lubricant is introduced at constant pressure 
Po at the inlet edge 6 = 6, . For this problem, the boundary conditions to be satisfied 
by the pressure are 

pla, ’ w) _ Po ’ 

p(6. ,w) = 0, (3.1) 

ple, Wo) _ plé, -— Wo) = 0, 
where w, = L/r and 2L is the total length of the bearing. In terms of f(@) and P(@, w) 


we require that 


{(0;) = De (3.2) 
{(@.) = 0 
and 
P(6, , w) = P(@, , 0) = 0 — 
(3.3) 
P(6, wo) = P(@, — wo) = — f(A). 


With the use of these boundary conditions, we can determine the coefficients K, , Kk. , 
A,,B;,C;, and D, as well as the eigenvalues \; . From (2.6) and (3.2) it is easilyyshown 


that 


: [ if (@,) — 6,) 
K, = Sig 2) — K4)] + p (3.4) 
k(6.) — k(6,) 
and 
B[g(0,)k(0.) — g(02)k(@,)| — pok( 62) 3 5) 
— en maaan _ . (3.0) 


K 
" ke G2) —- k(0;) 


In order to satisfy the first two conditions of (3.3) we require that each Q,(@) satisty 


Q,(6,) = Q,(6.) = 0. (3.6) 
From (2.15) and (3.6) we have the two equations 
C,F‘"(0,) + D,; sin 0,F;"'(0,) = 0 (3.7) 
and 
C.F‘ (0.) + D; sin 6.F;"'(6.) = 0 (3.8) 


which may be combined to give 
C {sin 6,F‘'(0.)F.?(0,) — sin 0.F\'’(0,)F;"(6.)] = 0 (3.9) 


and 


D, {sin 0,F‘(6,)F{?(6,) — sin 6,F‘?(0,)F\(0,)] = 0. (3.10) 
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Since C; and D; cannot both be zero for a non-trivial solution we require that 
sin 6,F§"(0.)F;”(@,;) — sin 6.F$"(0,)F;?(0.) = 0. (3.11) 


This is an equation for the eigenvalues \; . The use of this result is reserved for Sec. 4. 
With the dj determined, it is now possible to reduce the two unknown constants C; and 
D, in (2.15) to only one unknown by making use of (3.7) or (3.8). In order to determine 
C, and D, we apply (3.6) to (2.18) to get the two equations 





6, 
Col agra t Do=0 (3.12) 
0 cos 6 
and 
™ @ _f° & o 
Co i oar a + Co [ cos? 8 + D, = 0. (3.13) 


Subtracting these expressions gives C, = 0 after which D, = 0 from (3.12). Consequently, 
there is no solution corresponding to 7 = 0. 

The remaiiing coefficients may now be evaluated by making use of the last two 
conditions of (3.3). From (2.12), (2.19), and (3.3) we have the two equations 


: [A,; cosh (A,;wo) + B; sinh (A,;w,)1Q,(0) = — f(A (3.14) 
and 


tC) 


>» [A, cosh (A;w.) — B, sinh (A,;w)1Q;(8) 


t=-1 


— f(@), (3.15) 


where it is to be understood that the single unknown coefficient in Q,;(@) has been in- 
corporated into the A, and B; . Before employing these expressions, we first note from 
2.11) and (3.6) that the Q,;(6@) are orthogonal with respect to the weighing function 
cos 6 over the interval 6, , 6, . That is, 


Os 
| cos* 6Q,(8)Q,(8) d@ = 0 (3.16) 


for i + j. If (3.14) and (3.15) are multiplied by cos*@ Q;(@) and the results integrated 
over 6, , 6, two equations in A, and B; result which, when solved simultaneously, give 


a@s 
| {(9) cos’ 6Q,(8) dé 
By eee -- (3.17) 


—_— ; 
cosh (2, w,) | cos® 6Q2(6) dé 
6, 


and 
B,; = 0. (3.18) 
The latter result makes the pressure solution an even function in w which is to be expected 


from the symmetry of the problem. This completes the evaluation of the coefficients, 
and the pressure solution 


p(0, w) = (0) + 3 Wy(w)Q.(6) (3.19) 


t=] 
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is now completely determined for the case where the lubricant enters at constant pressure 
through the inlet edge of the journal bearing. 

For the case where the lubricant enters through an axial bearing groove at 6, where 
6, < 0¢ < 6, the analysis is analogous to that preceding except that two pressure 
solutions of the type (5.19) must be developed, one for the range 6, , 9¢ and another for 
the range 6, , 6. . Once the pressure distribution has been determined, it is possible to 


compute such quantities as load carrying capacity, coefficient of friction, axial thrust, 


and rate of flow of lubricant. These will not be discussed here, however, since their 
definitions are readily available in any one of numerous lubrication texts. 

4. Evaluation of eigenvalues. In the preceding development, (3.11) gives the 
eigenvalues \? . In this section one technique for determining them will be described 
which employs Whittaker’s method for finding the root of smallest absolute value of a 
power series [4]. 

We first recall that the hypergeometric function may be expanded in the form 


oe a b a h a +- l b + l ° 
PG; Dees l + S++ 3 in 1.1) 
L<¢ le 2 c+] 
The functions F'’?(@) and F;°’(@) are easily expressed in this form if we note that the 


coefficients 11 1.1) eontan terms of the type 


a+n)(b 1. 9 ab + nia — b rn, LZ 
where v is zero or any positive integer. From (2.16 and (2.17) we have that ab N« i 
anda —_ r 2 so th t. tor re (@). 
A n 
+ ny(b +n - - — (Zn o 1.3 
1 5 
Similarly, for F;°’ (6 
v ] 
( ? D+ — +. Y\(D, { 
n > \n o)lZ ] 1.4 
The coefficients in the « ‘pansions of F 9) and rs (?) are therefore polynomials in 


> 


A; . These expansiol ; may now be substituted into (3.11) and the series rearranged 
to give a power series in \; where each coefficient is, in general, a rapidly converging 


infinite series. If we represent the result a 


then, by Whittaker’s method, the smallest root Xj is given by 


1 
é ( é 
Cn ¢ é ¢ é C5 
- é = = —_— 
é é ( é é L.6 
é é é ¢ : 
8) é €; | 


The second root may now be found by factoring out the first root and repeating the 
procedure, and similarly for all other roots. 

5. Asymptotic forms. Because of the complex nature of the eigenfunctions and 
Eq. (3.11) for the eigenvalues, it is convenient to develop approximate forms which 
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are valid for large values of \; . Watson has derived asymptotic expansions for F(a, b; ¢; x) 
when various combinations of the parameters are large [5]. When applied to F;" (6), 


the appropriate expansion gives 


2 cos [30(9 + 4d3)'?] TL + {9 + 4a5}'?) (5.1) 


FY? (0) ~ ar Le 2) 172 
: (9 + 4X2)" cos”? 9 TG[—1 + {9 + 4n7}"7) ’ 





where T° represents the gamma function. If all functions in \; are expanded and only 


terms of order 1 are ke pt, (5.1) reduces to 


: 08 A; 6 ' 
re ~ SSe. (5.2) 
cos 6 


Similarly, it can be shown that 


sei 4 sin \;0 - 

Fi; (0) ~ > ns a (5.3) 
A, sin 6 cos i] 
If these results are substituted into (3.11), it is found that 

lm re 

A, ~ he (5.4) 
6. — A 
Che sume result can, of course, be determined from the general theory of the asymptotic 
behavior of the eigenvalues in the Sturm-Liouville problem [6]. 


6. Example. In order to give some indication of the nature of the solution deter- 


mined above, hall consider the simple case where the lubricant enters under constant 
} 
pressure at @ = @, and we shall suppose that 6, —6, . That is, the clearance between 
the journal and the bearing surface will be taken to be symmetric about the line of 
enters 6 0. In addition, we shall take @, to be small. This corresponds to a small 
bearing are and makes \, large from (5.4), thereby justifying the use of the asymptotic 
torms deve loped wove. 
Krom (3.4 d (3.5) we have that 
. 28q(0) + Do _ 
kK, 7 oe (6.1) 
2k( 6.) 
Po 9 
(6.2) 
K 2? ) 
so that, from (2.6), 
28q(02) + po Po ‘ 
f(@) = Bq(@) — ——— "+ (0) + (6.3) 
2k( 02) 2 


It will be recalled that this is not only a portion of the finite length bearing solution, 
but is also the pressure distribution in the corresponding infinite length bearing. 
From (5.4) the square root of the first eigenfunction is given by 


(6.4) 
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and (2.15), (5.2), (5.3), and (6.4) give 
cos [x 26.) | 26. sin [7r6/(26.) 
0,10) ~ C —75 CSOD =. (6.5 
: . cos’ ’~ 6 wr cos” 6 
From (3.6) it is found that D, = 0 so that 
cos [76/(26.) = 
@),( 6) ~C ne ; ¥ . (6.6) 
: cos ‘~ ¢ 
Similarly, 
| 
6, sin (76/6.) ” 
Q.(0) ~ D —353— (6.7 
tT COS 0 
and the pressure solution from (3.19) is 
9Anvl ft i 
Z200\ G5) _ Pp p 
} 3qa(@) — | = a 
pe, u = ee ies 91-( a I 6) i 9 
St (6.8 
( us \ cos [76/20,] , 0. x \sin (70/0,) 
+ A, cosh |=— w]) ——5—— + A, cosh w ) —a ——— + 
26, cos ’~ @ i 0 / cos a] 
where g(@) and h(@) are given by (2.7) and (2.8), and, from (3.17), 
— . : r 6 P 
A, = eal I : | cos ~ 6 cos v dé (6.9 
- \ zO6 
a“ e - 
6, cosh (<= u 


and 

















yt 
O 


>O 
“15 A0 -5 re) 5 10 15 DEG 


Fia. 2. Circumferential variation of pressure. 
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—29 
6, cosh (rwW/ 62) 


6, 
/ [8g(6) + K,k(6)] cos”? 6 sin x . dé. (6.10) 
0 2 


The integrals in these expressions are easily evaluated by successive integrations by 
parts if we recall that 7/20, was assumed large. 
We shall now apply (6.8) to the specific problem where 


radius of bearing = r = | in. 

length of bearing = 2L = 1 in. 
speed = 2000 r.p.m. 

eccentricity = e = .001 in. 

viscosity = » = 1.2 X 10°° reyn 

oil inlet pressure = py = 30 Ib./in.’ 
outlet edge of bearing = 6, = 15° 


With these conditions, the pressure solution to two eigenfunctions reduces to 


~| , tan 0 l sin 6 " 
p(d, w) ~ 1510 tan 6 — 3s7| 2 — + log (Ltsin 6) ] + 15 


cos 6 1 — sin 6 (6.11) 
0 cosh (Gw) £28,68 5o ¢ » Sin 126 
— 1.80 cosh (6w) cos’? 8 + .0452 cosh (12w) cos’ 6 


This result can be expected to represent a good approximation to the complete pressure 
solution for all w not close to w. = + 3 since, in these neighborhoods, the hyperbolic 
functions in (6.11) are relatively large. However, in these cases, the behavior of the 
pressure solution is known from the boundary conditions. 


Equation (6.11) appears graphically in Figs. 2 and 3. Figure 2 illustrates how the 


p 
Ib/in* 


— ° 
30 62-15 





25 


20 











+ > tW 
0 1/8 \/4 3/8 72 


Fra. 3. Axial variation of pressure. 
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pressure distribution at various axial sections differs from that for the corresponding 
infinite bearing. While the distribution at the central axial section w = 0 does not differ 
greatly from the infinite distribution in this case, the difference would be considerable 
had the bearing length been taken shorter or the bearing arc larger. Figure 3 shows the 
variation in pressure as a function of axial distance for various fixed values of 6. The 
pressure decreases from a peak value at the central axial section to zero at both end 
sections. The corresponding pressure for the infinite bearing would be, of course, constant 


for all w and each fixed @. 
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A CLASS OF PROBLEMS ON LONGITUDINAL VIBRATIONS* 


BY 
VACLAV VODICKA 


Plz n, Czechoslovakia 


1. Introduction. 


Consider a bar of length 1 and constant area q formed of n parts 
vith le ngths l densities Py and Young’s moduli E, (k = I, 2, * , BR). Denoting by 
the displacement in longitudinal vibrations at the point x of the kth part and 
at any time /, the general form of the Laplace transform 


l °,p) =p | e”'u,(a, t) dt; k=1,2, +++ 2 
) = 4, cn > aks Whe aes be BR a 
72 = ‘ é. Kv, Pp), k Vp.’ \ >=) y fs 
Phere exists a vast class of cases, important from the technical and physical points 
O ew, where 
i D i} a F Wills, , p) e,Wiails 5 p)3 
‘ ies , _ aW, (2) 
oe ( Wy —= k= 1,2,--- - f. 
a KE.’ k ag ’ n 
(ll problems of this gr 


group may be treated in the following manner. 
2. General theory. 


The subsidiary equations corresponding to the well-known 
conditions at the 


sections separating adjacent parts of the rod are (primes denote the 
-derivatives): 
Pp l i UL(s » P) e.Ui.,(6 P) k=1,2 


,2,---,n—1. (3) 
Substituting from (1) and using (2) we have 


. Ss - 8, 
| cosh § + B, sinh # A, 
qd a 


IS, . IS; 
+1 cosh I te B, , sinh I , 
a a 


k+1 


+1 
: IS Ss ‘ IS; > IS). ‘ 
{, sinh § + B, cosh © h(a , sinh = + B,., cosh Bes.) (3.1) 
a a Gas a 
Apa k+ 
hy wo —etieet . k=1,2,---,n-1, 
Ay. pi. 


ived April 16, 1956; revised manuscript received July 26, 1956. 
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or, in the matrix form 


Aves |, M.(p) | a oe ee 


Ss IS; a ae IS 
cosh = cosh __ sinh — sinh 4 
a a h, a , 
, S S IS, IS 
sinh f cosh P cosh Prt sinh p 
v P a Qeei a, ) are 
Mp 
8 — 8 
cosh — sinh + aa sinh cosh © 
a a, ly a a 
; Oy. ps l 8 8 
— sinh — sinh £ t cosh — cosh 4 
a a 1 } a, a, 
k ee n— 1 
It follows from this that 
=F 
= (),(p) a | eam 1.2, 2-2 n— I, 


B B, 
VM (p)M,-,(p) «-- M.(p)M,(p); k i ee oe 


Q,(p 
Writing 
Q),.(p 1 (p My 1,2 i i n— 1 
we obtain from (4.1 
1, pa "'(p)B, , B q. (pA: + qe (p)B, ; 
k=1,2 


and the expressions (1) become 


cosh & | B, sinh =. W (a, p 
a, 


a 


, ee Af y sl pr 5, ODE al } Pr { 
PP Z,D : E p) cosh = r Qk i aoa |. 


k+ 


a 
k l,2,*** ,"m = 1. 


Transforming the conditions at the ends of the bar leads to the equations 


a, U,(0, Dp) +P Us(O, P) Yo; 
a Ut. P) +P At. Pp t. . 


n 


where a, , Bo, *** » Bn» Yn Gepend usually on p. 
By means of (6.1), (6.2) relations (7) yield a syste 


| 2 ° IX , 
a] iy cosh P r @ Pp sinh = |, oa WW k+iht, DP); 


(Vol. XVI, No. 1 


(4) 


(6.1 


m of two equations for the un- 
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knowns A, , B, . Solving it and substituting the resultant values A, , B, into (6.1) and 
(6.2) leads to the following complicated expressions: 


eo ne (4, cosh 2 + A, sinh ve) + Wi(z, p) (8.1) 
A a, a, 


U,.:(2, p) = 1 ya af" (p) cosh + G ze | 


+ a4] a" (p) cosh = + q.” (p) sinh abe + Wi+i(2, p); 


= 1, --,n—1 (8.2) 
A | («. cosh pl + 28 sinh Bl) 223) + (a, sinh pl 
a, a,, a, a, 


+ . 8, cosh a 200) lv — aW,(0, p) — BWi(0, p)] 





7 = BolYn = a,W,(L, P) ad B,W (1, p)] 


1 


A aly, — a,W,(l, p) — B,.Wi(1, p)] — | («. okt 


n 


+ . 8, sinh ea) 4. («, sinh D 


n 


+ 4 B, cosh ol) 29) |r =— aW,(0, Pp) — BoWi(0, p)| 


) naa) l . l 
A= E 10-(p) — - Bod » 1) |(a, cosh o + “* 8, sinh et) 


+ | ana 220) _ - aa‘) |(a sinh e + = 8,, cosh et), (8.3) 
The formulae (8.1) to (8.3) determine the Laplace transforms of the functions 
(x, 1) which are to be found by the methods of operational calculus. The most difficult 
part of the problem lies in finding an explicit form of the matrix (4.1). 
3. An interesting class of problems. We now proceed to treating problems 
characterised by the equations h, = 1 (k = 1, 2, --- ,n — 1), or, by (1) and (3.1), 


E.p, = const; g = 7,3, ++>., 8. (9) 


In this case it is possible to calculate the matrix (4.1) and our former results (8.1) to 
(8.3) become simpler. Physically the simplification arises from the fact that there is 
no reflection at the joints of the different parts. 

Using the notation 


a-a(t-—-), n=Yes ka1 nal (10) 
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we easily obtain from (4) and (4.1) 


cosh po; ; sinh po 
M Ap I I : I 1,2,°*- ,n—I] 
sinh po, , -osh po; 
Por CC 17 (11) 
S ! 
cosh pro , sinh p72 || 
Q.(p . 
sinh pr. , cosh pro 
However, 
cosh PT, sinh PT, cosh PTv+1 5 sinh PT v+1 
M,.:(p 2 
sinh PT» ; cosh PT, sinh PT v+1 » cosh PT, 
so that (4.1 vie lds the closed forms 
cosh Tk y sinh IT 
Q,.(p iii f ; k 1.2,°::  n— I (12 
sinh PTk cosh pT 
of the fundament: | matrices ()).(p ; 


Equations (8.3) then give 


, f ] l : 
A, = P sinh o\ ) - P Bn cosh of ae oT Je — a,W,(0, Pp 
3 a 


a A, 


- B,W10, p P B.ty, — aW,(l, p) — BW, p 


a 
A; AolY¥n — @ WD = B,Witl, p E cosh p{— + 7, | 


) l ' 
P 8, sinh ol + 7 ) Jo - aW,(0, p) — Bo Wi, p 


a, a 


/ / / . 
y . l \ AnD, a,,0 ( l 
i 2 q ae ~ Sook -— - ‘OS T g >) 
A= # Qa, See ) inh pi a. 7 } rj pe r- cosh p a. 1 } 
and the formulas (8.1), 8.2) become 
1 J ; . fl a 
[ Pe, A $14 a,W (0 Pp BW 3(0, P E sinh p\- + Tn-1 a } 
4+ PB cosh (! } 7 Cy) + fy. —a,W.(l, p) — BWI, 
- > I . 7 a,/ ! I >| Mr rf 
- - - nh = _— P b cosh e) le W, oe Pp ’ \ | . J 
a ay, a J 
, ; ; | a x 
l L,p - a0 aoW,(0, p) — Bo WO, p | sinh \ - ae 
\ Pp l A ' r r 
- fT, — tT) 7 3, COS = - 7 — — ivy — &@ (l,p 
A Tx | ” co h Pp a. " | .)] T y a 
2/7 ° i 1 x = pP 2 . m j x 
— B,Wi(L, p | sinh of: 5 ae 4 4, BO cosh aC - )|} 


+ Wiel, pd); &=1,2,+%+,n-1. 
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[t is easy to verify that the expressions (13) to (14.2) satisfy in fact the boundary con- 
ditions (7 

Formulae (14.1), (14.2) represent the Laplace transforms of u,(x, t) and from this 
we obtain in each concrete case the completed solution by the methods of operational 
calculus. Several examples are given in the following section. 

4. Special problems. We illustrate the preceding theoretical considerations by 
typical examples generalising the results given occasionally in the literature [1]. In order 
to have reasonable expressions we assume (9) to be valid. 

All the caleulations needed in each case proceed by five stages. First, one applies 
the Laplace transformation to the fundamental equations governing the problem and 
to the initial conditions. This yields the general form of U,(x, p), as indicated in (1), and 


especially also the expressions W(x, p). 
When the relations (2) hold, we write down the subsidiary equations corresponding 
to the boundary conditions at « = 0, « = land obtain the values ag , By , Yo , Gn » Bay Yn + 


The third and fourth steps consist in calculating the expression A by (13) and the 


Laplace transforms U,(x, p) according to (14.1), (14.2). Finally, we deduce the wanted 
solution , !) by the methods of operational calculus. 
4.1. Vibrations of a bar under its own weight. The upper end x = 0 is fixed, 
the lower ! supported so that the displacement is zero at all points. At ¢ = 0 the end 
lis released. 
Putting s (0 fundamental equations of the problem with belonging initial conditions 
ure 


Ou; Ou = 
a + 9, iji< <2 Sk t> 0; k=1,2,-*> ,n (15.1) 


ol ox” ia - 
du,(x, 0) 7 
u(x. 0 0, —— = 0, mi <2 Ce; k= 1,2, --> a. (15.2) 
c 
Solving the corresponding subsidiary equations with respect to U,(x, p) and comparing 
with (1) gives 
=, qg - 
W(x, p) = 9? . k= 1,2, --° ,n. (15.3) 
) 
Thus, the relations (2) hold. 
The boundary conditions at x = 0 and x = Il are 
du,(l, t) 
u(0,) =0, “—“=0; t>0 (16.1) 


OX 
and this leads to 
U,(0, p) = 0, Us(l, p) = 0. 


We therefore have 


aA, = B,, = r; Bo = V7 =a, —- %2 = 0 (16.2) 


and (13) gives 


A = © cosh r(_ + tat). (16.3) 
a Qn 


n 
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we easily obtain from (4) and (4.1) 


cosh po; , sinh po; | 
Mp I Perl. =k = 1,2,+++,n—1 
sinh PCr ; cosh po (11) 
-osh prs sinh pr. 
Q. p ” h PT ; PT2 
sinh prp , cosh p72 
However, 
| “ ! . ' 
cosh PT» sinh PT, cosh DPT y+1 5 sinh PTv+1 
M,.+1(p - 
sinh pr, , cosh pr sinh pr,+: , cosh pr, +; 
so that (4.1) vields the closed forms 
cosh IT i sinh DT 
0.(p sil eh ba 1S +s gel (12) 
sinh PTk ; cosh PT 
of the fundamental matrices Q,(p). 


Equations 8 3 then clive 


l ‘ 
A, = E sinh of ge: ) + P B,, cosh nes + Tx Ie ~ aol i(Q, P) 


a a,, 


- BoWi(0, p)] — © Bol, — anW,(l, p) — B,W(1, p 


if 


rs 
A, = aly, — a,W,(l, p) — B,WiCL, p)] — E cosh pl +7 ] 


a, 


) l = 
8, sinh of + 7 | — aW,(0, p) — BW, p 
a 


Dp l (a re) A,B \ l 
— : i > as ] ( eo } j _—_ r , s] 4 a 2) 
= Yn, On0 lI Ve ( 1 7 . 
- (. 1, =e a, e a a ree a, . 


( 1 
and the formulas (8.1 8.2) become 
U,(a 2p W,(0 ~ BW, p)]| a, sink €: + 7 ~2) 
1\0, DP A y Qot » P) 2) »pP ; a il p rs Tn-1 ay) 
I x\] : e 
F 8, cosh o + 7-1 }| + fy, — @.W,(1, p) — B.WiL, p 
a a a 
. (a sinh 22 — 2 g, cosh @)> + W,(z, p), 14.1 
a a; a ) 


; l x 
U44:(2, p z 1 aoW,(0, p) — BWi(0, p | a. sinh pf - — 
a, a, 


+ t-1 — tT) + F 8. cosh ol — + T,-1 — 7, )] + lyn — a, W, (1, p 


a \a, a 


— BW(l. nv - sinh T. 7 Z ) ane 8, eos (; -+ . > 
7 L I | re. 3 Fs cosh Pp aa +) J 


- (14.2) 
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[t is easy to verify that the expressions (13) to (14.2) satisfy in fact the boundary con- 


ditions (7 
(14.2) represent the Laplace transforms of u,(x, ¢) and from this 


Formulae 14.1 
we obtain in each concrete case the completed solution by the methods of operational 
calculus. Several examples are given in the following section. 


4. Special problems. We illustrate the preceding theoretical considerations by 
tvpical examples generalising the results given occasionally in the literature [1]. In order 
to have reasonable expressions we assume (9) to be valid. 

\ll the calculations needed in each case proceed by five stages. First, one applies 
the Laplace transformation to the fundamental equations governing the problem and 
to the initial conditions. This yields the general form of U(x, p), as indicated in (1), and 
especially also the expressions W(x, p). 

When the relations (2) hold, we write down the subsidiary equations corresponding 


to the boundary conditions at 2 0, « = land obtain the values ap , Bo , Yo, @n 5 Bn y Yn + 


Che third and fourth steps consist in calculating the expression A by (13) and the 
Laplace sforms U(x, p) according to (14.1), (14.2). Finally, we deduce the wanted 


solution u,(v, 2) by the methods of operational calculus. 


4.1. Vibrations of a bar under its own weight. The upper end x 
upported so that the displacement is zero at all points. At ¢ = 0 the end 


= () is fixed, 


the lowe! 


11s rele: sed 


Putting (0 fundamental equations of the problem with belonging initial conditions 
re 
ra) Ou, , , ~ 
a + 9, ina XS Ce, i> 0; k=1,2, +--+ ,n (15.1) 
; og 
du,(a, QO) - 
0 0 “a = 0, ia <2 <€e@2 k=1,2,--+,n. (15.2) 
c 
Solving the corresponding subsidiary equations with respect to U;,(x, p) and comparing 
ith (1 fives 
: q _ 
W(x, p) =; k=1,2, +> ,n. (15.3) 
p 


Thus, the relations (2) hold. 
The boundary conditions at x = 0 and x = lare 


u,(O0, t) = 0, dun(l, t) = 0; t>0 (16.1) 


OX 
and this leads to 
U,(0, p) = 0, Us(l, p) = 0. 
We therefore have 
a = 6 = 1, Bo = ¥o @ & = Ye = 0 (16.2) 


and (13) gives 


A = E cosh r( t ad me | (16.3) 
a, An 
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It then follows from (14.1), (14.2) that 
\ 


l x 
cosh r(- + Tr-1 — } 


\G, a 


U,(z,p) = 4/1- 
onl ( l oe 
ce 1p a, Tn 


l x 


% | 
gd \a,, Ay +1 4 4 
t- M2__fee —— | hie i, Sse 2 OD 


P cosh »(— +r ) 


and from this we find (with the aid of a table of Laplace transforms or by contour 
integration) the complete solution in the form 


gz il x 
u(z,)j == (! + 7,.5—-=— 
a, \d, : 2a,, 


16q (1 2. {2 r 
~ (- " teil 2 Ont ( l oe 
7 Tn-1 


Te ( \(_! x) l6g £ y 
u 1 2 a + FT Tins = a 1 = 
ci a \a, - |e 2 x \a ‘ 


# x 
yi ) - 
— (—])" m+ 1ye( e % Mees " (2n + 1)xrt 
° ba —a_e ] dean Ceeaeeieeemmmmee ker } ; 
: ts fie) 


a eg a, 





k 12, --- nm — l. tee: 


4.2. Bar with a mass at the end. The end x = 0 is fixed, a mass m is attached to 
the other end x = 1. The bar is initially stretched by the axial force gS and at t = 0 
its end x = | is released. 


The fundamental equations 


An, OU 
3 a = ; i = eS h, t > Bs b= 1,2, --> 1 18.1) 
ol oO 
with 
/ tat ] r- s : Ou r. U 
u(x, 0) = sl — + - — ), mn = VQ, tus eS ee 
\ 2, E. EB, al 
b= 1, 2. , ft 18.2) 


and with the boundary conditions 
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u(O = 0 
0°u, (1, t) , ou,(l, é) = du, (1, 0) = . 
eS — 9B, Eu = 8 De, a = 03 «E> O (18.3) 
give, as In the preceding case 
; k-1 l. — &-, 
W,(4 S\ Ug + =e a =1, Bo = Yo = O,7 
c=1 4/¢ “k (19) 
a, = a »s = gE , = S ; Le 
mp — 
Thus, we find without difficulty 
A = sinh rf te rat) - = cosh (_ + 7, ) (20) 


and formulae (14.1), (14.2) yield 





seteap = ($ehen tz) (21.1) 


sinh of 7 + 2) 


qS ae p+ 


m p| p sinh (2. + 7,- ried + ss cosh (4 + T,- | 


k=1,2,--- .n—1. 








The solution is , 
Sx n E 2 S l 
u(r, tl) = E. (1 -_= 1) + 245 ( de. rat) 
‘A a, E, m \a, 
5a 
: a, 6,t 
— —3 oO OE 


/1 a l l 
5 ( l i ~ ) gE, + (! + r, J (2: ab) + i | — a Fane — + Ta-1 
| he. Ma, a, . ma, a, a, 


a s a,,T x a, E;,, 
sy i ia, eee ale | ee (1 — Hs) 
= | & E, | Ore E, t Es Ak+1 E, 


2qS [ l ' sec 6, 
[@+~)e+ G+) eee 


ma, a, ma, 


sec 6, 











a,(+ 2) . 
Ars? ogg Of = 1,2,02,m— 1. (21.2) 
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) are the positive roots of 


The 6,(p = 1, 2, --- 


l gE, 
6tan 6 = & + Ta-1 . 
\An Ma, 


Two equal rods of the kind treated above move longi- 


4.3. Collision of two rods. 


[Vol 


XVI, No. 1 


(21.3) 


tudinally in opposite directions with equal speeds v and collide at ¢ = 0, at the origin 
xz = 0. By symmetry we consider only the rod 0 < x < 1 and so long as the rods are in 


contact we have the equations 


Ou; 5 Ou 


—— = @. — Bea eS ee | 
ot OX 
with the initial conditions 
—_—" Ju,(a, 0) 
u,(z, 0) = O, =a —= —f i. <2 es 
ot 
and with 
yu, (Ll, t) 
u,(0, t) = 0, — = 0), 
OX 
From this we obtain, as before 
W,(z, p -- l, 2 nN); @ B,, l 3 
p 
) l 
A = i, cosh p( “p T, ), 
a a y 
and formulae (14.1), (14.2) become 
l 
cosh of + T 
] \ad 
[ p | 
Pp / 
cosh p - T 
a 
l 
cosh P| 
’ a 
- Pp l 
p 
cosh | 
a 
Hence, the solution is 
] 
—r | Qr + 1)a| 
7 j —s { al 
u(r. ft 1. + }- ps COS 
7 \a 2r+ | ] 
9 
A" 
re) / 
ri t ( + T . 
(! | 
‘ 2r + 1)r + T T 
| a a 
2 +T } 
\a 
k 


L. 2, n 
Ries n 
a = ¥, 0 


sin 
9/ 
~| 
(2r + I)rl 
l 
(14 ns) 
a, 
nwo I 


i? “ee 
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5. Final remarks. For n = 1 the formulae for u,(z, t) in (17.2), (21.2) and (24.2) 
change into the well-known elementary results given in the literature [1]. 

Our preceding deductions may serve as a kind of guide in treating other problems 
on vibrations of composite continua. A special example has been discussed in detail in 
the author’s paper given in the literature [2]. The method of matrix analysis can be used 
also in treating other important questions of engineering such as torsional oscillations 


[3], conduction of heat [4, 5], etc. 
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BOOK REVIEWS 


Plane waves and spherical means applied to partial differential equations. By Fritz John. 
Interscience Publishers, Inc., New York, 1955. viii + 172 pp. $4.50. 


The present tract contains a collection of various results on partial differential equations, their 
unifying thread being the use of certain elementary identities for plane and spherical integrals of an 
arbitrary function. The author’s principal aim is to show that many results on fairly general! differential 
equations follow from those identities. The book begins with a concise introduction, which illuminates 
the paths to be followed. Chapter I deals with the decomposition of arbitrary functions into functions 
of the type of plane waves, i.e. into functions having parallel planes as level surfaces. The decomposition 
used here consists of expressing a function by spherical means of integrals over hyperplanes, and is due 
to J. Radon. This kind of decomposition is to be compared with that furnished by Fourier analysis 
(into plane waves of exponential type). The ‘“‘Randon transform’ has the advantage over the “Fourier 
transform’”’ that the integrals appearing in it involve the given function itself, and not its Fourier trans- 
form. In Chapter II the Randon transform is applied to the solution of the Cauchy problem for the 


homogeneous hyperbolic equation with constant coefficients L(u) = Q(0/dx, , +++ , 0/dx, , d/at) u = 0, 
where Q(m , °** , m, A) is a form of degree m in its arguments, with constant real coefficients; with the 
initial data (d*u/dt*) so = f(a, , *** » In) fork = 0,1, -+- ,m — 1. Chapter III gives a construction of a 
fundamental solution for a linear elliptic equation, and for a system, with analytic coefficients. The 
problem amounts to finding a solution of the symbolic equation L(u) = 6, where 6 is the Dirac function; 
and the method of solution employed is that of decomposing the 6 function into plane wave functions, 
thus reducing the problem to that of solving L(u) = f, where f is a plane wave function. In Chapter IV 


are found formulas for an arbitrary function in terms of spherical integrals of the function. These formulas, 
which form the principal tool used in the later chapters, can be regarded as generalizations of the ex- 
pressions of Chapter I, which give the arbitrary function in terms of its integrals over planes. They 
can also be viewed “as the analytic counterpart of the geometrical fact that spherical shells can be swept 
out by spheres in two different ways, just as the identities of Chapter I are connected with the fact that 
the exterior of a sphere can be swept out by planes’’. Chapter V is concerned with the mean value theorem 
of Asgeirsson for the ultra-hyperbolic equation, together with a somewhat more general identity, having 
to do with the ellipsoidal means of a function, due to A. S. Howard. Chapter V deals mostly with the 
problem of the determination of a function from a knowledge of its integrals over a sphere of fixed 
radius, which can be solved by means of the identities of Chapter IV. Chapter VII contains the major 
application of the identities on spherical means in Chapter IV, to the proof of the differentiability of 
solutions of linear or non-linear elliptic equations or systems, provided that the coefficients are sufficiently 
regular. Chapter VIII extends the results of the preceding chapter to non-linear elliptic equations, save 
that here the regularity not of the solutions themselves, but that of certain integral transforms of them 
is established. The lucidity of the exposition combines with the attractive format of the book to make 
it the delight of mathematicians interested in the subject. 
J. B. Draz 


Recent advances in science (Physics and Applied Mathematics). Edited by M. H. Shamos 
and G. M. Murphy. New York University Press, 1956. x1 + 384 pp. $7.50 


The book is an outgrowth of the First Sy mposium on Recent Advances in Science held at New York 
University in the spring of 1954; it contains the following articles: Methods of Applied Mathematics 
(R. Courant); The Future of Operations Research (P. M. Morse); Atomic Structure (I.I. Rabi); Micro- 
wave Spectroscopy (C. H. Townes); Nuclear Structure and Transmutations (H. A. Bethe); Elementary 
Particles (V. F. Weisskopf) ; Electronuclear Machines (L. J. Haworth); Neutron Physics (N. F. Ramsey) ; 
Transistor Physics (W. Shockley); Ferromagnetism (R. M. Bozorth); Cryogenics: Very-Low Tempera- 
ture Physics and Engineering (F. G. Brickwedde); Physics and the Engineer (E. U. Condon 


(Continued on p. 32) 
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APPROXIMATE SOLUTIONS FOR A CLASS OF INTEGRAL EQUATIONS* 


BY 
RICHARD LATTER 
The RAND Corporation 


1. Introduction. Exact solutions for integral equations of the type 


f(z) = g(z) +2 [ dyk(x—y)fy), O<2<a, (1) 
70 

are, in general, exceedingly difficult to obtain, except for particularly simple forms for 
k(x — y) such as a sum of polynomials times exponentials. In the present discussion, 
it will be shown that by an extension of the Wiener-Hopf technique, it is possible to 
obtain approximate solutions as well as approximate eigenvalues for this equation. 
The accuracy of the solutions will be shown to improve exponentially with increasing 
order of approximation. Application of the method is made to the cases wherein k(x) = 
e'*! for which case an exact solution is obtained, and k(x) = 1/(2r)'” exp (—1/2z’), 
for which case the lowest eigenvalue is obtained approximately. 

2. Fourier transform of (1). The first step in the present method consists in obtaining 
the generalized Fourier transform (Titchmarsh [1]) of (1). For this purpose, it is assumed 
that f(x) = O(e°'*') (ce < 1), g(x) = O(e*'*') (d < 1), k(x) = O(e™"*'), and k(x) = k(—z), 
k(x) real; and that these functions belong to L,(— ©, ©). The following arguments 
treat (1) as if it holds for all x. This constitutes no limitation since if (1) is assumed to 
apply only for 0 < x < a, then the conditions k(x) = 0 for | x! > a and g(x) = 0 for 
xz < —aand for x > 2a may be arbitrarily required. From these conditions it is observed 
that f(z) = Ofor2 < —aand forz > 2a and therefore f(x), g(x) and k(x) automatically 
satisfy the assumed growth conditions. To simplify the subsequent formulas, the following 
quantities are defined and their region of analyticity indicated: 


(7) f.(w) = | dx e*“ f(x) (Imw > oe), 
(77) f(w) = dz e*" f(x) (Imw > — 2), 
70 
(7it) f(w) = | dx e*" f(x) (Imw < —o), 
iv) Kw) = | dx e*“k(x) (—1 < Imw < 1), 
v) g.(w) = | dx e'*"g(x) (Imw > a), 
(vi) q.(w) = [ dx e*“g(xz) (Imw > —»), 


(viz) g_(w) - | dx e*"g(x) (Imw < —d). 


*Received October 1, 1956; revised manuscript received January 21, 1957. 
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In terms of these quantities and by the Fourier inversion formula, (1) becomes 


| dw ¢ [f_(w) g:(w)] + [ dw ¢ [f .(w) 
: ; (2) 
— g,(w) — g.(w) — (AK(w) — 1)f,(w)] = 0, 
where tr > —1;c¢,d r <1,and uw < —c, —d. By a theorem due to Titchmarsh ({1]) 


p. 200), (Z yields directly for the transform of (1) 


[AK (w l}f.(w) = f.(w) + f_(w) — g.(w) — g-(w) — g.(w), (3) 
where in the region | Im w | < 1, the functions f_(w) — g_(w) and f,(w) — g.(w) — 
ga(w) — (\K(w) — 1) f,(w) are regular and tend to zero uniformly as | Rl w | — © in 


any interior strip. Since f,(w), g.(w) and K(w) are also regular and tend to zero uni- 


formly as| Rlw|— © inany strip in | Im w| < 1, the same properties hold for f,(w) — 
g.(w). 

3. Integral equations for f.(w) and f_(w). In order to evaluate the functions 
fa(w), f-(w) and f.(w) which determine the solution of (1), we shall need the following 


easily proved properties: 


(i) f,(w) and g,(w) are entire functions of w, 


(ii) f.(w) — g.(w) is regular for Im w > —1, 

(iii) f_(w) — g_(w) is regular for Im w < 1, 

- . X¥-( Ww) » 4 ae rc] » 

(iv) AK(w) — I] o(w) for —1< -a<a<1 ({I], p. 339 and [2], p. 51) 
xX w 


where 


(a x Ww is regular and free of zeros for Im Ww < a, 
(b) m. < | x_-(w)w”* | < M_ for Im w < a, where m_ and M_ are positive con- 
stants, 
(c) x.(w) is regular and free of zeros for Im w > —a, 
(d) m. < | x.(w) w'"* | < M, for Im w > —a, where m, and M, are positive 
constants, 
(e) o(w) = e(we(—w) and e(w) | ] (w — w,), where n is the number of roots 
w,; and —w, of \ K(w) — 1in| Imw| < a. 
Using the property (iv), we may rewrite (3) as 
x Ww) : a 
£ e(w)e w)f(w) = f.(w) + f_(w) — g.(w) — g_(w) — g.(w). (4) 


The present method depends upon converting (4) into two simultaneous integral equa- 
tions for f.(w) and f_(w). To obtain the first of these equations, we re-express (4) as 


e( —w) f.(w) — g.(w) — g,(w) f(w)—g (w) . 
_ L(w) — g-(w) — gw), , ae (4’) 
x. (w e(w)x_(w) e(w)x_(w) 


where from the properties (i), (iii) and (iv) the term on the left-hand side is regular 
> —a and the second term on the right-hand side is regular except at the 


for Im w > 
zeros of e(w) for Im w < a. We now consider a contour C, which connects the points 
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ia — © and ia + © by a curve which lies in the region —a < Im w < a@ and passes 
below the zeros of e(w) and a contour C_, which connects the points —ia — © and 
-ta + © by a curve lying below C, but above Im w = —a. Then we split the term 


fw) = g.(w) — g.(w) 
e(w)x_(w) 


’ 


and hence the whole of equation (4’), into two parts, one regular in the half-plane above 
C_, and the other below C, . Thus by Cauchy’s theorem 


f (pw) — q.(w) — g (w) ] f iat+@ a-ia-« 
! t= 5 —| + + | 
e(w)x_(w Qmi \Ic_. Je +e Jia-« 


v-1ta 


dw’ f.(w’) — g.(w’) — g,(w’) 


w —w e(w )x_(w’) 
where w lies within the closed contour. From the remark following (3) it follows that 
the third and fourth integrals on the right-hand side of this equation vanish. Thus 


, , / 
— w e(w )x_(w ) 


f Tos - g.(w q,(w) . 1 if / = dw’ f.(w’) a g.(w’) — g.(w’) 
u 


e(W)x (Ww 2r2 a JCa 


(5) 
= J(w, C_,) — I(w, C,). 


rom the definition of J(w, C_.), we see that this function is regular for Im w above 
C_., which results from the uniform convergence of the integral with respect to w in 
any region whose closure is in the region Im w above C_, . Similarly, J(w, C,.) is regular 
for Im w in the half plane below C, . Thus (4’) may be written 


i Be I(w,C.) = f.(w) {—o) _ I(w, C_.). (6) 

e(w)x _(w) , xX. (WwW) 
rom the preceding remarks, the right-hand side of this equation is regular for Im w 
above C_,, and the left-hand side is regular for Im w below C,, . Thus (6) defines a function 
regular in the whole w-plane and thus it defines an entire function. From properties (iv) 
b) and (d), 1/e(w)x_(w) and e(—w)/x,(w) are bounded for Im w below C, and above 
C_, , respectively. Further, from Theorems III and V of [2], f-(w) — g-(w) is bounded 
in every properly included half-plane below C, . Analogously, f,(w) is bounded in the 


half-plane above C_, . Finally, from the uniform convergence of I(w, C_,.) and I(w, C4) 
with respect to w, it follows that these quantities are bounded and approach zero as 
w|— © in the interior of their respective half-planes of regularity. Thus, the two sides 
of (6) are bounded and, hence, (6) defines a polynomial p(w) of degree at most zero. 
However, since f_(w) — g_(w), f,(w), I(w, C_.), and I(w, C,) vanish as | Rl w|— © 
in their respective regions of regularity, p(w) must be identically zero. Thus 
x.(w) 1 / dw’ f.(w’) — g.(w’) — q,(w’) (7) 
w) . —— ; 4 
e(—w) 2m7. w —w e(w )x (w’) 
for Im w above C_, and 
f_(w) — g_(w) l [ dw’ f.(w’) — g.(w’) — g,(w’) (8) 
. . ee ee eee , , a acerca 
e(w)x_(w) 2mri J, : w — w e(w )x_(W ) 


for Im w below C 


,* 
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In analogy with the calculation just described, if we start with (4), rewriting it as 


ic -, (w) — 
x-(w)e(w)e “*" f,(w) = x e "IT f.(w) — g.(w)] 
e(—w) : 
(w) ; e 
+4 +e" [f-(w) — g-(w) — g.(w)], (4””) 
eet | 


then split the second term on the right-hand side into two parts, one regular for Im w 


below C; and the other regular for Im w above C_; , where C_, connects —78 — © to 
—i8 + © by a contour passing above the zeros of «(—w) and contained in the region 
—8 Im w < £6, and Cs, connects 78 — to 78 + © and lies above C_, in the strip 
—8 < Imw < 8, then we have 
; ” ea 1 f¢ dw’ f_(w’) — g_(w’) — g,(w’) ia 
f.(w) = — - | ee xe (we (9) 
¢ x-(w)e(w) Qri J cg, VU — Ww Kc—wW } 
for Im w below Cz and 
; (w) if 
Lf -(w) — g.(-o)] A err = —H | 
e(—w 2m2 J, 
° : : , (10) 
dw’ f_(w’) — g_(w’) — g,(w’) a 
Ar a oe - ei a 100 Ol 
Co — @ ée(—W ) 


for Im w above C_s; . The presence of e **” and the specific arrangement of g,(W) in 
(4’) are required in order to bound the two sides of (4’’) according to the argument 
in particular, for this purpose, it is necessary to have the easily 


presented after (6): 
= as Im w — —o and|e ““[f,(w) — g.(w)]| = 


derived results that | e “*’f,(w) | = o(1) 


o(1) asImw-— o~, 

The important result is the pair of integral equations (8) and (10). The solutions of 
these equations, when combined with (7), determine f,(w) and finally f(x). Unfortunately, 
exact solutions can only be obtained for the case that K(w) is a rational fraction in w. 
However, approximate solutions can be readily obtained. 

4. Approximate solutions of (8) and (10). To solve (8) and (10) approximately 
we let the contour C, be shifted toward +7 and the contour C_, be shifted toward 
—io, The contours sweep across the singularities of the integrands of (8) and (10). 
The singularities are easily seen to be poles, except possibly for essential singularities 
arising from x.(w) and x_(w). If such essential singularities exist, they may be removed 
from the finite part of the plane by a simple device. Namely, since only values of k(x) 
in the interval —a < x < a enter into (1), k(x) may be chosen zero outside this interval 
as remarked in Sec. 2. As a consequence, K(w) becomes an entire function and, therefore, 
x.(w) and x_(w) have no essential singularities in the finite part of the w-plane. Thus 
if the contours are moved to infinity, the only singularities crossed will be poles. In 


this manner, (8) and (10) may be written 


f_(w) — (w) F (w’ I 
, i oN r(w) + > Res. (. J Au 7 = -) + R(w, a’), (11) 
; efens , am ee 


e(w)x_(u e(w )x_(w ) u 
, x.(w f_(w’)x.(w’) e**”’ 5 . E 
[f.(w) — g.(w)]- é = s(w) + ae, (es —+—— ] + S(w, 8’), (12) 
‘ e(—U e(—wW ) =< 
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1958 
where 
I dw’ q.(w’) + g.(w’) 
\W) = aT es i irae ames ict engage porte ; 
2riJc, W — we ew)x_(w) 
l [ dw’ [g_(w’) + g.(w’)]x.(w’) -scu? 
u eee oe =— SS a ~~ 9 
QriJc_,W —w e(—w’) 
; ; SF ax dw’ f.(w’) 
fi w, a J > 9 "e be 7 “6 x ager ee ’ 
271 Jiat-o W — Welw )x_(w ) 
; | [ aT! dw’ f_(w’)x.(w’) ~iaw’ 
S(w, 8’) = -== te 
2m J_i9 w—w e(—w ) 


In (11) and (12) the summations are over the residues at the poles of the integrands 
crossed by the contours as C,, is moved up to za’ and C_, is moved down to —7@’. Approxi- 
mate solutions of (8) and (10) are therefore provided by (11) and (12) if R(w, a’) and 
S(w, 8’) are neglected. The residues in (11) and (12) are obtained in this approximation 
in terms of the approximation itself by solving an obvious system of linear equations 
derived from (11) and (12). In the case that the contours C, and C_, sweep across 
multiple poles, the residues of (11) and (12) involve derivatives of f,(w) and f_(w). 
These derivatives may be evaluated by differentiating (11) and (12), and including 
the resultant equations into an extended system of linear equations for determining the 
residues. 

\ question may be raised with regard to the uniqueness and convergence of the 
residues obtained from the solution of the linear system derived from (11) and (12). 
[It has not been possible as yet to establish rigorously and generally this uniqueness 
and convergence. In a number of simple applications of the present method, in par- 
ticular to cases wherein K(w) is a rational fraction, the linear system has been found 
to have unique solutions; that is, the determinant of the system was non-singular except 
at eigenvalues. Moreover, for K(w) a rational fraction, convergence was automatic, 
since the determinant was of finite order and an exact solution was obtainable. However, 
since the determinant of the linear system involves only x,(w), x-(w) and e(w), and 
hence depends only upon K(w), and, moreover, is non-singular except at eigenvalues 
for some forms of K(w), it is a reasonable presumption that the determinant has no 
special degeneracy and that only exceptional forms for K(w) will lead to non-uniqueness 
or non-convergence. 

Finally an estimate of the accuracy of the approximation for f,(w) and f_(w) may 
be obtained by estimating the order of magnitude of the terms R(w, a’) and S(w, 8’). 


Using the easily proved results that | f.(w)e**” | ~ 0 as Im w— © and | f_(w) | ~ 0 
as Im w > — ©, we find that 


| R(w, a’) | < Ke _ 
and 
S(w, 8’) | < Me7*’ 
where K and J/ are constants. Thus the accuracy of the approximation improves expo- 


nentially with increasing a’ and 8’. 
The corresponding approximation for f{,(w) is determined by substituting (8) and 


(10) into (4), after neglecting R(w, a’) and S(w, 8’). 
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5. Approximation for f(x). From the Fourier inversion formula and from (4), we 


we find that, for 0 < x < a, f(x) is given by 
; DS a ae 
T(x) = on | dw’e ‘*" JalW) 
a i et é f .(w’)x.(w’)e" **"" 
ak. Ee 
On J—ss e(w )x_(w ) e(—w ) 
x .(w’) ( f_(w’) ) x.(w’) g.(w’) + g_(w’) + g,(w’) 13 
; “a ; - ie 7 i. , (1) 
€ u E(w )x_(W ) X-(W ) e(W )e(—wW ) 
where — © < 8 < o, For our purposes, we restrict 6 so that the contour in this equation 


lies below the zeros of e(w) and x_(w), below the essential singularities, if any, of x_(w) 
and below the singularities of g_(w), but above the zeros of e(—w) and poles of x,(w), 
above the singularities of g,(w), and above any essential singularities of x,(w). [For 
convenience it is assumed that w = 0 is not a root of «(w)]. Now in the first integral 
on the right-hand side of (13) we allow — 8 to increase to, say, a’; in the second integral 
we allow —8 to decrease to, say, —8’. This procedure leads to the following expression 


for f(x): 
‘patie f .(w’)x .(w’)e" **""’ 
aN ais ® Dae ae (. X ) i asl 
Mr) =2 7 Res. cia 7 : — ; + I(x, a’) 
_ e(w )x_(w ) \ e(—w 
ru , ° , 7 
é X.( Ww) j_(w ) : 
— 7 > Res. ; ( = 7 + J(z, B’) 
e(—wW ) e(w )x_(w ) 
Dy tails g.(w’) + g_(w’) + g,(w’) 
aaa | dw’ Sy Lae’, 7 x.(w’), (14) 
2r J-is e(w )e(—w )x_(w ) ; 
where 


: a pie Pam hiatal tr (eo' vl") — sau" 
I(x, a ; dw 7 a oe ee 7 e 
2r J; ' e(w )x_(w ) e(—w ) 


, citi , .(w’) f_(w’) 
J(z, B’) : | dw’e **" x mJ 2 a 
ae Seer ee e(—w ) \e(w )x_(w ) 


and 0 < x < a. The summations in (14) are over the residues at the poles crossed by 
the contours as — is increased to a’ in the first integral of (13) and is decreased to —p’ 
in the second integral. 

The final approximation to f(x) is obtained by neglecting J(z, a’) and J(z, 6’) in 
(14) and by using the expressions (11) and (12), neglecting R(w, a’) and S(w, 6’), for 
f.(w) and f.(w). A straightforward calculation shows that the error resulting from 
neglect of I(x, a’) and J(z, 6’) is 


| T(x, a’) a mae 
and 
| J(z, B’) | < Be*’*, 


where A and B are constants. Thus for large a’ and @’, the error is small except in the 
neighborhood of x = 0 and x = a, which for increasing a becomes less and less important. 
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6. Eigenvalues of (1). In the case that g(x) = O in (1), we have an eigenvalue 
problem. Then (11) and (12) become a system of linear homogeneous equations for 
the determination of the values of f,(w) and f_(w) at the poles of the integrands of 
(8) and (10). The values of \ admitting of solutions of (11) and (12) under these cir- 
cumstances will clearly then constitute approximate eigenvalues for (1). 

7. Example: 


{(@) =e "'" +A | dy e"'*~*' f(y). 


“0 
To illustrate the application of the present method, the above simple integral equation 
will be treated. This equation admits of an exact solution since the transform of the 
kernel is a rational fraction and, consequently, the series expansions in (11), (12) and 


(14) terminate. 
rom the definitions (i)—(vii) of Sec. 2 and property (iv) of Sec. 3, we find: 


] 
q_(w) : 
‘ a + iw : 
i] (w) ~~ 
= oP 
=e a 
g.\w) = a 
a— Ww 
2—-l—w 
, tol 
A\AK(w) — 1 = ; , 
w+ ] 


ew) = w — (2A— 1)”, 
x. (w) w+, 
l 


x (Ww) 
w=— 3 


Substitution of these quantities into (11) and (12) and evaluation of the residues leads 


to the result 


| (QQ — 1)'* -—i ie ais ” l 
w + < 5 Qi — 7- .({2x — 1]'”?)e? ——., 
J ; at ww a-=7> 1(2r — 1)’ ™ + (( r 1) vl) (| r 1] } wo-— 
f.(w) = ——— + 4 [-(24 — 1)? + 3G] 
‘ a iw \ 
| l 4 1 — exp {—alfa + 7(2d 7D) 
a — 1(2\ — 1)'” a +22 — 1)” 


-exp [ia(2\ — 1)'7] + [(2\ — 1)? — a] f_(—[2a — 1]'”’) 


-exp [ia(2\ — 1)’ | ee’ 
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Forming from these equations the linear system for evaluating f,((2\ — 1]'””) and 
f-(—[2\ — 1]'”’) and carrying out the calculation gives 
ita «ifs a LS 2 te { tO exp fia(2a — 1)"?] 

Mle , Bias R-—1 kha —-p7—-1 } 


l—a 


iT — D4 1' a 


f_(—[2n ed 7 : 1, a — D+! » [ta(2X — 1)*”*] 
A-ha — Ts ne eee ex -- 
JAF Di’ +S — 5 ae —- 1 ) 
a+ 22 — 1)’ 2 a(2r - 1)' 7] = —_# es 
3 Lb aa _— = a ‘ aa 
‘i - 5751 sb ie (2. — 1)? — 1 
where the determinant, D, of the linear system is given by 
p = &P [—ia( 2) - — ae _ exp [ta(2\ — o. 
(i[2x — 1]'” )* (i[(2\ — 1] — 1)° 
The values of \ for which D = 0 are the eigenvalues. For the present case, \ = 1/2 
is the smallest eigenvalue. 
For 0 < z < a, evaluation of the residues in (14) yields 
f(x) = — — f.({2. — 1]'”*) exp [—7x(2A — 1)'”7] 
* , , (2r ] « 7 +\ Lo? 4 I ‘ - a 
rr F , +1 2a 1 : 1/27 
~ 7a §f-(-—[2A — 1] °°) - = ? exp [iz(2\ — 1)°”*] 
(2) — 1) \J-\—! , Jta- yom iz. 
ai — 1 - 
aL é 


If D # the above expressions for f,({2\ — 1]'””) and f_(—[2A — 1]'”) may be sub- 
ales into ‘this equation for f(x), which after rearranging gives the final result 


J 1 + @ + 46 1/2 ' i= a = -e«\ 
Ke) = Vix — 1 — 1 oP OA — I aK Fi J 
exp [=12(2A — 1) Ty l+a ; 1/2 
~ exp [—ia(2d — 1) 
Die +2x—1) + @r—-1*+1 Pl-* 
= ae ,a\, exp [ix(2 — tt lin went — 
a—-jp7—-1i" | Be+a-D “etn -—! 


This simple example is typical of the application of the present method to equations 
for which K(w) is a rational fraction. For other functions K(w), the procedure is com- 
pletely analogous, but has the same difficulty as with the usual Wiener-Hopf method 
that the functions x.(w) and x_(w) are in general not evaluable analytically but must 
be treated either numerically or approximately. 

8. Example: 

" 0a 


_—. “ ” eX a | 21¢ 
Rt es (2m)? i dy exp [—3(z y) \f(y). 
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For the case of the Gaussian kernel, the present method has been applied to the deter- 
mination of the lowest eigenvalue. Previous consideration [3] of this problem has led 
to the calculation of upper and lower bounds for this eigenvalue. The present calculation, 
while illustrating the above method and indicating its accuracy and limitations, provides 
an extension of the earlier work on this problem by obtaining an accurate estimate of 
the eigenvalue in the region where e bounding values differ significantly. 

Taking the transform of k(x) = 1/(2x)'” exp (—2x’/2), we find 


\K(w) — 1 = \ exp (—43w”’) — 1. 
The roots of this expression occur at 
w, = + [2 Ind + 4nni]'”, 
where n is a positive or negative integer or zero and it is known [3] that \ > 1. For the 
approximation described below, we limit \ so that only the two roots with n = 0 lie 
within the strip | Im w! < 1. This requires that \ S 20. With this restriction, we may 
write 
a tile x_(w) , 172 7 172 
[AK(w) — 1] = = [-—w — (2 Ind)" ][w — (2 Ind)"”’)]. 
x.(w) ; 
From ([1], p. 339) it follows that 


exp [— 7_(w)] 








x-(w) = eee oe 
ind 
, 1 
x.(w) = may er : 
where 
E dw’ an “ 1 
iat Qri iz £2 = nfl — d exp (—}u Ny wtih 


and the contour lies in the strip 0 < Imw < 1. 
From (11) and (12), it is found that 





p(w) f(2innf) 
e(w)x_(w)  x-({2 In A] *7 (2 Ind)? — w’ 
. x.(w we _ f$-(—[2 In d]*)x.(—[2 In sai ) exp [ia(2 In \)'”7] 
"? ¢ (w) - “0 In A)'? + w ; 


where residues at the roots w, , n # 0, of \K(w) — 1 have been neglected. These terms, 
as argued at the end of Sec. 4, are of the order of exp[—a | Im (2 Ind + 4ni)'”*). For 
a ~ 1, it may be shown from [3] that A ~ 2.5 so that | Im (2 In \ + 4mi)'” | ~ 2.3. 
The error is then of order e~*** ~ 1/10 and decreases rapidly as a increases. 

The linear system for determining f.({2 In \]'”) and f_(—[2 In \]'”*) may now be 


written 





f-(—[2 In a} .({2 In A]? -0 
x-(—[2 In A]? Le =a , 
j [2 In d]'””)x.(—[2 In d]'”) exp [ta(2 In d)'”?] 


+ f.({2 In AJ'”)x,({2 In AJ”) exp [—7a(2 In d)'”?] = 0. 





The eigenvalues are determined by the 
the vanishing of the determinant of the 


_x+([2 In ] 


5 (— fa in np) OP [| 


From the definitions of X+(w) and xX (Ww) 
(2 In A)? + i)" 
; (2In dA)” —7 


As a further simplification the expression for r 


exp [22a 2 In\X 


—ih 


r_({2 In A] *~) 


where 


Py} 


njil —A 


\ 


J({2 In A] 


and the integration is along the real axis, P 


above expression becomes 
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' 7 a, 
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solubility of this linear system which leads to 
system, namely, 


x.(—[2 In A]’”*) 


2 0. 
x_({2 In A]”) 


exp [ia(2 In ))' 


we may re-express this equation as 


exp [—2{7_([2 In A]'’*) — 7_(—[2 In A]’”*)}]. 


({(2 In A]’’~) may be written 


‘ 1/2 
-* t 3 J((2 In A]’*), 


5 


a +1 i dx 


: Fe iS ' 
exp (9901 -2\IndXJ x — (2 Ina 


standing for “‘principal value.’’ Then the 





oe (2 In d)'”? + 7\’ —4.J({2 In d]'”’) 
exp [7a(2 In \)**] 7s] «exp ; —. 
(2 In A)" — 2/ 2m 
A simple argument shows that this yields 
9 l I l ‘ i 
2tan = 73 J({2 In A] + n3 
(2 Ind vi 
a : ° i72 
(2 In A) 
where n is a positive or negative integer or zero. It may be shown that as \— 1, 
J ({2 In d]'’*) — 0 as 2 In X. Thus, since for the lowest eigenvalue a > © as \ — I, it is 


easily seen that n 0. The larger eigenvalues correspond to n > 0. 
: Ws Bis 
An evaluation of J({2 In \]’””) was carried out numerically so as to exhibit the values 


of \ as a function of a. The results are contained in the following table. 


( r I; Ty 
6.03 1.10 1.00 1.15 
4 03 1.20 1.05 t Zo 
2.65 1.40 2 1.43 
2.05 1.60 1.44 1.62 
1.70 1.80 1.65 1.81 
1.46 2.00 1.87 2.01 
1.10 2.50 2.39 2.50 

891 3.00 2.91 3.00 

747 3.50 3.43 3.51 

641 +.00 3.98 +. 04 

559 $50 $.54 4.60 
.492 5.00 5.14 5.19 

390 6.00 6.47 6.51 

253 8.00 9.93 9.93 
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The quantities [, and I, are 


9 na/2 =% 
I, 9,172 | exp (— 32°) ix} ; 


| . [ exp (— 42°) dx ~— ete 1 — exp (—40°)} 
* sa ites a(2m)'’* 2 , 


1 


D2) 
“7 } J0 


and correspond to the lower and upper bounds of the eigenvalue deduced in [3]. It is 
observed that to three figures the eigenvalue obtained above lies within the bounds 
I, and I, for a 2 0.6 which is consistent with the earlier error estimate. For smaller 
values of a the present approximation begins to break down as expected. 
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BOOK REVIEWS 


(Continued from p. 20) 


Gas dynamics. By Klaus Oswatitsch. English version by Gustav Kuerti. Academic 
Press, Inc., New York, 1956. xv + 610 pp. $12.00. 


The rapid development of supersonic aerodynamics during the past twenty years has created a 
need for authoritative text books. Since most people working in this field have been engaged very fully 
on highly specialized problems and since so many and varied contributions have been made to the 
subject it was inevitable that the first leading works of reference should have been prepared by teams 
of writers rather than by individuals. In the circumstances it was unavoidable that these works, excellent 
though they are, should suffer from a non-uniformity in style and, owing to difficulties in liaison, should 
have been partially outdated at the time of their publication. 


Dr. Oswatitsch has therefore performed a great service in preparing the present volume single 
handed. He provides us with the opportunity of covering the complete range of Gas D mies in a 


style which is even, clear and easily followed throughout. The work is of considerable scope and covers 
almost every important aspect of the subject. While the treatment is thorough and rigorous the emphasis 


many interesting phenomena occurring in Gas Dynamics, without introducing 


is‘on explaining the 
complications due to Mathematical details or laboring particular practical applications. As the author 


himself states this is not an Engineering handbook, nor, it might be added, is it a Mathemat treatise. 
: £ 
The aim is clearly to give n understanding of the phvsical character of the subject and this ertainly 
. Ss , . 


accomplished. 
The book follows a logical order. Thermodynamic principles and one dimensional models of gas 
j £ 
flows precede chapters on general equations and theorems. Three chapters follow on steady subsonic, 


st since 


steady supersonic and steady transonic flow respectively. This last chapter is of especial int 
it not only contains a much tidier account of recent work in the United States than has appeared else- 
where but also present the author’s own approach to this subject with a clarity which the reviewer 
missed in earlier versions. After a fairly complete chapter on linearized theory the book concludes with 
a chapter on viscous effects and a chapter on experimental methods. 

Dr. Kuerti is to be congratulated on his English translation of the book. Unlike some earlier English 
versions of German text books, which were often harder to follow than in the original, the present work 
reads as though it has been prepared in English. This is evidently the result of the translator taking the 
trouble to understand the whole contents and of his close collaboration with the author. 


The book is very well set out; it is liberally illustrated with diagrams and remarkabl) ir photo- 
graphs and the list of references is excellent both in its completeness and in being up to da It would 


serve as a first rate text book in a graduate course but is also readily understood by the student working 


without supervision. It is strongly recommended to all those engaged in work on high speed aerodynamics. 


M. Hotr 


Abacs or nomograms. An introduction to their theory and construction, illustrated by 
examples from engineering and physics. By A. Giet. Translated and revised by J. W. 
Head and H. D. Phippen. Iliffe & Sons, Ltd., London; Philosophical Library, New 
York, 1956. $10.00 


This book was originally published by Dunod in 1954 under the title ‘““Abaques ou nomogrammes”’. 
The scope is indicated by the following chapter headings: Relations between two variables—Cartesian 
abacs—Alignment charts—Alignment charts not based on parallel coordinates—Relations between n 
variables. The mathematical treatment is elementary, avoiding even the use of determinants, but all 
important types of nomogram are thoroughly discussed and illustrated by numerous examples. 


W. PRAGER 


(Continued on p. 46) 
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ON THE HEAT TRANSFER TO CONSTANT-PROPERTY LAMINAR 
BOUNDARY LAYER WITH POWER-FUNCTION FREE-STREAM 
VELOCITY AND WALL-TEMPERATURE DISTRIBUTIONS* 


BY 


ISAO IMAI** 
Institute for Fluid Dynamics and Applied Mathematics, University of Maryland 


Summary. The heat transfer to constant-property laminar boundary layer with 
power-function variations of free stream velocity (u, = cx”) and of temperature difference 
between wall and free stream (7, — 7, = bz") is studied by means of an improved 
version of the WKB method developed by the author. It is found that the local heat- 
transfer coefficient h can be approximately given in the form 


hal = en | an ta + nce — oytean — |, 
where 8 = 2m/(m + 1), a is the non-dimensional velocity gradient at the wall (usually 
expressed as a = f’’(0)), o is the Prandtl number, k is the thermal conductivity, and » 
is the kinematic viscosity. 

As examples of application of the formula two cases of special interest, i.e. the Blasius 
flow and stagnation flow are considered for several values of n, and the results are com- 
pared with the purely numerical computations-by Levy and by Chapman and Rubesin, 
as well as with the analytical approximate formulas obtained by Schuh, by Fettis and 
by Lighthill. It appears that the above formula is the most accurate among those hitherto 
proposed. 

1. Introduction. The calculation of the heat transfer to the laminar boundary 
layer with arbitrary distributions of main flow velocity and wall temperature is very 
difficult even for the constant-property flow, because the problem is governed by a 
system of non-linear partial differential equations. But, as was shown by Fage and 
Falkner [1] the system of equations can be reduced to that of ordinary differential equa- 
tions if both the main-stream velocity and the wall temperature are assumed to vary 
as power functions of distance from the start of the boundary layer, namely u, = cz”, 
7, — T, = bx”. In fact, the momentum equation has been exactly solved by Hartree [2], 
and it remains only to solve the energy equation with variable coefficients given in 
terms of the Hartree solution. The exact treatment of the latter problem is naturally 
only possible by numerical means just as in the case of the momentum equation. Thus, 
Chapman and Rubesin [3] obtained the numerical solution for several values of n for 
the case of the Blasius flow (m = 0) by use of the Runge-Kutta integration method, 
assuming the Prandtl number o to be 0.72. Later, Levy [4] carried out extensive numer- 
ical computations for a number of combinations of m, n, and o by reducing the differential 
equation to finite difference form and employing the IBM machine. From his numerical 
results, Levy proposed an empirical formula for the local heat-transfer coefficient. 
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In view of the considerable amount of labor involved in the numerical computations 
with three parameters m, n, and o, Schuh [5] attempted to obtain the asymptotic solution 
of the energy equation when either n or a is large so that the thermal boundary layer is 
very thin compared with the velocity boundary layer. Making further assumption that 
a/n — 0, he succeeded in obtaining the solution in a closed form using the Bessel function 
of order 1/3 and hence a very compact formula for the local heat-transfer coefficient. 
Also he showed that the formula gives numerical results in good agreement with Levy’s 
numerical computations. Dropping the second assumption ¢/n — 0, Punnis [6] and 
Fettis [7] succeeded in finding the solution in a closed form in terms of the confluent 
hypergeometric functions. Especially the latter author derived an explicit formula for 
the local heat-transfer coefficient, but he did not make any numerical discussion of the 
formula obtained. 

In a paper which probably escaped the attention of the above-mentioned three 
authors, Lighthill [8] had developed a method of approximate calculation of heat transfer 
to the laminar boundary layer with arbitrary distributions of wall temperature and 
main-stream velocity by making the simplifying assumption that the velocity profile 
in the boundary layer can be approximated by its tangent at the wall. As an example, 
the method was applied to the Blasius flow and the result was found to be in good agree- 
ment with Chapman-Rubesin’s numerical calculation. 

On the same assumption as made by Lighthill, Tamaki [9] independently developed 
an approximate method of calculating thermal boundary layers, especially for the case 


of the Blasius flow; his formula for the local heat-transfer coefficient is essentially the 
same as Lighthill’s for this special case, but he treats a number of interesting cases 
using his formula and also extends the method so as to cover even the turbulent poundary 
layer. Further, he carried out numerical computations for our problem for the cases of 
B = 2m/(m + 1) 1, 0.5, 0, —0.14, —0.1988 and the range of n from —1.0 to 1.0 with 
o = 0.7. 


In the present paper, the basic energy equation is solved by use of the 
version of the WKB method proposed by the author 110, 11, 12] which has been shown 
by experience to give accurate results for various problems governed by the linear 
ordinary differential equation of the second order. The local heat-transfer coefficient is 


1 
} 


given in terms of algebraic functions of the three parameters ™, n, and oc. Although the 
full expression obtained is somewhat complicated, it can be approximated by a rather 
simple expression. Applying the formulas, both full and simplified, to the cases of the 
Blasius flow (m = 0) and the stagnation flow (m = 1), it is found that they are in very 
good agreement with Levy’s numerical computations. 

2. Basic equations. For the laminar boundary layer on a body placed in an incom- 
pressible constant-property fluid, the continuity, momentum, and energy equations 


can be written as 


+) a 
—+—=40, (2.1) 
OX OY 
Ju yu dp Vu : 
Au +o) = 24 2, (2.2) 
Ox OY dx Oy 
oT OT | oT _ 
pU, ( —-+0 —) = k—;, (2.3) 
Ox OY oY 
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where x and y are measured parallel and perpendicular to the surface of the body, u and v 
are velocity components, p is the pressure, 7' is the temperature, p is the density, C, is 
the specific heat at constant pressure, k is the thermal conductivity, and u is the viscosity. 
The boundary conditions are 
T = T,(2), u=v=0 for y=0 (2.4) 
T=fT,, u—u(z) as yuo (2.5) 
where 7’,() is the surface temperature and 7’, and u,(x) are respectively the temperature 
and velocity of the main stream. 
In the following, it will be assumed that 


T, —T, = b2", (2.6) 
wu, = cx”, (2.7) 


where }, c, m, and » are constants. In this case it is known that the system of partial 
differential equations (2.1), (2.2), and (2.3) can be reduced to a pair of ordinary differ- 
ential equations. Thus, introducing the variable 


1 = (5 1) y (2.8) 


2 vx 
and assuming that the stream function y defined by 
u = Oy/dy, v = —dy/dx (2.9) 


and the temperature 7 can be expressed in the forms 


9 1/2 
y = (- + = vu) f(n), (2.10) 
= — A = 6(n), (2.11) 

Eqs. (2.1), (2.2), and (2.3) are reduced to 
"+ ff" + al -— f”) =90, (2.12) 
6’ + of 0’ — on(2 — B)f'6 = 0, (2.13) 
where 8 = 2m/(m + 1) ando = uC,/k is the Prandtl number. The boundary conditions 
f=f' =90, 6=1 for »=0 (2.14) 

nd 

f'>1, 090 a q-> >. (2.15) 


Equation (2.12) is well known as the Falkner-Skan equation and has been solved 
by Hartree [2] for several values of 6. Hence Eq. (2.13) can be regarded as a linear 


differential equation of the second order with known variable coefficients. 
3. The improved WKB method. Tor the integration of Eq. (2.13) it is convenient 
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to make use of a refined version of the WKB method as developed by the present author 
[10, 11, 12]. Thus, consider a second-order linear differential equation 


d’® 





—=+a°P(x)b =0, am, (3.1) 
dz 
where a is a parameter and P(x) has the form | 
P(x) = a,x(1 + dix + bow? + ---), a, ~ 0. (3.2) 
Then, it can be shown that the solution of Eq. (3.1) can be expressed as 
& = P~'* exp (Q/4a’){C, exp [r(az — Q,/2a)] 
+ C, exp [—7(az — Q,/2a)|} + O(a") for rx #0 (3.3) 
and 
ma\'’? 41/6 5/6 , — >) 
=({—) Pe °e{ 0,0 H,)3(e) Ca ite) 
& (72) (C0 H 2) + Ci" H2u)} _ 
+ O(a *) for «= O(a”), 
where 
s= | Paz, (3.5) 
4 went. «we 
4 dz dx 
5 | se (3.6) 
ala 5 + ie he + he” +--+), 
(2. “¢ * ; 
- \ - 
a ] y J 
Q | (0+ 3) dz + (3.7) | 
‘ iii 36 2° 36 2 
w = KE+ BAK (3.8) 
¢ = (é - LK : ae (3.9) 
: + Xk 3.10) 
K" a +X, 3.11) 
and the constants XA, A, , A, are given in terms of a, , b, , «++ as follows: 
3 /3 pare 
r 35 (3 a, ] (5b, — 3b), 
4 ] . = 
A ———— (25b om 35b,b, + | th ¥ 
io a, 3..o) 
27 (3 -* , 
= ((1 95 — 4250p’ 
26950 (5 a)  (6125b, _— 


— 9800b,b; + 12080b7b, — 3624b)). 


In order to consider the behavior of ® as x — 0, it is convenient to use J.1/;() 
instead of H{);”(w). Thus, we have 
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HR = gal ia — BS), (3.13) 
His = si r* *(a' a a = "Sand (3.14) 


and 


POT y;2(w) 


3 —-1/6 2 1/3 wil 1 2 " 3 1/3 
= (3a) Reape [i Be 2 Ga) eft], Gap 
wenitilie fe 3 -1/6 gs | {i 
2° ¢J-1,:(~) = (3 a,) T(2/3) K 1 — 5 b, 


eT ee ee 
+ (3a +2)Ga.) lt adie | 


4. Solution of the energy equation. In order to apply the above formulas, we have 
first to bring the basic equation (2.13) into the standard form (3.1). This can be done 


by putting 
1 ” 
6 = ep| —he [ san]-0 (4.1) 
“a J0 


6” — P(e = 0, (4.2) 


(3.16) 


so that Eq. (2.13) becomes 


where 


P(n) = j of? + o! + n(2 — aby (4.3) 


Now, Eq. (4.3) should be expressed in a form similar to Eq. (3.2). For this purpose, 
differentiating Eq. (2.12) with respect to 7 three times, putting » = 0 each time, and 
using the boundary condition (2.14), we have 
f =f’ =0, f{" =a, f’'’ = —8, 

y ‘Vv 2 -V (4.4 
j** =0, f* = (2 — Ba’, f*' = 2(2 — 38)Ba (4.4) 


for » = 0. Here we have put f’’(0) = a, whose values have been calculated by Hartree 


[2] for several values of 8. Hence 
° ] 2 1 3 ] 9 2 5 1 ‘ 6 ~ 
j= Zan’ — 5 Bn? + =, (28 — Na’n + 5, 22 — 38)Ban' + --- (4.5) 
Substituting Eq. (4.5) into Eq. (4.3) we have 


8 28 — 1 rat ' 
Poa) = Anf - Sent tet ate (4.6) 


- 
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where 


A = ag{4 ™ n(2 — B)}. (4.7) 


Comparison of Egs. (8.1) and (3.2) with Eqs. (4.2) and (4.6) shows that 


a e** e (4.8 
B 
a { bh . = b = 0, 
2a 
» ol * (4.9) 
} oa (1 — 2B)a cn (2 — 38) 8 og ap | 
164A 24 ; es 60 24A} 


The formulas (3.3) and (3.4) have originally been obtained on the assumption that the 
parameter a is very large. At first sight Eq. (4.8) might seem to conflict with this assump- 
tion. Since, however, the condition for the validity of the formulas is essentially that 
| az | should be large, the condition is readily seen to be satisfied if A is large. This latter 
condition is satisfied if either the Prandtl number o or the exponent n is large as was 
assumed by Schuh [5]. Furthermore, it should be mentioned that the formulas have 
been found to give accurate results in various cases of application even when a is not 
very large. 
Let us assume A > 0. Then P(n) > 0 for 7 > O Consequently 


z= | P’? dn > 0 (4.10) 


70 


so that 


102 =—2 < v. 
Therefore, by the boundary condition (2.15), which can also be expressed as 


6—0 as zo 


we have to take C, = 0 in Eq. (3.3). Thus 
6 = C,P™* exp [—(z + 30, + 1Q)] for » #0 (4.11) 
and 
mi\'’? 
6 = 0,(% Pee eH) for 4 = 0. (4.12) 


The constant C, can be determined by the boundary condition (2.14) as 


C, = —(2/mr)'’*[H(0)]™’. 1.13) 
Hence Eq. (4.12) becomes 
6 = H(n)/H(0) for 7» = 0, 1.14) 
where 
Hq) = Poe 2 gt" Hw). (4.15) 


Now the heat flow rate at the wall is given by 


‘ (22) | 4 ae (“.) (2) . (4.16) 


\ OY 
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Therefore, if the local heat-transfer coefficient h is defined by 


qo 
Awe ——., 
re - is 


we have 


ha/k (de 
—L—5 = —(2 — fp)" ( s) 
(U,2/v) dn 


7=0 


But, by Eqs. (4.1) and (4.14) 
(42) (22) _ HO) 
dn/ ,~0 dn/,-0 HO) 
Now, combining Eqs. (3.6), (3.8), and (3.9), we have 


2°61 — »,)A]"'”* 2-*/*1(2/3) will } 
e Dae ME wm b> 
HO 378 1(2/3) L+ “a73) 


(0) 2\"" va L — ef(9/8)A? + $d] — do)" 
H"(0 l As (: ) (a — ajay” 1 — e[(9/8)A" + 5) Kd = Xo) 
H(0) 9D \3 c+ XI — Ad) 
where 


oor 4 fa) 
> - 1.04682 
: (2/3) 1.046 


and the constants b, , A, Ao , A, are given by Eqs. (3.12) and (4.9) as follows 


2 \ -2/3 + 
Y (5 : 
| a ' { 
\ — A '(25b, 4 145;), 
io 
| 
97 5) \ 1 
r -—— (3 A} (6125b, — 9800b,b, — 3624b'), | 
26950 \2 J 
vhere 
8 } (j — 23)a ca 
<a 24 1GA ’ 
(2 — 38)8 oO a3 : 
b, = . ’ A = {i n(2 — B)\oa. | 
60 24A 2 + cous 


If it is assumed that A is very large while b, , bs , and b, are O(1) (this is true 
very large so that A = O(n)), Eqs. (4.23) give 
A = Q(A~? a No _ 0(A~'), Ay _ o(A~* 3 


Hence Eq. (4.21) can be approximated by 


H’(0) B (1)" T(2 3) 1/3 -1/3 
it A I ES a = An”), 
H(0) 10a 3 T'(4/3) ts ) 


(4.17) 


(4.18) 


(4.19) 


(4.20) 


(4.21) 


(4.23) 


(4.24) 


if n is 


(4.25) 
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5. Comparison with previous investigations. The temperature distribution in the 
boundary layer is given by Eqs. (4.1), (4.11), (4.13), and (4.14) in the form 


A(n) = _(?) . - exp —\3o [ fdyn+2+34Q, + 10) for ~ 0 (5.1) 
AN \0r H(0 I 2 ! J an Te 21 T 4 n 7 io. 
and 
H(n) , 1 os P = n 9 
H(0) exp | 50 | J an| for ” = 0. (5.2) 


In the following we shall restrict ourselves to the heat-transfer coefficient h, which is 
given by 
he/ke | 1v2 H’(0) 


—(2 — £8) (5.3) 
H(0) bo 


(u,x/v)? ~ 
Here H’(0)/H(0) is given by Eq. (4.21) or approximately by Eq. (4.25). 

First let us compare our result with the results of previous investigators. 

Schuh [5] considers that the thickness of the thermal boundary layer becomes small 
compared to the velocity boundary layer if either n or o becomes large, and that it 
would then be sufficient to replace the velocity profile by its tangent at the wall, since, 
for calculating the temperature field, only that part of the velocity profile is important 
which lies within the thermal boundary layer. Thus he puts - 


df/dn = an (5.4) 
so that Eq. (2.13) reduces to 
6’’ + tacn 6’ — ona(2 — B)n6 = O (5.5) 
which, on putting 
¢ nn (5.6) 
becomes 
os + - rs — ga(2 — B)c@ = 0. (5.7) 
Further he assumes that 
a/n — OV 5.8 
so that Eq. (5.7) simplifies itself to 
d*6/dé* — oa(2 — B)¢6 = O. 5.9) 


This can be integrated in terms of the Bessel function of order 1/3. The final result for 


the heat-transfer coefficient is 


hz/k T'(2/3) (aan)' 


—"- == ~<a" 5.10) 
(u,x/v)”* ~~ 3? T(4/3) ( ». 10) 


| 


bo 


= 8) v 


Next he tried to obtain an asymptotic solution of Eq. (5.7) which is valid for ¢ — 
with no restriction on n. For this purpose, he introduced a new variable ~ defined by 


f= (ca)? 9 (5.11) 


Ss 
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to write Eq. (5.7) in the form 


d’o/de? + 42° dé/dt — &0 = 0 (5.12) 


with 
e = n(2 — 8). (5.13) 


However, he could not find any solution of Eq. (5.12) in a closed form and made no 
attempt to solve it numerically. As examples of application of the formula (5.10), he 
calculated the values of (ha/k)(u,x/v)~‘’” for the Blasius flow (6 = 0) and for the stag- 
nation flow (8 = 1) for several values of n and compared the results with the purely 
numerical calculations by Chapman and Rubesin [3] and by Levy [4]. 

Later, Punnis [6] and Fettis [7] independently showed that Eq. (5.12) can be inte- 
grated in a closed form using confluent hypergeometric functions. Furthermore, the 
latter gave an explicit expression for the heat-transfer coefficient at the wall. His result 
for the temperature distribution function 6 is expressed as 


r(1 + 2/3) , +3/12\ 147 +3 /a ~ 
0 T1733) ‘exp(— 8°") W_ 1 + 20/3) 1 0(€'/6), (5.14) 


denotes the Whittaker function and & is defined by Eq. (5.11). The 


where W, am 
heat-transfer coefficient is given by the formula 
he/k___ (oa/6)"* P(A + 2¢/3)1(-1/3) dil 
u,a/v)'”” (2 — g)'” (2/3 + 2/3)T(1/3) cs 


Now it should be mentioned that Lighthill [8] had previously obtained a general 
formula which gives the heat transfer at the wall in terms of the skin friction for arbitrary 
main-stream velocity and wall-temperature distributions, based on the simplifying 
assumption adopted by Fage and Falkner [1]: 

u = (1/p)r(x)y, (5.16) 
where r(x) is the skin friction so that r(x) = u(du/dy),.o . This assumption is of course 
the same as Eq. (5.4). Lighthill’s formula can be written in our notation as follows. 


ate) = Meee) Hay Uf wor a) 


Lighthill himself applied this formula to the case of the Blasius flow in order to test its 


—1/3 


dT ,(2,). (5.17) 


accuracy. 
Now the formula (5.17) will be applied to our problem defined by Eqs. (2.6) and (2.7). 
First we have, by Eqs. (2.8), (2.9), and (2.10), 


9 ott’ ees . 
(x) = u( 2) = pu, (t+ Hs) , (5.18) 


VI 


Introducing Eqs. (5.18), (2.6), and (2.7) into Eq. (5.17) and making some calculation, 
we obtain 


TUL + 2/3) M(-1 3) (5 19) 
(2 — pg)’ 1(2/3 + 2/3) T(1/3) sais 


_ _ (oa/6)"* 


q (x) ; 
T, — 7;)k(u,/vx)'” 
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which is nothing other than the formula (5.15). Therefore, it may be stated that Fettis’ 
result is included in Lighthill’s general formula. This is perhaps what should have been 
expected, because both results have been obtained under the same simplifying assumption 
about the velocity profile. Further it may be remarked that some idea of the accuracy 
of Lighthill’s formula (5.17) may be obtained by comparison of the result (5.19) with 
our and other calculations. 

Extensive numerical computations have been performed by Levy [4] by reducing 
Eq. (2.13) to finite difference form and employing an electronic IBM machine; the results 
cover a range of n from —2.5 to 4.0 and o from 0.70 to 20.0 for four values of m/(4, 1, 0, 
—0.0904). For the range of the calculations, it was found that the local heat-transfer 
coefficient can, with the exception of large negative n values, be expressed within 5 


per cent as 
_ hz k ; . 


WE Bim, no", 5.20) 
(U,2/V) 


where the function B(m, n) can be approximated by 


9 \ 0-104 9 \ 0°374+0°12m/(m+1 

f = my am sles an os 

Bim, 7 0.57(—2"_ +. 0.205) (— + 1) = 91) 
m -+- ] / m-+ | 


and the exponent y of the Prandtl number varies as shown in the following table: 
3 1.6 1.0 0 —0.199 
4 0.367 0.355 0.327 0.254 


6. Blasius flow. For the case of the Blasius flow we have 6 = 0, so that Eqs. (4.24) 


and (4.23) give 





0 b=—-=, w&=0 
16A 24 ; 
] oa | a 
r 0 he = —— a — 1 =0 
12 A’ 18 A 
Hence Eqs. (5.3) and (4.21) become 
ha/k 1 (2/3) gis 1 1 o’a*\'”" : 
5 = oa mA’ ( —-—~4—455 (6.1) 
(u,x2/1 2°° 3°°°T(4/3) IS A 12 A 
which can be approximated by 
ha/k 1 1(2/8) az 
777178 = Sift arian » (6.2) 
(u,x/v) 2° 3° T (4/3) 
where 
A = (4 + 2n)oa, a = 0.4696. (6.3 


In this case, Schuh’s formula (5.10) becomes 


ha/k_ _ Y(2/3)(noa)"*” (6.4) 
(ue)? = 2°°3?* 7 (4/3) —- 
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and Fettis’ formula (5.15) becomes 


_ha/k tag-1/3 T'(2/3)V. + 4n/3) 


172 = 1/3 6.5 
(U2 /y) "4" J 1(4 /3) r(2 /2 4 4n/3) (oa) ( ).0) 





which is essentially the same as the formula given by Lighthill. Table 1 shows the values 
of (ha/k)(u,2/v)~'’* for several values of n as calculated by the various formulas. In 
Table 1 the first column shows the values calculated by Chapman and Rubesin by the 
use of the Runge-Kutta integration method. The Levy values were computed by the 
present author from Levy’s values for « = 0.7 by using Levy’s empirical formula (5.20). 


TABLE 1 


Local heat-transfer coefficients for Blasius flow (8 = 0). Prandtl number = 0.72. 


Chapman- Lighthill 
n Rubesin Levy Fettis Schuh Imai Imai* 
0 0.296 0.2901 0.304 0 0.3011 0.2850 
] 0.489 0.4814 0.490 0.4524 0.4845 0.4874 
2 0.597 0.5865 0.594 0.5700 0.5903 0.5929 
3 0.684 0.6634 0.671 0.6525 0.6680 0.6702 
4 0.744 0.7258 0.734 0.7182 0.7310 0.7329 
5 0.799 - 0.787 0.7737 0.7847 0.7864 
10 1.006 - 0.987 0.9747 0.9817 0.9829 


The Lighthill values have been reproduced from his paper [8] and are the same as those 
given by Eq. (6.5). The “Imai’”’ values are those calculated by use of the author's full 
expression (6.1), while the ‘“Imai*”’ values are those calculated by the simplified expres- 
sion (6.2). It will be seen that the values calculated by the author’s formulas, both full 
and simplified, are in agreement with those computed by Levy within 1 per cent, except 
for the case of n = 0, where the difference is less than 2 per cent and 4 per cent, respec- 
tively, for the full and the simplified formula. By performing detailed examination of 
the accuracy of his numerical computations, Levy has estimated the accuracy of his 
results to be within 2 per cent. Therefore, our formula (6.1) and (6.2) may be considered 
to be perfectly satisfactory. Moreover, because of the simplicity, the simplified formula 
may be especially recommended for practical purposes. 

Further it may be observed that Lighthill and Fettis’ formula is better than Schuh’s, 
which, however, is fairly good except for the case of n = 0. This rather high accuracy 
of the formulas of Schuh and of Lighthill and Fettis is naturally due to the fact that the 
simplifying assumption u « y is very accurate for the Blasius flow. 

7. Stagnation flow. Similar comparison is made for the case of the stagnation flow 

3 = 1), taking o = 1 and 10. 

Table 2 shows the case of o = 1. Here the “Levy” values have been reproduced from 
his paper, while the ‘““Levy*’’ values are those computed by the present author using 
the empirical formula (5.21). The “Fettis” and “Schuh” values were also computed by 
the author by use of Eqs. (5.10) and (5.15) respectively. The ‘Imai’ and “Imai*” 
values are those calculated by the author’s formulas (4.19) and (4.23) respectively. 

It will be seen that both the Fettis and Schuh values are not in very good agreement 
with the Levy values. It is interesting to note that the Schuh values are better than the 
Pettis values except for small values of n, in spite of the fact that the former is an approxi- 
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TABLE 2 


Local heat-transfer coefficients for stagnation flow (8 = 1). Prandtl number = 1. 

n Levy Levy* Fettis Schuh Imai Imai* 
0 0.5603 0.581 0.661 0 0.584 0.539 
0.25 0.6345 0.640 0.735 0.492 0.645 0.629 
0.5 0.6979 0.692 0.799 0.620 0.705 0.701 
1 0.8116 0.783 0.905 0.782 0.809 0.814 
2 0.9647 0.932 1.065 0.985 0.971 0.980 
3 1.089 1.055 1.189 1.127 1.096 1.106 
4 1.192 1.161 1.292 1.24] 1.200 1.209 
5 ~ 1.256 1.382 1.337 1.289 1.299 
10 - 1.630 1.693 1.684 1.622 1.630 


mation to the latter. As has been remarked by Schuh, this is probably due to the fact 
that the following two influences cancel each other: (a) approximating the velocity 
profile by its tangent at the wall increases h for 8 > 0; (b) neglecting the middle term 
in Eq. (5.7) reduces h. 

On the other hand, both the simplified and full formulas of the author give results 
in satisfactory agreement with the Levy values. It will be observed that the approximate 
“‘Tmai*’’ values are even better than the ‘‘Levy*” values calculated by the rather com- 
plicated empirical formula. 

Table 3 shows a similar comparison for the case of rather high Prandtl number of 
o = 10. The situation is very similar to the above case of o = 1. Here the superiority 
of the Schuh formula to the Fettis for n > 1 is rather more conspicuous. However, our 
simplified formula is seen to be in closest agreement with Levy’s computation. 

8. Conclusion. Summarizing the result of the above comparison it may be stated 
that our approximate formula 


ha/k l (ty * 1(2/8) , ; oe B 
5 : ——_— 35 + n(2 — B)} (ona) ° — 
t.2/p - ) B) é | 3 “(4 ‘ ) 2 + 4 10a 


TABLE 3 
Local heat-transfer coefficients for stagnation flow (8 = 1). Prandtl number = 10 
n Lev) Levy* Fettis Schuh Imai* 
0 1.286 1.316 1.424 0 | .255 
0.25 1.451 1.449 1.584 1.061 1.449 
0.5 1.590 1.567 1.721 Lay 1.603 
] 1.818 1.773 1.949 1.684 1. 847 
2 2.159 2.111 2.295 2.122 2.204 
3 2.418 2.389 2.562 2.429 2.476 
j 2.630 2.629 2.784 2.673 2.699 
5 2.844 2.976 2.880 2.891 
10 - 3.691 3.649 3.628 3.606 


is the best of the formulas hitherto proposed in view of the fact that it shows the effect 
of the parameters @ (or m), n and o in a simple functional form. Especially it will be 
interesting to note that the exponent of o is invariably 1/3 and that the effect of the 
curvature of the velocity profile at the wall appears as an additive constant —8/10a 


in the formula. 
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Elements of partial differential equations. By Ian N. Sneddon. McGraw-Hill Book Co., 
Inc., New York, Toronto, London, 1957. ix + 327 pp. $7.50. 


In the author’s own words: “the aim of this book is to present the elements of the theory of partial 
form suitable for the use of students and research workers whose main interest 


differential equations in 
solutions of particular equations rather than in the general theory.’”’ An 


in the subject lies in finding 
attempt has been made to touch upon subjects from the entire field of partial differential equations. 
There are six chapters in all. Chapter 1 is introductory in character and deals with ordinary differential 
equations in more than two independent variables. Of special interest here is the treatment of Pfaffian 
differential forms, Carathéodory’s theorem on integrability, and its application to thermodynamics. 
Partial differential equations of the first order occupy Chapter2. Among the subjects covered are Cauchy’s 


problem for first order equations, linear and nonlinear equations, Cauchy’s method of characteristics, 


Lagrange’s and Charpit’s methods and applications, for example to stochastic processes and to birth 
and’ death processes in bacteria. Chapter 3 is concerned with a preliminary discussion of second (and 
higher) order equations, preparatory to the last three chapters, which consider in more detail the three 
main types of second order equations. One finds a treatment of linear partial differential equations 
anonical forms for second order equations, characteristics of equations in 
two and three independent variables, Riemann’s method of integration, and integral transforms. Chapter 
4 deals with Laplace’s equation, the prototype of the elliptic second order linear partial differential 
equations. The Dirichlet, Neumann and mixed boundary value problems are dealt with, Green’s func- 
tions,.Kelvin’s inversion theorem, and allied topics. The typical hyperbolic equation, the wave equation, 
t matter of Chapter 5. One, two, and three space dimensions are considered; the 


constitutes the subject 
Riemann-Volterra solution of the one dimensional equation, and Helmholtz’, Weber’s and Kirchhoff’s 
method 


general solutions in two and three dimensions. There is a brief introduction to Marcel Ries 


ficients, 


with constant coef 


of solution, based on Riesz’ generalization of the Riemann-Liouville integral. The concluding Chapter 


6 has the heat equation, the typical parabolic equation, as its topic. There is an appendix on systems of 
surfaces, meant to provide a brief outline of some of the properties of systems of surfaces used in Chapter 
2. The use of this book as a text is facilitated by the presence of a large number of problen s, together 
with the solutions of all odd numbered ones. 


The plastic methods of structural analysis. By B. G. Neal. John Wiley & Sons, Inc., New 





York, 1956. xi + 353 pp. $7.50. 

“The fundament defect in the orthodox elastic de sign t chnique can now be exposed. It has been 
seen that : simply supporte 1 beam carrying a uniformly distributed load .. . has a load f tor against 
collapse of 1.75, whereas 1 similarly loaded put fixed-ended beam were designed elasti illy to the 
same working stress its load factor against coll ipse would be 2.34.”’ With these words midw iv thro igh 
the introduction, Professor Neal embarks on a careful and detailed study of the status of beams and 
frames. His own work and the work of former colleagues at Cambridge and at Brown pl i central 
role in the text along with reference to the many other workers in this country and abroad. The book 
is intended Aas § step by tep zg tide through most ol the complexiti s of the actual calculation of ollapse 
loads, deflection estimates, effects of shear and axial force on limit moment, minimum weight design, 


and problems of shakedowr 
Mathematics is ke pt to the barest minimum even in the proofs of the theorems which are releg 
of buckling, 


ited 


to the Appendices. Except for providing the necessary warnings about the possibility 
fiects, nor homogeneity of material and similar difficulties, the author is n 


3auschinger effects, leflected 


in his purpose. He demonstrates the power ind the relative simplicity of plastic methods of analysis 


red the prac ticing enginee! who should read all of the text and check his solution 


for frames encounte 
to the probl ms with the nswers in the back of the book. 


D. C. DruckER 


(Continued on D- 72) 





ASYMPTOTIC SOLUTIONS OF TOROIDAL SHELL PROBLEMS* 


BY 
R. A. CLARK 


Case Institute of Technology and Massachusetts Institute of Technology 


1. Introduction. A few years ago a method of asymptotic integration was developed 
by I. Reissner and the author for a non-homogeneous differential equation involved in 
certain problems of the theory of thin elastic toroidal shells. Using this method solutions 
were found for the problem of bending a curved tube (considered as a sector of a shell 
[4] and for the mathematically similar problem of a toroidal expansion 
joint subject to an axial force [2]. The main purpose of this paper is to indicate how the 
method used before can be refined and how the approximate solutions previously obtained 


f a mere. 
or revolution 


for these problems can be generalized and extended. 

It is shown, in particular, that the previous solutions are the leading terms of asymp- 
totic expansions in inverse powers of a large parameter. Higher-order terms or corrections 
of two different types are obtained for toroidal shells of circular cross section and uniform 
thickness. Some of the results are found by specializing solutions derived first for toroidal 
shells of arbitrary cross section and possibly varying thickness, assuming only that the 
shell has 

2. Basic equations for shells of revolution. The problem of determining rota- 
tionally-symmetric stress distributions in thin elastic shells of revolution may be reduced 
te that of solving two coupled second-order differential equations. Using the formulation 

tially the notation of E. Reissner [9] the differential equations for small deforma- 


has a plane of symmetry perpendicular to the axis of revolution. 


anda @ssel 


tions may be written in the form 


X" — 0X + woY 
Y’’ — 6.Y — nweX 


F, 
(1) 
G. 


Here primes denote differentiation with respect to an independent variable ~ which is 
the parameter in a representation of a meridian of the middle shell surface by equations 


of the form 
r = r(§), z = 2(g), (2) 


where r, 6, z are the cylindrical coordinates of a point on the middle surface. Representa- 
tion (2) is assumed to be such that r and z have continuous derivatives up to at least 
third order and also such that the quantity 


a = [(r’)’ + (@)*]'” (3) 


*Received February 15, 1957. Sponsored in part by the Office of Ordnance Research, U. S. Army, 
under contract with Case Institute of Technology and in part by the Office of Naval Research, U. 8. 
Navy, under contract with Massachusetts Institute of Technology. 
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be 
i?) 


never vanishes. The coefficients and loading terms in Eqs. (1) are given by 


»!\2 (r’D/a)’ 1] (rD/a)’ |° 1} (rD/a)’ |’ 
) = — 7S as a 7 c u 7 a" ’ 
7 (" ” “(rD/a) r {2 a) | ’ AS a) | 
) ik 
Li ( 


Pas (“") : (r’/aC)’ r/aC)’ : , 1} (r/aC)’ |’ 
= i ik ae 1\ oy , 2627, 64s 
. I (r/aCl) *~ r/aCl) 21 (r/al 
) » az’ i a \'”? . }) 
u® = [12(1 — v*)]'? —, f= r(- - (rV), 
rh wr D 


/ ¥\1 l 2 roe al Y\/ Pm 
- Lat | r Ji re (2° /aCl) ‘ Z _— 
G = r'\—) k + ( .) 9] oe te IV) += I 
Es al oe yg (r/al) r 


where C = Eh, D = Eh*/12(1 — v°), E is Young’s modulus, »v is Poisson’s ratio and / 
is the wall thickness of the shell. Note that the quantity u is inversely proportional to 
the thickness A and so is large compared to unity if h is sufficiently small compared to 


other dimensions. 

Quantities pz, and py denote radial and axial components of surface load intensity, 
while V is an axial stress resultant determined by overall static equilibrium. Normal 
stress resultants NV; and N, , transverse shear resultant Q, and stress couples M; and 
M, are all defined in the usual way. The components of displacement of a point of the 
middle shell surface in the three coordinate directions r, 0, z are denoted by w, v, w, 
respectively, while 6 denotes the angle of rotation, due to deformation, of a tangent to 
a meridian. Although all other quantities are assumed to be independent of @, the cir- 
cumferential component of displacement v is allowed to have the form 

v = kré (k = constant) 5) 
which occurs when a sector of a shell is subjected to pure bending in a place perpendicular 
to the axis of revolution (see [8]). Once differential Eqs. (1) have been solved subject 
to appropriate boundary conditions, the various quantities may be determined from the 
following relations, 


ry | rapy dé 
rN; r(s } 4 +26V, rQ = 2’(- ) } re 4 rV) 
a a ra ax 
aN, = (%)"] vr — Lol" | rope 
ais 0)" a [2-4/2 ' 
ems = (2) "ux + |F - Sie bp 
B { >) As u - Ve — yN 
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3. Asymptotic integration when @ has a first-order zero. Approximate solutions 
of differential equations (1) can often be obtained by methods of asymptotic integration 
provided the parameter u» in these equations is sufficiently large compared to unity. 
The form and completeness of the results depends on the nature of the coefficient func- 
tions 0, , 9, , and ®. For instance, if these functions are all bounded and if # is also 
bounded from zero, asymptotic expansions in » can be found which represent solutions 
of either the non-homogeneous differential equations or the corresponding homogeneous 
equations and which involve only elementary functions. Such expansions have been 
given for general shells of revolution by F. B. Hildebrand [5]. If ® vanishes at some 
transition point or if 6, and ©, have at most second-order poles, as they usually do 
on the axis of revolution, asymptotic solutions can be obtained for the homogeneous 
differential equations using methods developed by R. E. Langer [6a, b] which involve 
Bessel functions. In general, Langer’s methods yield only the leading terms of asymptotic 
expansions of solutions but may be used to find higher-order terms if the shell is such 
that Eqs. (1) can be reduced to a single second-order differential equation. Reduction 
of Eqs. (1) to a single second-order equation is possible whenever the shell thickness 
and the meridianal curvature are constant or, in general, whenever the ratio 
(0, — 0,)/® is constant, as has been shown recently by P. M. Naghdi [7] who also 
lists the form of Langer’s asymptotic solutions. 

The present paper is concerned with asymptotic solutions of non-homogeneous 
equations (1) assuming that all coefficient functions and loading terms are bounded 
functions of £ for the shell or portion of shell being considered and that ® vanishes at 
some isolated point. According to Eqs. (4), ® vanishes whenever z’() vanishes or when- 
ever the tangent to a meridian is perpendicular to the axis of revolution. Define £ so 
that £ = 0 corresponds to the point where ® is zero and let the range of £ be restricted, 
if necessary, so that ® is bounded from zero away from & = 0. Expanding about = 0 


P(e (1211 — p?)]'/%Xe az" yo] az 2(2 ) — # ” 
ee hs ”)] (| i. + E = rh r<0 23 M4 


It follows that the coefficient of the first power of in expansion (7) is different from 
zero only if z’’(0) + 0 or only if the radius of curvature of a meridian is finite at § = 0. 
The function u® then is said to have a first-order zero and this is the only case we will 
consider. The parameter u may be conveniently defined now as the value of the coefficient 
of € in (7). Using a zero subscript to denote the value of a quantity at £ = 0 we may 


we have 


write 


, 3 
11 — 1/2 AkRo a 9 __ 42) 71/2 _ & . ” 
mn (12(1 v)] {12(1 v)] Thole,’ (8) 


Toho 


where FR, is the radius of curvature of a meridian at § = 0. Dividing expansion (7) by 
formula (8) we now have 


(é) t+ At? + --- (9) 
where 
ia 0 o ho 
oe ee (10) 
2 “0 ay To om 
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If there are no “rapid changes” in the shell surface or in its thickness, functions 
6, and 0, are of order unity outside some neighborhood of the axis of revolution and 
so are small, generally, compared to u® provided ux is sufficiently large. This suggests 
that approximate solutions of differential equations (1) can be found by neglecting 
terms in 9, and 0, . The resulting differential equations may be interpreted as the real 
and imaginary parts of a single complex differential equation of the form 
Z" — ipoZ = ff (11) 


where 


Z=X4+1Y, f=F+7. (12) 


Thus, if Eqs. (1) cannot be reduced to a single second-order differential equation in an 


exact manner, they may be approximated by such an equation. The nature of the error 
in the approximation is discussed later on. 

We consider now the problem of obtaining an approximate solution of differential 
equation (11) when uz is large. If we assume that a solution can be represented by an 
expansion in inverse powers of uw, we obtain the following formal solution 


, if, 1 (f\" < 

Z=- - + in +: eee (13) 
uP m Pm \® 

viously is not valid in some neighborhood of £ = 0 where ® vanishes. However, 

terval where © is bounded from zero a finite number of*terms of (13) does 

solution of (11). The first term of ( 13) will be called the membrane 


which ol 
over al 
furnish an approximate 
solution of (11), or of differential equations (1), since retaining only this term leads to 
the same results as the so-called membrane theory of shells. An approximate solution 
of Eq. (11) which is valid at € = 0 and which reduces to the membrane solution away 
from & 0 may be constructed as follows. First, in the immediate neighborhood of 


£ = 0 differential equation (11) may be approximated by 


n 


Z”’ — iwtZ = f(0), (14) 
where we assume that f(0) ¥ 0. A solution of (14) is given by Z = yp *’*f(0)T(u'”*E) 
provided the function T(x) is a solution of 

T’'"(x) — 1xT (x) = |, (15) 


(Henceforth, primes are used to denote differentiation with respect to the written or 
rgument of a function.) Let 7(x) stand for that (unique) solution of (15) 


implied 
» o and consider the expression 


which has the behavior T(x) ~ 7/z as | a 
Z(t) = 2/3 £(E)T (7). (16) 
S/ K J\S sd 


Near ¢ = 0,  ~ é£ and f() ~ f(0) so that Z, approximates a solution of (14) or (11). 
< YU) 
Away from é = 0, ’°®@is large for large u so that Z, ~ zf/u&. Consequently, Z, furnishes 
d § L . 1 ? 
an “interpolation” between a solution of Eq. (14) and the membrane solution of (11). 
The real and imaginary parts of the required function 7'(x) and its derivative are tabu- 
lated briefly in [2]. Also the following representations hold, 
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T(r) = | exp (— iri — 14°) dt 


a ( 


mitoses 1\ 2’ a 
=o = > (~ 7g" (2+) 2 (17) 


r)! 


"(By 
~ » 3'n 


1 
Ix 


It can be shown that the real and imaginary parts of Z, are truly asymptotic solutions 
of differential equations (1) in the sense that the relative error approaches zero uniformly 
as u — ©. Even so, Z, may not yield sufficiently accurate results when yp is only mod- 
erately large making it desirable to have higher-order approximations to solutions of 
Eqs. (1) or (11). Recently [8] complete asymptotic expansions were found which repre- 
sent a solution of a non-homogeneous second-order differential equation such as (11). 


They are of the form 


[2n/3] 
Z. t) = a? 3 Zz. bu n/3 : Gon on—ae($)T" 2n 3h) (4) (18) 
n=0 k=0 


ZA) = w??T(x*) YS wa, (&) +? Dd wd, ) 
n=0 n=0 
(19) 
+p *?T"(2) > ue, (€), 


n=0 
vhere T'*’(2) denotes the kth derivative of 7 with respect to x and 
z= p'*(8w/2)””, o = [ %'’*(&) dé. (20) 
70 


The first term of expansion (18) is similar to the expression for Z, , while expansion (19) 
analogous to expansions obtained by Langer [6c] for solutions of a homogeneous 
second-order differential equation. The various coefficients in either expansion may be 
successively determined so that the expansion formally satisfies differential equation (11). 
The expansions do not converge, in general, but a finite number of terms of either fur- 
nishes an approximate solution with an error of the order in yu of the first term neglected. 
Depending on the circumstances one expansion may be easier to use than the other but 
one may also interpret expansion (19) as simply a “rearrangement” of the terms in (18) 
accomplished by using differential equation (15) to express higher derivatives of T in 
terms of 7 and 7”. If Eq. (11) is only an approximation of Eq. (1), it would not be con- 
sistent to retain arbitrarily many terms in Z, or Z; , considered now as approximate 
solutions of differential equation (11). In [3] it is shown that terms of relative order 
u ’* in either Z, or Z; depend on terms such as those in 0, and 0, in differential equa- 
tions (1). This fact leads to the result that when uw is large the relative error in approxi- 
mating a solution of Eqs. (1) by a solution of Eq. (11) is of the order u~*”*, making it 
eenerally inconsistent to retain more than the first two terms in the expansions for Z, 
ind Z, . But if Eqs. (1) can be reduced to a second-order equation in an exact manner, 
then theoretically asymptotic solutions with arbitrarily many terms can be found. 
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\. 


For later reference we now list the explicit form of the first three terms in the two 
expansions for a solution of differential equation (11), namely 


i Ds reaaiseeseis t As niemiees 
ZAE) = a SO T(x) + -OV(OS)'T’'(x) 
M rm 9 
(21 
+ tq (Q'Q"FT"(x) — 4Q°(Q7(Q"f)''xT’"(2)}, 
bu 
agit -_ dine gis f( 0) t f(0)T’(x) [ a” 
iit) = - fOQOT (x) Geet = ey ee: Beek, ——sE Cir, 22) 
Z3(E a O)QOT (x) + : E ay 2 i 200 J, Qe” dé ( 
where 
QE) = (8w/2)""/8”, Q) = 1. (23) 


Note that definitions (23) and (20) imply x = u'*@Q* and also that dx/dé = y'*/Q’. 
If 7'(x) is replaced by its asymptotic form 7/xz when z is large both (21) and (22) reduce 
to the membrane solutipn 7f/u@ plus terms of the same order as the second term in 
expansi 3). When zx = O(1) it follows that £ = O(u~'”*). Then if f, &, and Q are 
expanded in powers of & and terms of relative order » ** are neglected, expressions 
(21) and (22) both reduce to 


; 7 oe , ee. . ee cai 
Z T (us) + = $= sO) [WET (ut) — wT (u'*8)] 
i ei \ Vv D4 
\ (Z 


| 
+ 2 f’(0) — 2A,f(0)]?, 


where A, is defined by Eq. (10). Although expression (16) for Z, also reduces to the 

membrane solution when z is large, for x = O(1) it reduces to an expression that agrees 

with only the first term of (24), showing that the relative error in Z, is of the order 
73 in general. 

The asymptotic solutions described above are applied to shell problems in the following 
sections. 

4. Toroidal shells of arbitrary cross section. We shall consider two particular 
problems involving toroidal shells. One problem is the bending of a curved tube assuming 
that the tube forms a sector of a shell of revolution. The other is the deflection of a 
toroidal expansion joint due to an axial force, where the expansion joint consists of a 
toroidal shell slit around its inner circumference and joined at the edges to sections of 
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Fic. 1. Cross section of toroidal-shell expansion joint. 
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cylindrical pipe (see Fig. 1). As will be seen, the two problems are mathematically quite 
similar. Asymptotic solutions of both have been given previously for toroidal shells of 
circular cross section and uniform thickness [2, 4]. For the present we assume only that 
the shell has a plane of symmetry perpendicular to the axis of revolution. 

If P denotes the axial force applied to the expansion joint the loading terms in differ- 
ential equations (1) become 


ee -~ i r’e! (2’ /aC)’ Lay" se 7 
F=1 (3) 2r’ G= E +» (r/aC) | al on (25) 


Generally G is small compared to F and if coefficients 6, and ©, are neglected in Eqs. 
(1) it is consistent to neglect G also. The right side of differential equation (11) is then 
given by the first of Eqs. (25). 

For the problem of a tube or sector of a shell of revolution bent by end moments 


applied in a plane perpendicular to the axis of revolution we have 
C 1/2 
F=0, G= 7(2¢) k. (26) 


The constant k is the same as that in Eq. (5) and is related to the change in curvature 
AK of the neutral surface of a tube by k = r, AK, where r, is the radius of the cylindrical 
neutral surface. For tubes with closed cross section the applied moment m is given by 


m= [ — 7,)Noea dé. (27) 


After differential equations (1) have been solved and stress resultants obtained in terms 
of k, the relationship between k or AK and the applied moment m may be established 
using (27). 

Let the plane z = 0 be the plane of symmetry for either problem and consider only 
the lower half of the shell, as shown in Fig. 1. Let £, and £ correspond to the inner and 
outer points, respectively, where a meridian meets the plane of symmetry. For the tube 
bending problem it follows from symmetry that 8 = 0 and Q = 0 at both é, and £, . 
These conditions also hold at £ for the expansion joint but may not be the actual edge 
conditions at £, . However, we shall assume that 8 = 0 and Q = Oat &, for the expansion 
joint also, since it can be shown that axial deflection and maximum stresses are not 
significantly affected by a change in the conditions at £, . Then, according to Eqs. (6) 
and (12), the boundary conditions for either problem are equivalent to 


Zi) = 0, Zé) = 0. (28) 


The closed meridians of the toroidal shell are assumed to have continuously turning 
tangents which implies that r’ vanishes at é, and & . Since the assumed loading term 
f(£) for either problem involves r’ as a factor, it follows that the membrane solution 
if/u® vanishes at £, and £, so that asymptotic solutions Z, , Z, , and Z; all satisfy bound- 
ary conditions (28), at least to within terms of the order of the error involved. This 
means that no solution of the corresponding homogeneous differential equation is required 
and that approximate expressions for the various stress resultants and couples and 
displacements can be obtained from the real and imaginary parts of Z, , Z, , or Z; together 
with Eqs. (6). To evaluate the stresses and displacements as functions of one must 
first tabulate the functions ®, Q, and x given by Eqs. (4), (20), and (23). These functions 
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in turn depend on the exact shape of a meridian of the middle surface and on any varia- 
tion in the shell thickness. Nevertheless, maximum stresses and other quantities of 
interest can be obtained explicitly without assuming any specific shape for the shell 


cross section. 

Since T(x) ~ i/zx for large x the most significant values of Z as given by (16), (21 
or (22) occur where x O(1) or £ = O(u"*”*). For this range of € expression (24) may 
be used as an approximate solution of differential equation (11) with a relative error 
of order »~*’*. When the real and imaginary parts of (24) are substituted in [nqs. (6) 
it becomes apparent that for large values of u the most significant stresses are the cir- 
cumferential direct stress op and the meridianal bending stress o;, , defined 


No 6M, 


(99) 


The maximum numerical values of these stresses and their location ean be determined 
in a straightforward way using the following information about the real and im iginary 
parts of T(x), denoted by 7,(x) and T,(x). The maximum value of 77(x) occurs at x 0 


where one has 


7 (O) = 1.288, TiO) = 6, T’(0) 0, THO) = .939. (30) 
The maximum value of 7'’(x) occurs at x 1.225 where 
T (1.225) = —.692, 7,(1.225) = 817, 77(1.225) = .753, 7T%(1.225) = .214 (31 
Also, 7.(x) is an even function of zx, 7;(x) an odd function of x. 


For the expansion joint problem we find after various calculations that 


way 2 
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As before a zero subscript denotes the value of a quantity at = 0. The plus and minus 
signs in the formula for o;, correspond, respectively, to the signs in €, and the constant 


A, in &, is defined by Eq. (10). 


1958] TOROIDAL SHELL PROBLEMS 55 
The total axial deflection or spread of the expansion joint due to bending only is 
given by 
fs 
6 = Qlw(E.) — w(t] = —2 [ r'B dé. (35) 
Jt, 
Since the integral extends over the full range of £, expression (24) is no longer appropriate. 


If the first two terms of Z, are used instead and 8 is determined according to Eqs. (6), 
the expression for 6 takes the form 


i ( Q ) ™ [" ( a ) iis ; 
é —~— |(— r\—) Q°T,(2x) dz, 36) 
mu \t%oDo, Ta rD é : ( 
where x, and x, are values of x corresponding to é, and &j . Since T(z) ~ 2/z* for large z, 


, ’ . - . y+ -1/3 
the major contribution of the integrand in (36) occurs for x = O(1) or & = O(u"'”*). 
Also, interpreting 7'(x) as a Fourier transform according to the integral representation 


given in (17) it follows that {% 7,(~7) = —2/2. Using these facts and introducing the 
explicit form of u as defined in (8) it can be shown that 
‘ 2 nv2 oF ] a 
6 = [12(1 — v*)]'* =a E + o( :)|. (37) 
Ehs mM 


This expression and the leading terms in expressions (32) are exactly the same as the 
corresponding results for a toroidal shell with a circular cross section of radius R, and 
« mean shell radius ry , showing that when uz is sufficiently large the axial deflection and 
maximum stresses depend on the exact shape of the cross section only near the point 
where the slope dz/dr of a meridian is zero. 

For the tube bending problem results similar to (32) can be obtained with Fk in 
place of P/h5 . However, it is more convenient to obtain first the relationship between 
the applied moment m and the constant k. Substituting for aN, from Eqs. (6), with 
Y taken as the imaginary part of Z,; when f = ikr’(aC/r)'”’, and integrating by parts, 
finally reduces formula (27) to 


Ne - m OM 1/2 ate C 1/2 a , 
m = — Hrs (a ) | (2) Q°T,(x) dz, (38) 
M To da ¥ 


a 


which is analogous to expression (36). Evaluating the integral in (38) in the same way 
as in (36) and introducing u from Eq. (8), yields the result, 





2r 1 
= To a RohSER 1 + (2) |: . 
m 2d — >) 3 Roh x O Pa (39) 
We recall that k = r, AK, where AK is the change in curvature of the neutral surface 


and r, is the radius of the neutral surface. Writing r, = r(&,) the value of &, is by definition 
i solution of the equation N, — vN; = 0. From expression (24) and Eqs. (6) one obtains 


p= 1388] 4, 4a 498] 4 of), (40) 


Le To 


where A, is given by Eqs. (34). It follows from (40) that r, = rofl + O(u7?”)], so that 
expression (39) can also be written in the form 


™m 
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This result may be compared with the equation m = EI AK which holds for a straight 
tube according to elementary beam theory, where J is the moment of inertia of the cross 
section about an axis through the centroid. Writing 

m = pEIAK (42) 


for a curved tube serves to define the quantity p as a rigidity factor which has the value 
unity for a straight tube. Comparing (41) and (42) it follows that 


or roRoh? ( 1 )] 2 rate | ( 1 )] 
p [--- =o i ———_ - ) 3 — — “SO: (43) 
HFT |1 + o( 4s) a P| + Olan 13 


12(1 — » 


Formula (43) is remarkably simple and shows that when uz is sufficiently large the effec- 


tive stiffness pEI of a curved tube depends on only the dimensions rp , Ry and hy associated 
with the point of the cross section where £ = 0. For a tube with circular cross section 
of radius R, and uniform thickness hy , J = tRoho . Also the quantity a may be identified 
with R, so that the rigidity factor then becomes simply p = 2/u when z is large. 

The maximum stresses which occur in the bending of a curved tube may be compared 
with the maximum stress in a straight tube of circular cross section with radius Ry and 
uniform thickness hy) given by 

- mR _ m 14 
eT te. si 


Introducing the expression for m given in Eq. (39) and determining the maximum 
- 1. \ : 


stresses in the same way as in the expansion joint problem we find 


_ — + §62(1 —yv ( 2 |) E 1. 2.144, -; sae 


o f TiNel 
1\] 
(ts (45) 
i i ye ) ) 


where now 


a — % . (46) 
429 (J 

—“33 + O|-) 
m u/. 


The constant Ag is defined by Eq. (10), A, and A, by Eqs. (34). The leading terms in 
(45) are the same as those for a tube with circular cross section of radius R, . Thus, 
for sufficiently large 1 the maximum stresses due to bending of a curved tube of arbitrary 
cross section can be obtained by considering an ‘“‘equivalent”’ tube of circular cross section. 

5. Toroidal shells of circular cross section. We now consider the form which the 
results in the previous section take for toroidal shells or curved tubes of circular cross 
section and uniform thickness h. We shall also list additional terms of relative order up ~’° 


For shells of the special type now being considered parametric equations (2) may be 


taken as 
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r=act bsiné, z= —bcosé, (47) 

so that rp = a, Ry = b, a = b. The parameter uw and a new parameter \ are given by 
b re 

A= - = [1211 — »’)]'? —- 48 

> #*OR eT 2 (48) 


Previous asymptotic solutions of the special problems being discussed did not take 
into account the effect of terms in \ which was assumed to be negligibly small. 
For the expansion-joint problem the largest values of the meridianal bending stress 


Trp are now given by 


bh 1/3 > 
o:n(fs) = +.482 ( *) P [ += 86 a Ae o(4 :)], (49) 
a h mn u 
where we have taken vy = .3 and where 
225 ” 1 ; 
te = tte — 79-5 + o( +). (50) 
BM rm “ 


It follows that the largest of the two values listed occurs at the “inner’’ point where £, 
is negative. For the tube-bending problem the only stress formula which involves X is 
the following for the maximum circumferential direct stress, 


: er Pee 1 . 
dep\tp) = +.377p ‘sou I + .58 38 + o( 1). (51) 
a M 
where again we have taken vy = .3 and where 

1.225 ~ A ] _ 
Ep = +5 — O97 -m + o(*). (52) 

ra rn M 
From Eqs. (49) and (52) it is easily seen that terms in \ play a minor role when d is 


small compared to u’’*, which does not necessarily mean that \ is small compared to 
unity as had been assumed in earlier work. However, the correction due to retaining 
terms in \ may certainly be significant. For instance, when \ = 1/4 and uw = 8 the 
correction to the maximum bending stress in the expansion joint, as given by Eq. (49), 
amounts to 10 per cent. 

The relative size of the correction terms in \ given above does not furnish any estimate 
of the accuracy of asymptotic formulas since at least part of the error term is independent 
of \. Theoretically, exact expressions can be obtained for the terms of relative order 
uw (or higher) but the necessary calculations become quite involved. However, if all 
terms in \ are neglected higher-order terms in » can be obtained without great effort. 
Neglecting all terms in X is equivalent to approximating differential equations (1) by 
the single complex equation 


Z’’ — insin §£Z = fy, cos, (53) 
where {, = (b°/aD)'” P/2m for the expansion-joint problem and f, = i(b*hE/a)'” k 


for the tube-bending problem. Neglecting terms in A is also equivalent to assuming that 
the various stresses and displacements are symmetrical with respect to = 0. An approxi- 
mate asymptotic solution of (53) which includes terms of relative order u~*” is given 
by either expression (21) or (22). To determine the most significant stresses one again 
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58 . 
f ‘ E es : ; ; 
restricts attention to the range where & = O(u ) and expands all functions in powers 


of £. In this way the following results are finally obtained. 
For the expansion-joint problem 


yal bh\' , 42 yf 
Fiien ~*~ 60 = 3401 — 7 (%) - [ _ 2 4. of )| 
\a h be \ 
_ UT ie ia (2h) ) Ph _ 054, o(*) | 
OrtB.ma ( om ~ =; a y?)'46 a’) h? ue Ma u ’ 
995 7 
‘i | 1335 _ O17 + of |. (55) 


(54) 


rn KM Me 
) bP 4 1 
§ = [121 — )P* 5114+ Of )|. (56) 
Eh \M/ 
Analogously, for the tube-bending problem 
- 054 4 
06D.ma Seaplfm) = + 4 (kh OsMm l — ” “t- o( - 
bh Me (R-7 
o6) 
814 - 342 mi 
OrB.ma o:3(0) = - mse osm) | — —z + o(+) 
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where £,, is again given by Eq. (55). Note that there is no “correction” term of order 
uw” in formula (56) for the axial deflection 6 indicating that in an exact solution this 
curved tube the analogous 


term is probably proportional to \. For the bending of 
result for the rigidity factor is p (2/u) [1 + O(1/y)]. 
In [4] where the leading terms in the above results for the tube-bending problem 
were first obtained, the asymptotic formulas are compared with numerical results 
calculated by L. Beskin [1] using a strain-energy principle in a manner equivalent to 
finding an exact trigonometric series solution of differential equation (53). In Tables 1 
and 2 we have listed the numerical values obtained by retaining either one or two terms 


in formulas (57) and also the values obtained by Beskin. 
TABLE 1. 
Va »f Cop max/Osy 
377? 377 p?/3-.020 Be 
1 OO 31/4 357 949 
2.50 694 674 3350 
3.54 875 R55 RS 
5.00 1.10 1.08 1.010 
7.07 1.39 1 $7 1.277 
10.00 1.75 1.73 1.680 
12.50 2.03 2.01 1.99 
17.68 2.56 2.54 2.54 
25.00 3.22 3.20 3.26 
35. 36 + 06 4.04 4.08 
50.00 5.11 5.09 5.09 
70.71 6.45 6.43 6.43 
100.00 8 13 8 1] Ss. ie 
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TABLE 2. 


Values of —~ ‘1 —- py? OrB ,max/TSM- 
8142/3 814y?/3-.278 Beskin [1] 
1.00 814 536 .432 
2.50 1.50 1.22 1.06 
3.54 1.89 1.61 1.47 
5.00 2.38 2.10 2.01 
7.07 3.00 2.72 .2.69 
10.00 3.78 3.50 3.51 
2.50 4.38 4.10 4.11 
17.68 5.52 5.24 5.25 
25.00 6.96 6.68 6.67 
55.36 8.77 8.49 .42 
50.00 11.0 10.7 10.7 
70.71 13.9 13.6 13.5 
100.00 17.5 17.2 17.2 


We observe that for o¢p max the values found using the two-term asymptotic formula 
coincide with the values obtained by Beskin for » > 50 and differ from Beskin’s results 
by less than 3 per cent for » > 10. Formula (55) for the location of the maximum yields 
values which differ from those given by Beskin by less than 2 per cent for u > 5, and 
by less than 1 per cent for u > 10, except for » = 50 when there is a difference of more 
than 3 per cent. 

From Table 2 for cez.max We see that although the one-term asymptotic formula 
does not agree as closely with Beskin’s results as does the one-term formula for o¢p, max ; 
the two-term formula yields values which differ from those found by Beskin by 1 per 
cent or less for u > 7, which is an even closer agreement than was found for o¢p, max - 
In fact, the two-term asymptotic formula for o¢,m.x differs from Beskin’s results by 
less than 5 per cent for uw as small as 5. 

Finally, it is of some interest to compare expressions (49) and (54). Unless \ < .1/y'” 
the correction to the one-term asymptotic solution due to retaining terms in \ may be 
much more significant than that due to terms of relative order u~*”’, at least as far as 
the dominant stress o;, for the expansion-joint problem is concerned. A similar remark 
can be made about the circumferential direct stress o»p for the tube-bending problem, 
although this stress is not the dominant one for this problem when uz is large. Also, the 
“-correction” to £, or €p in formulas (50) and (52) may be considerably more important 
than the second term in formula (55). 
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ON LIMIT ANALYSIS OF PLATES* 
BY 
WALTER SCHUMANN 
Brown University 


Summary. The bending of thin perfectly plastic plates of arbitrary shape under 
transverse load is studied. The collapse load for a concentrated force on such a plate is 
found to be 27 times the yield moment. 

1. Introduction. In the following an attempt is made to generalize the work of 
Hopkins and Prager [1]} on the load-carrying capacities of circular plates to non-sym- 
metrical cases. The basic equations for this purpose have been given previously by 
Hopkins [2] for a plate composed of a material obeying Tresca’s yield criterion. We 
now discuss in detail the different types of plastic regimes that can arise. After a brief 
consideration of the discontinuities in the moment field, the technique of limit analysis 
[3] is used to determine the limit load for a concentrated force. Upper bounds on the 
limit loads in non-symmetrical cases have also been obtained by Rzhanitsyn [4]. Only 
in very artificial examples has it been found possible to obtain the actual moment 
distribution and deformation mode during collapse. 

2. Equations of equilibrium in curvilinear coordinates. Let M, , M; be the principal 
moments; Q, , Q, the shearing forces; p the load per unit area; p; , p2 the radii of curvature 
and s, , 8, the arc-lengths of the principal stress trajectories (Fig. 1). Then the equations 
of equilibrium are 


aah! 1 eS 

om 4 aM _ 9, =0, 

Os, Pe 

OM, , Mi —- M _o uo, a 
O82 Pi 


dQ, IQ. Q 1 Q. 

: a a 

Os; O82 P2 Pi 
We now introduce orthogonal curvilinear coordinates u, v such that the principal stress 
trajectories are given by u = const. and v = const. Let 


ds = Edw’ + 2F dudv+ Gav (2) 
be the first fundamental form of the net of stress trajectories.** We have then 
P te) 1 @ es 4.¢ 
Os, - E' . Ou F OS> ai qG'” Ov R 
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Pi 1 I (3) 
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*Received March 28, 1957. The results presented in this paper were obtained in the course of re- 


search sponsored by the Ballistic Research Laboratories of Aberdeen Proving Ground under Contract 
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tNumbers in square brackets refer to the Bibliography at the end of the paper. 

**Because F = 0, it is not necessary to use here the tensor notation g;; . 
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Fria. 1 
where ¢ is the angle of inclination of the curves v = const. to the z-axis, and subscripts 


denote partial differentiation. Further, it is convenient to introduce the quantities 
» = (M, + M,)/2 and 6 = (M, — M,)/2. The first two of Eqs. (1) are then 


(w + 5), 25- 

<a 7 -—@Q Q, 

a Slag sl 

(+) 
©) 9s 

(w — 6 26°¢, 

a + ——- — g. = 0, 

qi’ Fi ais 


and the third equation of equilibrium becomes 


Q (Q.) G E 
< — t+ (), oz -t- (), — — 5 = —). (5 
= + — > 1 ar" “asc d ; 
To eliminate Q, , Q2 we introduce (4) in (5) and get 


Wen, Yr , __ 3 | (_ G ’ 7 EE} Buu Sy 

B+ G+ ca [aa + (aqe)o|+ F- ¢ 
7s 1 ( G ) i (. : E . 6 
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3. The yield condition. Tresca’s yield condition is shown in Fig. 2. For this con- 


dition only four essentially different types of regimes occur, namely two stress determined 
regimes, represented for example by the points A and B, and two kinematically de- 
termined regimes, represented for example by the sides AB and BC of the yield hexagon. 

regime A. As 6 = O, the angle ¢ is indetermirate. This regime can occur in a finite 


region only if the pressure p is zero (see also Prager [5)). 
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M x 
*4 ° Straight lines 
Hencky-Prandti net— ra 


for p=0 ‘ 
ee 





Fic. 2. 
Regime B. In this regime w = M,/2 and 6 = —M,/2, where J/, is the yield moment. 
Equation (6) becomes 
2) : myy\ 1/2 (7 
2M eu. = p (EG) ° (7) 
For the special case p = 0, we have 
us = 0 (8) 


This is the condition that the stress trajectories form a Hencky-Prandtl-net [6]. 
Regime AB. Yor this regime it is necessary to consider the velocity field. The principal 
curvature-rates are denoted by x, and x, . For moment states represented by points on 
0. Thus the surface into which the plate is 
deforming has parabolic curvature, and therefore one set of the principal lines of curvature 
is straight (see for example, Blaschke [7], p. 112). Because of the isotropy condition 
follows that one family of the principal stress trajectories consists of straight lines. 


the side AB, the flow rule requires x, = 






Asw— 6 M. M, the second of Eqs. (1) becomes 
Q, = 0. (9) 
V 
y 
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| ae 
Filx 
| = 
| 
| 
oe 


Fic. 3. 








64 WALTER SCHUMANN . [Vol. XVI, No. 1 


If we introduce p = p, and g as new coordinates (Fig. 3) the first of Eqs. (1) gives 


0(pM,) 
ee — My = 0Q, , (10) 
Op 
and the third of Eqs. (1) becomes 
O(pQ,) , 
SPE = 0. (11) 
ap - pp 1] 


This is taken from a manuscript by R. T. Shield. The integration of (10) and (11) leads 
to the expressions 


M, M, + Ay) + Bi) 4 [ ppdp — [ vp dp, (12) 
p p J d 
. A(g) 
Q, = Ay) _! | ve dp, (13) 
p p 


where A and B are functions of g only. The evolute e and the functions A and B are 
to be determined from the conditions of the problem, such as joining a region in regime 
AB to regions in other plastic regimes. 

Regime BC. In this case we have x, = — xz which is the definition of a minimal 
surface. On such a surface the lines of principal curvature form an tsometric net (Weather- 
burn [8]) and isometric coordinates u, v may be introduced for which FE = G. If z, y are 
rectangular Cartesian coordinates, then z = x + zy is an analytic function of x = u + w 
(see, for example [7] p. 179). Also ¢ = arg (dz/dx) = arg z’ = Im (log z’) is a harmonic 
function (Fig. 4). 

The equilibrium equation (6) becomes with EF = G and 6 = — M,/2 (from the 


vield ondition 

Aw — 2My,, = —pE, (14) 
where A is the Laplace operator referred to u, v. As E = x2 + yZ = | 2’ |°, andas gy 
is the real part of the analytic function 


Eq. 14 may be written 
gly"! ae (2’’)? : of 
Aw — 2M Re | nz "| = —ple’|’. (15) 


We remark that instead of ¢ the conjugated function y may be used. Equation (14) 
may be easily transformed to the principal shear-lines, which also form an isometric 
net. 

For the special case p = 0 one can find immediately two particular solutions of (14), 


w, = up,M, = uM,-Re =| 
. (16) 


Ww, = vy,M, = vM,:Im Fak 
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The general solution is therefore 
®W=w, tw, OF w= We tw, (17) 


where w, is any harmonic function. 

4. Boundary conditions. We suppose that the plate is at yield in bending near the 
boundary and discuss the three following types of boundary conditions: a) simply 
supported edge, b) clamped edge, c) free edge. 

a) The simply supported edge. For the velocity field we have W = 0 at the edge, 
where W denotes the rate of normal deflection of the plate. For the moment field the 
moment M,, , acting in the direction of the normal n to the boundary, vanishes. This 
shows that at the edge either the two principal moments have opposite signs, or at 
least one of them is zero. We are therefore in regime B when the boundary is a principal 
stress trajectory, or in regime BC when the trajectories meet the boundary at an angle. 
In the latter case w = 4M, cos 2a, a being the angle between the direction 1 and the 
outwards-drawn normal. 

b) Clamped edge. Along a clamped edge where deformation of the plate occurs the 
edge is a principal stress trajectory. 

To prove this let us suppose first that no hinge occurs along the edge. If u, v are 
orthogonal parameters such that at the edge u = 0, and if 


Il = Ldw + 2M dudv+ N dv’ (18) 
is the second fundamental form of the surface into which the plate is deforming, we 
conclude from W W, = W,, = 0 along the boundary u = 0, that 

M=N=0. (19) 


This shows that the edge is a line of principal curvature and therefore a line of principal 
stress. In the case L = 0 it follows from the fact that as u — 0, N/L — 0. If a hinge 
occurs, then | 17, | = M, is necessarily a principal moment. If the plate is bent down- 
wards, then in the latter case w = — M,/2. 
c) Free edge. As for the simply supported edge we have the condition » = 
1/,, 2) cos 2a. In addition the shearing force vanishes at the edge, 
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or 


M, 0 ; Oa 
Q,=- sin 2a) = —2w A (20 
. 2 al 7 ~ Ot ) 


where Q, can be calculated from (4), and where 0/dt denotes differentiation along the 
edge. 

5. Discontinuities. The question of discontinuities in the moment field has been 
discussed previously by Prager [9] for the von Mises yield condition. In the following 
we investigate a discontinuity separating two regions which may or may not be in the 
same plastic regime. Figure 5 shows the plane of Mohr’s circles and the image of Tresca’s 
yield condition which is an oval of diameters M, and 2M, . Looking at the possible 
intersection of two circles, we see that three essentially different possibilities of strong 
discontinuities can arise, namely a jump of type “BC — BC’’, which can be BC — BC, 
BC — EF, EF — EF (see notation of Fig. 1), or a jump of type “BC — AB”’, or finally 
a jump of type “BC — B”’. 

a) Jump of type “BC — BC’’. Denoting by M, , M,, the moments acting on an 
element of the line of discontinuity and by a, , a; the angles of the stress trajectories 
on both sides with respect to the normal, we conclude from Fig. 5 that 











Fig. 5 
T 
a; = 9g ~ Ga (Or a, = FT — Gas (21) 


because .V/, , M,, are continuous. When there is deformation of the plate near the discon- 
tinuity and assuming that a, ¥ 0, 7/2, both of the curvature-rates x, and x,, must be 


continuous. This requires 


‘ ‘ 
¢ cos 2 = kK; cos 2a, , 
Ky a; l Xo (22) 


\ 
k, sin 2a; = x; Sin 2a, , ) 


because x; = — kz on both sides of the line. From this and from (21) it follows that 
kK, = k» = 0 on the line of discontinuity which would therefore be a rigid fiber and a line 
of principal curvature. Thus a; , a, = 0, r/2, but then x, , x. can have any value. 
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A moment discontinuity between two regions which are both in the plastic regime 
BC is not permissible, but can take place between BC and EF. 

b) Jump of type “BC — AB?’’. It is sufficient to consider the two cases of a jump 
from BC to AB and the case of a jump from BC to AF, all other cases being similar. 
Let us start with the first one. 

The continuity of M, and M,, requires 


M} cos’ ai + M, sin? a? = Mi + M, sin’ a! | (23) 
f se 


(M, — M;) sin 2a, = My sin 2a%,, ] 
index 7 belonging to AB, index a to BC. The velocity field satisfies the conditions 
1 2 , a , 
Kp COS @i = ky COS 2a 
° | (24) 
kp Sin 2ai = 2x3 sin 2% .) 
From (24) we get, x; , x; supposed non-zero, 
9 , az , 
tg oq tg a; ’ (25) 


2a, = at + mr, 


where m is 0 or 1. Equations (23) are now 


‘ i a ; 
M; cos’ a + My sin’ a, = Mi + 7 (1 + cos a’), | (26) 
2(M, — M;) cosa’, = + M,, 
where the upper sign refers to the case m = 0. Taking either sign we get from (26) 
M? = M,,/2, which is not true. The assumption «x; , «2 * 0, is therefore false. Again we 


conclude that a! , af = 0, r/2. Similar considerations for the angles can be made for 


jumps between BC and AF. 

From regime BC we pass to AB over a line of principal curvature which is in regime 
B. If this line is not straight it must be a line of the set 2 of the stress trajectories and 
M, = M, = 0. The discontinuity is then a weak discontinuity in the terminology of 






Regime AB 


(straight lines) 


Regime BC 


(isometric net) 


o>’ 


Fia. 6. 
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Prager [9]. Figure 6 shows the behavior of the fields when we suppose that regime B 
does not hold throughout a finite region. 
The evolute e is connected with the regime-B-line by the condition 
1 E a 
u (2 ) 


9 7372 
pi I lila 
where £ is the scale factor of the isometric net at the regime-B-line. 
For the moments we have the conditions 


3 5 9 A 0 oO 
M,+Act E oa . | pip dp, — [ pip dp, = 0, o= = : (28) 
Pi Pi o 


d 


and for the normal component of the shear-force to be continuous, 


l . “a ME, 
1 | = | PP: do. | st qi ~ oR?" (29) 


Until now, no non-trivial examples of such a weak discontinuity have been found. 
Further the jumps of type “BC — B” have not yet been studied. 
6. Concentrated loads. The following consideration is not a direct application, 
but shows how useful the discussed fields are for the determination of the limit load. 
Let us take a plate with a convex boundary which is simply supported at its edge and 


which is loaded by a concentrated force (Fig. 7). 


» 





Fia. 7. 
4 
a) In order to get a statically admissible stress field [3] let us cut out a small circular 
region of radius e around the force and consider the isometric net indicated in the figure. 


The corresponding isometric coordinates u, v are chosen such that v = g = const. are 
straight lines, and u = log r, r being a measure of distance from the point of application 
of the load P. At the simply supported edge w = 3M, cos 2a, and at the small 
circle w = 4M, . If we assume that the plate is in regime BC then Eq. (14) becomes 
Aw = 0. (30) 
As — M,/2 < w < M,/2 on both boundaries, we conclude from a principal theorem 


on harmonic functions (Courant and Hilbert [10]) that w satisfies the inequality 
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Mo My 
—“— 2°35 (31) 


in the region between the boundary and the small circle. The stress point therefore lies 
between the two points B, C on the yield curve, and as the harmonic function w exists 


we have found a statically admissible stress field. 
Let us now calculate the shearing force across the small circle. From (4) we get 


"tna ell, oat, Pe [ E'? Q, dv = [ sio-on,, 
0 70 
As w, is of the order 1/log ¢ the limit of this equation, when e — 0, will be 
P = 2xM, , (33) 


which is a lower bound for the limit load. 
b) In order to find a kinematically admissible velocity field we take the isometric 
net of Fig. 8 which is given by the analytic function 


z= fe). (34) 





The function f is such that the transformation ¢ = f ‘(z) maps the boundary of the 
plate into the unit circle and the point of application of the load into the origin ¢ = 0. 
The parameter u is then zero on the boundary and negative inside, and u has a log- 
arithmie singularity at the loaded point. The velocity field W = — uW, satisfies the 
boundary conditions and the conditions of regime BC and is therefore a kinematically 
admissible velocity field. As the rate of work done by the concentrated load would be 
infinite, we replace the area inside the small “circle” u = uw. = const. by a flat area, 
for which the rate of dissipation is zero. It remains to calculate the dissipation outside 


this ‘circle’ and at the hinge on u = wu . As 


+My? _ -We +e)” 
xk; = E = E ’ 
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(this net does not represent the lines of curvature), 





227 »0 
ao M.Wo] dw] ite)? du 
a ia (35 
p du 
At the hinge Fig. 9 shows that 
\ 
Fic. 9 
; —dW —VW W 
tg i ds, E 2 ] 2 
As As E’’* 
R a 
tg 7 W 
W As oe : 
Decne. = ———5- on 27M ,W, . (36) 
; As i 
The second fundamental theorem of limit analysis gives therefore 
. 24M ,W MW, f°" a . 
PW, < 2xM,W, + 4+ o] dv | y, | du, (37 
u u 0 Ju 
and in the limit 
P< 2rM,, (38 


which is an uppel bound on the limit load. The exact value of the limit load is therefore 
27M 

We remark that this result is also true for a concentrated load on a clamped plate. 
The same statically admissible stress field can be used, and for the velocity field we 
must only take into account the dissipation at the hinge produced at the clamped edge. 
This gives in Eq. (37) an additional term, which is of order 1/| w%) | and which is zero in 
the limit. However as R. M. Haythornthwaite pointed out, the statement can be proved 
in this case in a more elementary way. 

7. Saint-Venant’s principle in plasticity. The application of the concentrated load 
is only of theoretical value, because in practice the load will always be distributed over 
a small area. However the examples of circular plates given by Prager and Hopkins show 
that the following conjecture may be true: 

If the diameter of the area, where the load acts, is small in comparison with the 
distance from the edge, the limit load will be near the value 27M, , and 27M, is the lowest 


limit load for the plate. 
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Matrix calculus. By E. Bodewig. North-Holland Publishing Company, Amsterdam; 
Interscience Publishers, Inc., New York, 1956. xi + 334 pp. $7.50. 


The author feels that the usual matrix notation does not exhibit “the true building blocks of a 
matrix, the rows and columns’’, and therefore causes ‘discrepancy between thought and calculation, 
and lack of elegance”. To remedy this situation, the author introduces the following notations: A;, = 7-th 

{1 ; = j-th column of matrix A; e; = column vector whose 7-th element is unity while all others 


row, and z 
vanish; E;; = matrix having the element at the intersection of row 7 and column j equal to ind all 
other elements equal to zero. Thus, A.; = Ae; and A;, = e;A, where the prime indicates the transpose; 
a;; = e,Ae;. These notations enable the author to represent operations explicitly that otherwise would 
have to be described verbally. For instance, A(J + cE;;) is the operation that adds the c-fold 


i-th column and leaves all other columns unchanged. 


of the i-th 


column of A to the 
The book consists of four parts, the first of which provides an excellent survey of the general theory 

(85 pp.). The remaining parts are concerned with specific problems of computation, each p ontaining 

sections on direct and iterative methods. Part II (96 pp.) is devoted to Linear Equations, Part III 

(48 pp.) to Inversion of Matrices, and Part III (98 pp. to Eigenproblems). 

While often highly concise, the exposition is always clear and easy to follow, though the English 


is occasionally awkward. Examples of a new approach to a familiar problem abound throughout the 


book. In this connection, Section III, C on Geodetic Matrices merits particular mention. The author’s 
frequent remarks on the suitability of numerical methods for different types of computing equipment 
will be welcome 1 by many readers. 


W. PRaGER 


Nuclear reactor physics. By Raymond L. Murray. Prentice-Hall, Ine., Englewood Cliffs, 
N. J., 1957. xi + 317 pp. $10.00. 


This book is a very cle irly written text on the fundament il features of reactors with emphasis on 
temperature 


tron flux distribution, on transient effects in reactors, and on the ny 
nd the slowing 





neutron physics and n 


effects on reactor operation Included are chapters on reactor control, transport theory, :‘ 
down of neutrons 
The author has ided sufficient mathematical treatment to make the discussior iseful and 
interesting but not so n h that only the specialist would be interested. The book seems have ex- 
cellent balance between what 
tician might want to know without becoming a specialist. The book seer 


a student of reactor fundamentals might need and what a well trained 
physicist or mathem ns excellent 
for those who simply want to know what problems enter into reactor design. 

{0HN TRUELL 


The International Dictionary of Physics and Electronics. D. Van Nostrand Co., Inc., 

Toronto, New York, London, 1956. xvii + 1004 pp. $20.00. 

This book is an ambitious undertaking (even for fifteen editors) amounting to one thousand pages 
of laws, equations, definitions, principles and concepts of physics. The objective as stated by the editors 
has been to provide a book useful as a general reference in physics. This work will undoubtedly be useful 
to many people especially to those who seek information outside their field of specialization 

{0HN TRUELL 
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ON FINITE EXPANSION OF A HOLE IN A THIN INFINITE PLATE* 
BY 
P. G. HODGE, JR.t anp R. SANKARANARAYANAN 
Polytechnic Institute of Brooklyn 


Summary. A complete solution for all values of applied pressure is obtained for the 
expansion of a hole of zero radius in an initially uniform infinite sheet. The analysis is 
compared with those of previous investigators and found to be simpler and/or more 
complete. It is shown that the results are applicable to the expansion of a finite hole 
in a uniform sheet and to the expansion of a hole in a tapered sheet. 

1. Introduction. The problem of finite expansion of a circular hole in a thin infinite 
plate has been discussed repeatedly in the literature. Taylor [3]' treated the problem 
by using Tresca’s yield condition and the flow rule associated with Mises’ yield con- 
dition, but he considered a rigid perfectly plastic material. Hill’s [4] extension of Taylor’s 
work was also restricted to a non-strain hardening material. Prager [1] solved the problem 
for a rigid strain hardening material which satisfies Tresca’s initial yield condition and 
the associated flow rule, but his solution is restricted to a certain finite value of the 
pressure. Alexander and Ford [5] have treated the same problem, using Mises’ yield 
condition and the associated Prandtl-Reuss relations. 

The present paper deals with the same problem treated by Prager, but the solution 
is obtained for all values of the internal pressure. As in Hill’s treatment [4] the math- 
ematical analysis in the plastic regions is simplified by reducing the problem to the 
solution of ordinary differential equations by means of a substitution of independent 
variables 

Specifically, a hole is expanded from zero radius in a thin infinite sheet of work- 
hardening material which satisfies Tresca’s yield condition and the associated flow 
rule. An incompressible material is envisaged, and since finite plastic deformations 
are considered, elastic deformations are neglected. 

In the case of plane stress problems Tresca’s yield condition can be represented 
by a hexagon in the o, , o2 plane (Fig. 1). If the “stress point” with coordinates o, , o2 
remains in the interior of the hexagon, the material is rigid. For a point on the boundary, 
plastic flow may commence in accord with the plastic potential flow law [6]. According 
to this law, on a side of the hexagon the strain rate vector must be the unique normal, 
while at a corner it may assume any position between the limiting normals. 

For a perfectly plastic material the hexagon in Fig. 1 is fixed, and the stress point 
is not permitted to leave it. For the isotropic strain hardening material considered in 
this paper, the yield curve maintains its shape, center and orientation, but merely 
changes its size. 

It has been shown [7, 8] that under certain restrictive conditions the plastic flow 
laws can be explicitly integrated. These conditions, defined as a “regular progression”’ 

*Received March 29, 1957. The results presented in this paper were obtained in the course of re- 
search sponsored by the Office of Naval Research. 


tNow at Illinois Institute of Technology. 
iNumbers in square brackets refer to the references collected at the end of the paper. 
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Fic. 1. Tresca hexagon lor plane stress. 


of the stress point, state that it must never lose contact with a side of the hexagon. 
In other words, the stress point may stay on a side, stay in a corner, or move from a 
side to a corner; it may not move from one side to another, from a corner to a side, or 
back into the interior. 

Prager [1] considered the enlargement of a circular hole of radius a, to 1.48a, by 
the application of a gradually increasing uniform pressure p to the edge of the hole. 
The stress point of any particle starts at the origin for p = 0 and moves along the straight 
line OG as p increases (Fig. 2). It first becomes plastic at G and moves along GD to D. 
Next, the stress point remains a finite time at D, engaging in plastic flow. Finally, it 
moves back along DC to the accompaniment of further plastic flow. However, Prager 
did not consider the final stage of plastic flow, when “non-regular progression’’ of the 
stress point takes place from the corner D, back to the side CD of the yield hexagon. 
Therefore, as he pointed out, his solution is only valid for a/a) < 1.48. In particular 
Prager’s solution is not applicable to the expansion of a hole from zero radius 

Since the stress point progresses non-regularly in the final stage of plastic deformation, 
the explicit integrated forms of the flow law cannot be used. A similar condition was 
encountered in an earlier paper [2] concerned with the bending of a simply supported 
annular plate. There it was found that the complete solution was built up of zones and 
that appropriate initial conditions were always available from plastic zones already 
formed. The analysis of the present paper reveals that a similar condition holds in the 
finite expansion of a hole in an infinite sheet. 

For conceptual simplicity, it is first assumed that the initial radius of the hole is 
zero. Since the plate is infinite, this implies that there is no characteristic length in the 
plane of the sheet. It follows that the geometric solution at all times must be similar to 
itself. Therefore, rather than the solution depending upon time and space independently, 
it depends only upon a conveniently chosen space-time ratio. This technique was first 
used by Hill [4]. 

The solution for a hole expanded from a finite radius is obtained from that for ex- 
panding a hole from zero radius simply by discarding the part of the solution which 
is not required. This is so because it is immaterial whether the pressure at any radius 
is applied by an external agency or through the displacement of an inner annulus of the 
sheet. The stress and velocity in an element depend only on what happens beyond the 


radius. 
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The basic theory is presented in the next section. Section 3 gives the details of the 
solution. The paper will conclude with a discussion of some limitations and extensions 
of the analysis, as well as its relation to other investigations cited. 

2. Basic theory. The state of a thin infinite sheet of initially uniform thickness is 
fully defined by the principal stresses o, and oy , the radial displacement u, and the 
thickness h. The stresses must satisfy the equation of radial equilibrium 


d(ha-) 4 h(o, — ov») i (1) 
or r 


Where the material is rigid, it may be regarded as the limiting case of an elastic 
material as Young’s modulus becomes infinite. Therefore, in a rigid region the stresses 
must also satisfy the compatibility equation for an incompressible material 


(922 = ; $e) + 3 (os — §,) = 0. (2) 


\ or or 2 


A stress point in the interior of the hexagon must satisfy Eqs. (1) and (2). The flow law 
in a plastic region is expressed in terms of the plastic strain rates 


€ ov /or, €éy9 = 0/1, e. = (1/h)(0h/dt), (3) 
where v is the radial velocity. Since the material is incompressible, 
ov/dr + v/r + A(log h)/at = 0. (4) 


The plastic flow law assumes different forms for each side or corner of the yield 
hexagon (Fig. 1). For the type of loading considered here, the relevant plastic regimes 
are CD and D. On CD, the strain rate vector must be directed normal to the side, hence 


+e = —e, = 6. (5) 
Further, the stresses must satisfy 
¢, = 0, o5 > 0 (6) 


if the stress point is to remain on CD. 
In the corner D, the stress point must always lie on the o, axis; hence 


es = 0. (7) 


% 





Fic. 2. The stress profile. 
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The strain rate vector must lie between the values corresponding to the sides CD and 
DE; hence 
0 < €s < —€, o (8) 
It remains to consider the rate of hardening. As in Prager’s treatment [1] this will 
be assumed to be proportional to the energy dissipation. Thus 
dc/dt = aD, (9) 
where a is a constant, D the dissipation function, and 
o = max || g, |, | os |, lo, — oe |]. (10) 


On the side CD, as well as the corner D, Eq. (9) can be written 


da/dt = —age, , 
where 
c= 03 —- 70, (11) 
At r = ©, the radial stress must vanish, and at r = a, the current hole radius, it 
must be in equilibrium with the applied load. Thus 
a f= @, ho, = —97; (12) 
at r= o, o, = 0. (13) 


Finally, at the beginning of the expansion the thickness must be equal to its initial 


constant value 
at t=0, h=ho. 
3. Solution. When a hole is expanded from zero radius, the material at a sufficiently 
large distance from the hole is rigid because it is stressed below the yield limit. Therefore, 


h = h, and u = 0. Atr = ~, the stresses must vanish, while at the elastic-plastic bound- 
ary they must satisfy 

o¢,— 03 = —O.- (14) 
Therefore, Eqs. (1) and (2) determine the stresses to be 

= s —— <=, (15) 


—_ oe oe 
where a, denotes the yield stress of the virgin material and p is the radius of the elastic 
plastic interface. The solution given by Eqs. (15) is valid for r > p. 

For r somewhat less than p, the stress point will be in regime CD. For this portion 
Eq. (5) shows that h is still equal to h, . Since there is no thickening of the material, 
there can be no straining and hence it is still rigid, even though plastic. The stresses are 


now given by Eqs. (1) and (14), with Eq. (15) as a boundary condition at r = p. Thus 
l p 
o, = me + log =} ; 
< T (12% 
(16) 


d 
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It is clear from the second Eq. (16) that the hoop stress vanishes for 
r= pe”. (17) 


Therefore, it follows from inequalities (6) that Eqs. (16) are valid for nu < r < p, where 
p and yu are related by 


p = pe”. (18) 


The portion of the sheet just inside the circle r = yu corresponds to the regime D. 
It is evident that u will increase with the pressure, and hence it may be conveniently 


taken as the time parameter. The ratio of the original radius r» of a particle to the radius 
of this boundary is then taken as the space-time ratio mentioned in the Introduction. 


Thus the entire solution will be a function of 
L = To/ p- (19a) 


The equilibrium equation, the incompressibility relation, the law of strain hardening, 
and the flow rule may now be expressed in terms of this parameter x. The resulting 
equations will be ordinary differential equations, and the solutions will be obtained 
relatively easily. It is convenient to define the following dimensionless quantities 


£ r/ wu, 8. = o-/Ge; 83 = o9/00; 8s = a/a, 9=h/ho, (19b) 


all as function of x only. If ¢ is any function of x alone, then the following relations are 


seen to hold: 





dg/dro = (1/u)y’, 
dg/dp = —(1/u)zy’, 
, a(r/ 20 
dr/dro _ ud(r/u) = 4 3 = i; ( ) 
Oro Oro 
‘ dp/ATo 
) /or = at a 4/,,e 
ww arjar, Me» 


where primes denote differentiation with respect to x. Therefore, the velocity and strain 


rates are 
v = Or/du = = — zt’, (21) 
e, = du/dr = —xt’’/nt’, (22) 
€9 = v/r = (& — xt’)/uk, (23) 
= Eo = — en /un (24) 
At the boundary of the plastic region x = 1, the displacement and velocity must 


vanish, the thickness must equal its initial value, and the stress point must be at D. 
It follows that, at x = 1, 


gE =], = 1, 7= 1, s= 1. (25) 
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a | 
~ 


The strain hardening law and the incompressibility relation, Eqs. (10) and (5), 


respectively, may now be written as 


sir -7 
—IS& ANE g; 
—xé"’/E + 1 — xt’ /E — xn’/n Q. (26) 
Integrals of Eqs. (26) which satisfy Eqs. (25) are then found in the form 
we" l (27) 
gE'n w. (28) 
Finally, the equilibrium equation in the corner D becomes 
(Ens? ¢’ = (. (29) 
The solution of Eqs. (27), (28) and (29) satisfying the boundary conditions (25) are 
g — © = i 4 S 0 
2+a)/U+a ‘ 
;*2 |] + (1+ a)zx |, (30) 
2rT a 
(2 se QA 
uy] x (1 x 


This solutign will be valid so long as Inequality (8) is satisfied. It is readily verified 
that this is true for all 1/2 < & < 1. Substitution of Eqs. (19) into Eqs. (30) yields the 
results previously obtained by Prager [1]. 

For é < 1/2, the solutions given by Eqs. (32) to (35) predict a strain rate vector 
which points too far up, suggesting that the portion of the plate just inside the circle 
& = 1/2, corresponds to plastic regime CD. 

For the solution in the region  < 1/2, Eqs. (27) and (28) are still valid since Eqs. 


(10) and (5) hold in regimes D and CD. In addition Eq. (5) implies that 


= 2 = Sais ( 
— al | 31) 


Here, as in the following, integration constants are determined to give continuity with 


Eqs. (30) at x = 1/2 
3ecause 7 is a constant, Eq. (28) is easily solved, the result being 


| ee | | (32) 
. E 1(1 + a) ; mn 


With é known, Eq. (27) furnishes 


x l ite 
s = (n,/x)* 4+ — ; (3 
m/*) E ' 4A(L + =| pe) 


az r* | r 2 . a a 2 x) 
q 1) + dx — - (34 
J, @/s EB 41 4 5 E | ) 
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Fic. 3. Distribution of radial and circumferential stresses. 
Finally, since s) = s + s, , the value of sy is obtained by adding Eqs. (33) and (34). 


It is easily verified that Inequalities (6) are always valid for 0 < ~ < 1/2. Therefore, 
the solutions represented by Eqs. (33) and (34) give the complete stress distribution 
in the plastic region. It is further verified that s increases as the value of £ is decreased 
from 0.5, so that the strain hardening continues. 

4. Conclusions. The complete solution has been obtained for the finite expansion 
of a circular hole from zero radius in a thin infinite plate of initially uniform thickness. 
The material of the plate is rigid-plastic, satisfies Tresca’s yield condition and the 
associated flow rule, and hardens isotropically. 

As already stated, the solution for a hole expanded from a finite radius is obtained 
from the present analysis simply by discarding the part of the solution around the edge 
of the hole which is not required. This is a consequence of the fact that the outer annular 
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Fic. 4. Thickness variation in deformed sheet. 
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fibres of the plate receive radial pressure transmitted only from the inner fibres, and the 
inner fibres may be equally well replaced by some externally applied pressure. 

As pointed out by Alexander and Ford [5], the governing equations in the plastic 
region of a plate of initially uniform thickness are identical with those of a plate which 
varies in thickness proportionately with the radius. Therefore, the plastic region analysis 
of the paper applies equally well to such a plate. 

The great simplicity achieved by using Tresca’s yield condition and the associated 
flow rule will be obvious by a comparison of the present analysis with that of Alexander 
yield condition and the 


, 


and Ford [5], who treated the same problem by using Mises 
associated Prandtl-Reuss relations. Therefore, it is of some interest to compare the 
results of the two theories. Figures (3) and (4) show the stress and thickness distri- 
butions, plotted as functions of r/p, for two values a = 5.0 and a = 0.5, of the hardening. 

Figures (3) and (4) show that there is fairly good agreement between the two distri- 
butions of stress, though there is some discrepancy between the two distributions of 
strain. Experiment alone, of course, can decide which analysis predicts the real behavior 
of the material more closely. However, even if experimental evidence were to show 
the Mises-type solution to be more accurate, the present treatment might still be valuable 
as an approximate solution due to its much greater simplicity. 
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CROCCO’S VORTICITY LAW IN A NON-UNIFORM MATERIAL* 


By C.S. WU ann W. D. HAYES (Princeton University) 


In the derivation of Crocco’s vorticity law [1] there is a basic restriction, that the 
material in the flow field be uniform, that the flow be isocompositional. This restriction 
is not generally appreciated, and is only implied in most published derivations of this 
law. The purpose of this note is to express the corresponding form of the vorticity law 
valid with non-uniform material, such as a gas mixture with non-zero concentration 
gradients or a material in which a chemical reaction is proceeding. 

The basic assumptions which are needed are; (i) no body forces except gravity; 

ii) no frictional stresses; (iii) steady flow; (iv) thermodynamic quasi-equilibrium; and 
(v) continuous first derivatives of thermodynamic variables and velocity. With assump- 
tion (iv) the material obeys a differential state equation of the form 


dh = T ds + 2 + i by, Un , (1) 


where u, and-n, are the molar chemical potential and the specific molar concentration, 
respectively, of the kth constituent, p and T are the pressure and absolute temperature, 
and h, s, and 1/p are the specific enthalpy, entropy and volume. 

Since the material everywhere in the flow field obeys Eq. (1), we may write the 
corresponding gradient equation ‘ 


Vh =TVs + == + SmVn . (2) 
With assumptions (i), (ii) and (iii) we may write the equation of motion of the fluid 
vi-ax(vx@ + V+ va =0, (3) 
where Q is the gravitational potential. The total enthalpy A, is defined 
h=h+ r + 2. (4) 
Eqs. (2), (3), and (4) may be combined to give the desired vorticity law 
qX(V Xq) = Vio — TVs — Do mVn, . (5) 


If there is diffusion present the velocity is defined as a mass mean, as the total mass 
flow vector divided by the total density. 

In the absence of any real-fluid effects in an isocompositional flow, both h, and s are 
constant along streamlines, and Vio and Vs are both directed normal to q. These real- 
fluid effects are viscosity, excluded by assumption (ii), heat conduction, and relaxation. 


*Received January 10, 1957. This work was partially supported by the Office of Scientific Research, 
Air Research and Development Command, USAF, under Contract AF 18(600)-498. 
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With a non-uniform material >> 4, dn, is zero along a streamline if there is no diffusion 
v chemical reaction is reversible. Thus, since diffusion and irreversibility of 


and if 
real-fluid effects associated with varying composition, no new 


chemical reaction are 
feature appears in the generalized form of Crocco’s vorticity law. 

The vorticity law is most useful if h) and the n,’s are constant over the entire flow 
field; in this ease the entropy gradient may be obtained directly from the vorticity 
or the lateral components of the vorticity from the entropy gradient. However, with 
no real-fluid effects present and with sufficient knowledge of the flow field the vorticity 
law still gives a useful relation between vorticity and entropy gradient, as then h, is a 
only of the streamline and the n,’s may in principle be obtained, and the lateral 


functior 
gradients of these quantities may be calculated. 
In a binary mixture we may set 
c nM, , (Ga) 
l—c NoMa ; 6b) 


Pe. (6c) 


where 1/7, and J, are 
this notation we have 


2. an p de (7) 


and Crocco’s vorticitv law takes the form 


qx(V Xq Vho — TVs — pVe. (8) 
(EFERENCE 
[1] L. C1 0, E n neue olro funktion fri die E forsch ing de r Be wequng de r Gase mit Rotation, Z ANIM 
17, 1-7 (1937 


THE SCHWARZIAN DERIVATIVE AND THE APPROXIMATION 
METHOD OF BRILLOUIN* 


By AUREL WINTNER (The Johns Hopkins University) 


1. Let the coefficient function of the differential equation 
2’ + P(x = 0 (1) 


be the square of a positive function f(t) which is given, and possesses a continuous 


second derivative, for large positive ¢, say for f, S t < . Then the function 


g(t) = | f(s) ds (2) 


J to 
is strictly increasing and has a continuous third derivative, whereas the function 


*Received January 17, 1957. 
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€ 


4F(t) = 3f°()/f'°'(®) — 2f''(O/f* (0 (3) 


exists and is continuous. 
In view of (2), the definition of F(é) is equivalent to 


F(t) = —2f7*(O{e(d}, (4) 


if | } is the symbol of the differential parameter of Lagrange-Schwarz, 
, p99 2y 
le} = ’’/¢’)’ — &"/¢') /2. 
It is also seen from (3) that 


F() f(t) = k(Ok’’(0), where k(t) = [f()]'”’. (6) 


2. If t > ~, then the function (2), being increasing, tends either to a finite limit 
¢ or to ¢f = o. Suppose that the latter is the case, 


| sat = @ (f > 0), (7) 


. 


and that the function (6) is absolutely integrable, 


. 


| | F(t) | f(i) dt < @ (f > 0). (8) 


The following theorem was proved ten years ago’ (and, in certain particular cases of it, 
was repeatedly rediscovered” in the meantime): 

For large positive t, let {(t) be a positive, twice continuously differentiable function satis- 
and (8), where F(t) is defined by (3). Then the general solution of the 
(1) 7s the superposition of two linearly independent solutions which 

) to the real and imaginary parts of the function 


at 
y(t) = [f()]'” exp E | f(s) as| (9) 
“to 


In addition, differentiation of this asymptotic representation ts allowed. 
words, the assumptions (7) and (8) are sufficient in order that the general 


lying conditions (7 
differential equation 


are asymptott ast 


In other 
solution of the differential equation (1) be of the form 2(t) = ¢,2,(t) + e.2%,(t), where, 


f > a 


as 
a(t)~ yl), x2(t))~ y() and al) ~ yd, «20 ~ yd, (10) 

if the real functions y,(t), y2(t) are defined by placing 
y(t) + iy(t) = [f()]-'”? exp [ig(d]). (11) 


3. The point in the theorem just quoted was the generality of its assumptions (and, 
correspondingly, the comparative simplicity of its proof), since, under assumptions 


1A, Wintner, Phys. Rev. 72, 516-517 (1947). 

*See, e.g., M. Zlimal, Czechoslov. Math. J. 6 (81), 75-93 (1956) (in German, with a Russian abstract), 
vhere one of the books of R. Bellman (1953) is mentioned. The main result of Zlimal’s Chap. II (Theorem 
6 on p. 83, along with the proof on pp. 84-85) is the simplest of the corollaries of the note in the Phys. 
Rev., referred to in the preceding footnote; see the general theorem italicized above, which shows that 


the problem is independent of Zl4mal’s convexity assumption. 
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more stringent than the pair of conditions (7) and (8), the asymptotic form (9) was 
well known from the mathematical as well as from the physical literature of the subject.° 
The proof of the theorem depended on a combination of two steps: (a) the standard 
use of the formal substitution of the Liouville‘ and Riemann,’ and (8) the application 
of an “Abelian lemma,” which runs as follows.° 

If « = e(¢) is defined, for large positive y, as an arbitrary continuous (possibly not 
even real-valued) function which is absolutely integrable, 


[ | ee) | de < o, (12) 
then the general solution z = 2(¢) of the differential equation 
d’z/dy? + [1 + e(y)]ze = 0 (13) 


s the derivative, dz(¢)/dy, of every solution) has, as g — ~, the same asymp- 


(as well < 
vior as what results in the trivial case, e(g) = 0, of (13); so that 


totic beh: 
t+ ¢, cos ¢) > 0, dz(y)/de — (ce, cos gy — ce, sin g) > 0 (14) 


2) — (c, sin g 4 


W\} 


as ¢ — ©. In the application of this Abelian lemma, the independent variable g must 
be chosen to be the function defined by (2) or (4) above. The réle of step (a) was pre- 
cisely the reduction of (1) to (13), where e(¢) is “small” (for large g) in an appropriate 
sense, whereas step (@) consisted in assuring that condition (12) is a sufficient “degree 
of smallness’ for the asymptotic form (14) of the general solution of (13). 

In what follows, a more general and less artificial approach will be taken to the 
asymptotic behavior of the solutions x(¢) of (1) and of their derivatives. The improve- 


ments will depend on a formal fact pointed out by Brillouin.’ 





, Sansone, Equazioni differenziali nel campo reale vol. 2, Bologna, 1949, p. 53; 

see also pp. 101-110 of Titchmarsh’s book, referred to in footnote’. 

4See G. Sansone, op. cit., vol. 2, pp. 68-70. 

5B. Riemann, Gesammalie mathematische Werke, 2nd ed., Leipzig, 1892, pp. 335-336, or the presen- 
tation in E. G. C. Poole’s Theory of linear differential equations, Oxford, 1936, pp. 121-122; see also E. C. 
Titchmarsh, Eigenfunction expansions associated with second-order differential equations, Oxford, 1946, 
p. 102. 
64 proof, based on successiive approximations, was given by H. Weyl, Gruppentheorie und Quanten- 
mechanik, Leipzig, 1928, pp. 61-63. Actually, as pointed out in my note in the Phys. Rev., referred to 
same proof is contained already in an important paper of M. Bécher (1900). For another 


above!, the s 
arrangement of the proof, a proof based on the “Abelian” aspect, see A. Wintner, Amer. J. Math. 69, 
251-272 (1948 


4 


As shown recently [A. Wintner, Amer. J. Math. 78, 895-897 (1956)], the latter proof is such as to 
lead to simple criteria which are substantially more general than (12) and are, nevertheless, sufficient 
in order that (14) should hold for the general solution of (13). 

As pointed out at the end of the paper just mentioned, the question of asymptotic behavior is stricter 
than that of stability (boundedness). For an application of the Loiuville-Riemann transformation to 
uestion, see U. Barbuti, Annali di Pisa, ser. 3, 8, 81-91 (1954), where, with regard to the 
bility, recent papers of G. Ascoli (1947), L. A. Gusarov (1949) and U. Barbuti (1952 


For references, see G 


+ 


the latter 
question of st 


are cited. 
7L. Brillouin, Quart. Appl. Math. 6, 166-178 (1948) and 7, 363-380 (1950) (in particular p. 172 and 


p. 364 respectively); see also L. Brillouin and M. Parodi, Propagation des ondes dans les milieux pério- 
diques, Paris, 1956, p.216. Strictly speaking, in what Brillouin is interested is, not the asymptotic behavior, 
for indefinitely large ¢, of the solutions z(t) of a fixed differential equation (1), but a situation in which 
the t-range is limited but the differential equation (1), instead of being fixed, contains a variable parameter. 
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4. For large positive t, let f(t) be any, perhaps not even real-valued, twice con- 
tinuously differentiable function which is nowhere zero, 


f® #0 for tSt< (15) 


(which, if f(¢) is real-valued or purely imaginary, means that 
fi) >0 or f(t) <0 (16) 


throughout]. Then, with either of the continuous determinations of the square root 
occurring in it, Eq. (9) defines a twice continuously differentiable function 


y(t) #0 for tSt< om. (17) 
Two differentiations of the definition (9) show that 
yO) + 7°) — POFOlyO = 0 (18) 


is an identity, if F(t) is defined by (3) [that is, by (4) and (2) together or, equivalently, 
by (6)]. Starting with the identity (18), Brillouin observes that the function (9) can be 
expected to furnish an approximation to a solution z(t) of (1) when the coefficient, 
[ ], of y(t) in (18) is “sufficiently close” to the coefficient, f’(t), of x(t) in (1), that is, 
if f°(t)F(t) is “sufficiently small.’”* 

The issue is, of course, the specification of conditions under which such a “degree 
of smallness”’ (for large ¢) will actually be realized. But it turns out that such a specifica- 
tion can be obtained, in explicit terms, from a general criterion. The criterion in question 
is as follows.” 

If y /,(t) and y = y,(t) are two, linearly independent, solutions of a differential 


equation 
y’’ + gdy = 0 (19) 


in which g(t), where f, S t < ~, is any continuous function, and if h(é), wheretf, St < @, 
is any continuous function having the property that the corresponding differential 


equation 
x’ +h()x = 0 (20) 


is so ‘‘close’’ (as ti — ) to the differential equation (19) that 
| du P+] P19 — WO lat < @ (21) 


holds, then the two differential equations are asymptotically equivalent, in the following 
sense. The differential equation (20) will possess two linearly independent solutions, 
x = x,(t) and x = 2,(t), which, along with their derivatives, will satisfy the asymptotic 
relations (10) ast— o. 

5. In order to combine this criterion with Brillouin’s method, let (18) and (1) be 
identified with (19) and (20) respectively (so that | g — h| becomes | f*F |) and, under 
the assumption that f(¢) is positive, let y,(t) and y,(t) be those (real-valued) solutions 


8In this regard, see the comments made at the end of the preceding footnote. 
*A, Wintner, Quart. Appl. Math. 15, 428-430 (1957). 








86 NOTES [Vol. XVI, No. 1 


of (18) for which y,(¢) + dy.(t) becomes the (complex-valued) function (9). Then the 


condition (21) reduces to the convergence of the integral 


[ | y(t) |° f°) | FD | dt, 


e 


since | 7; + | Yo y|-and|g—h F | f’. But the convergence of this integral 

is equivalent to the condition (8), since | y |" f’ = f, by the definition, (9), of y, where 
exp | = 1. 

Accordingly, the general criterion of Sec. 4 supplies, as a corollary, not only the 


theorem quoted in Sec. 2, but also an improvement of that theorem. In fact, it is now 
seen that the assumption (7) of the theorem italicized above can be omitted, that is, only 
the other assumption, (8), of the theorem is needed. 

6. This reduction of the assumptions of the theorem could have been proved directly. 
The main point was, however, the possibility of a straightforward approach, one replacing 
the formalism of the Liouville-Riemann substitution by its simplified version, the formal 
identity of Brillouin. 

The spurious character of the restriction (7) is indicated even by Brillouin’s elemen 
tary case. In fact the latter results by assigning the following condition for the coefficient 
function of (1 


Fit) =O tor 4< t < @. 22 


In view of the definition (3), the assignment (22) represents a differential equation of 
second order for f(#), and Brillouin’s approximation theory (loc. cit.) depends on the 


circumstance that the general solution of the (non-linear) differential equation (22) is 


f(b C,/(t + C2 sf 23) 


(C, and C, are the integration constants), as is readily verified by direct substitution. 
Clearly, (22) means that the differential equation (18) is not only ‘‘close’’ to, but is 


identical with, the “perturbed” differential equation (1) and, correspondingly, condition 
(8) is certainly satisfied. But condition (7 happens to be violated. In facet, it vithout 


| 


loss of generality) ¢ l and C, 0, then the general solution, (23), 


f(t) = €°, and so 


of (2?) reduces to 


| f(t) dt | t 7? dt < t, > 0) 


holds, instead of (7). The exact solution of this case of (1) is supplied by the real and 


imaginary parts of the function 
t exp (i/t), 0<t, t < o), 


since (9) reduces to this function for the choice f(t) = ¢ (if 7 1s replaced by Zand a 


constant factor, introduced by the arbitrary integration limit t, > 0, is disregarded). 
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ON A ROUTING PROBLEM* 
By RICHARD BELLMAN (The RAND Corporation) 


Summary. Given a set of N cities, with every two linked by a road, and the times 
required to traverse these roads, we wish to determine the path from one given city to 
another given city which minimizes the travel time. The times are not directly pro- 
portional to the distances due to varying quality of roads and varying quantities of 
traffic 

The functional equation technique of dynamic programming, combined with approxi- 
mation in policy space, yields an iterative algorithm which converges after at most 

\ 1) iterations. 

1. Introduction. The problem we wish to treat is a combinatorial one involving 
the determination of an optimal route from one point to another. These problems are 
usually difficult when we allow a continuum, and when we admit only a discrete set 
of paths, as we shall do below, they are notoriously so. 

The purpose of this paper is to show that the functional equation technique of 
dynamic programming, [1, 2], combined with the concept of approximation in policy 
space, vields a method of successive approximations which is readily accessible to either 
hand or machine computation for problems of realistic magnitude. The method is dis- 
tinguished by the fact that it is a method of exhaustion, i.e. it converges after a finite 
number of iterations, bounded in advance. 

2. Formulation. Consider a set of N cities, numbered in some arbitrary fashion 
from 1 to V, with every two linked by a direct road. The time required to travel from 
i to j is not directly proportional to the distance between 7 and j, due to road conditions 
and traftic. Given the matrix 7 = (t,;), not necessarily symmetric, where ¢;; is the 
time required to travel from 7 to j, we wish to trace a path between 1 and N which 
consumes minimum time. 

Since there are only a finite number of paths available, the problem reduces to 
choosing the smallest from a finite set of numbers. This direct, or enumerative, approach 
is impossible to execute, however, for values of N of the order of magnitude of 20. 

We shall construct a search technique which greatly reduces the time required to 
find minimal paths. 

3. Functional equation approach. Let us now introduce a dynamic programming 
approach. Let 

{, = the time required to travel from 7 to N, 7 = 1,2, ---,N — 1, 
using an optimal policy, (3.1) 

with f 0. 

Employing the principle of optimality, we see that the f; satisfy the nonlinear system 
ol equations 


f; = Min [t+ f,), t=1,2,°---,N—1, 
jeta (3.2) 


fy = 0. 


*Received January 30, 1957. 
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This system differs from the usual systems encountered in dynamic programming 
in that we do not have a ready computational scheme. 
4. Uniqueness. Let us show that there exists at most one solution of the system 


in (3.2). 


Assume that {f;} and {F;} are two solutions, with fy = Fy = 0, and let & be an 
index for which f, — F, achieves it maximum. Then 
fe = Min [t,; + fi] (4.1) 
ik 
F, = Min [t,; + F;]. 
ixk 
Let the minimum in the first equation be assumed for j = r, and the second for j = s. 
It is clear, since ¢;; > 0 for all 7, j, that r # k, s # k. Then we have the equalities and 
inequalities: 
f lian T Te Ss le T 
— f (4.2) 
I t F,<t, +4 
These lead to 
f. —-F,.<f, —F 
ik OS mm 4 s (4.3) 
ohm Fs 
Since k was an index where f, — F, achieved its maximum, we must have 
n-F,=f.—F., t.4) 
which can only be true if in (2) we have 
fi. = te + f, - (4.5) 


Now repeat this procedure for the pair {f, , F,}. It follows, from the foregoing argu- 


ment, that there must be another pair {f, , F,} with f, — F, = f, — F, = f, — F,. 
Furthermore, p + s, and p # k, since we have 

he _ tis -+ ley + ft . (4.6) 

Proceeding in this way, we exhaust the set of points 7 = 1, 2, --- , N — 1, with the 

result that one of the terms in the continued equality above must be fy — Fy = 0. 


Hence f; = F; fort = 1,2, ---,N — 1. 

5. Approximation in policy space. Let us now turn to the problem of determining 
an algorithm for obtaining the solution of the system in (3.2). The basic method is that 
of successive approximations. We choose an initial sequence {f{”}, and then proceed 


iteratively, setting 


_ Min (t;; + f;°'], += 1,2,---,N—1, 
. Da) 
ty’? =0, 
for k = 0, 1, 2, Ve 
The choice of {f;’} seems to require some care. Let us then invoke the concept of 


approximation in policy space in order to obtain a sequence which is monotone increasing. 
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Perhaps the simplest policy that we can employ is to proceed directly from 7 to N. 
Define 
fi” = tin ’ t= 1, 2, Sos N. (5.2) 


It follows that f{” as defined by 
fi? = Min [é,; + f;°], ¢=1,2,---,N—1, 
~ (5.3) 
fy = 0, 
satisfies the inequality 


gs < ge i= 8 2, eee N. (5.4) 


This inequality is immediate when we realize that f{"? represents the minimum 
time for a path with at most one stop. It follows then inductively that the sequences 
{f!} as defined by (5.1), with f;° as in (5.2), satisfy the inequalities 


* < ay ” = l, 2, eee N, k = 0, l. 2, hs (5.5) 


‘ , 

It is important to note that there are many other policies we could employ to obtain 
monotone convergence.* 

It follows that 

lim f;° = f,, +=1,2,--- ,N, . (5.6) 
eo 
furnishing a solution to (3.2). 

It is clear from the physical interpretation of this iterative scheme that at most 
(N — 1) iterations are required for the sequence to converge to the solution. 

6. Computational aspects. It is easily seen that the iterative scheme discussed above 
is a feasible method for either hand or machine computation for values of N of the 
order of magnitude of 50 or 100.T - 

For each 7, we require only the column (¢;;), 7 = 1, 2, --- , N of the matrix T. Hence, 
the memory requirement far a digital computer is small. 

7. Monotone increasing convergence. Turning back to (3.2), let us consider the 
sequence of approximations defined by 


f;” = Min ¢,; , ~=1,2,---,N=1, 


jy = 9, (7.1) 
fi**) = Min [t,;; + f;’], i=1,2,---,N-1, 
fy =0. 


It is clear that f{**” > f{*. It is not, however, obvious that this method yields a 


- 


bounded sequence. To establish this, let us show, inductively, that 


rs tes i=1,2,---,N, k=0,1,2,-::, (7.2) 


uniformly 


where |{/,} is the solution of (3.2). 


*I owe this choice of an initial policy to F. Haight. 
tAdded in proof (December 1957): After this paper was written, the author was informed by Max 
Woodbury and George Dantzig that the particular iterative scheme discussed in Sec. 5 had been ob- 


tained by them from first principles. 
I 
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The inequality is certainly true for k = 0. Hence, assuming that it holds for k, we 
have 


Min [t;; + f;°] < Min [t;; + f;] 


| 
om | jms 


< fi 


(7.3) 


It follows that the sequence i } converges to {f;} ask— « , furnishing the desired 


monotone convergence. Once again, only a finite number of iterations will be required. 
It is to be expecte d that the first method will converge more rapidly. 
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ON THE DAMPED OSCILLATIONS EQUATION 
WITH VARIABLE COEFFICIENTS 


By E. V. LAITONE (U7 ty of Califor? Berkele 


A useful uppel bou dad can be civen for the oscillatory solutions ot the a nd-order 


linear differential equation with continuous differentiable coefficients 
Ned & tT Dp t u'(t t qoull U l 


which commonly occurs in vil 


ration studies and airplane or missile dynami S 
The solutions of I q 1) are oscillatory whenever 


(1) qe = p $ — p’ 2) =m > U 2 
see Bellman [1] or Kamke 2], and for this case it will be shown that as long as ¢ is mono- 
tonic the solutions are bounded in the following manner: 

nt 
< A&(O)u(O)~ 4 [u’(O + p(O)u(O)/2]"} pam 6 | p 2 dt) 3) 
Wt \ 0 


This result can be derived by introducing the unique transformation, see Bellman 
[1] or Kamke [2] 


u = VU exp # | p/2 dt) 4) 


70 


which preserves the same zeros in the oscillatory solutions. Then Eq. (1) is transformed to 


v(t) + o(év(t) = 0, (5) 


/ 


*Received February 5, 1957. 
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where ¢ is defined by Eq. (2). The inequality ¢ = m* > 0 insures that all the solutions 
are oscillatory, see Bellman [1] or Kamke [2], also it is known that if in addition ¢ is 
monotonic, then the amplitudes of v also vary monotonically, increasing when ¢ is 
decre: sing and vice-versa. 

Then the function 


F(t) = gv’ +0” (6) 


is introduced, see Kamke [2], in order to develop the inequality 


F(t 2qvv’ + ov + 2v'v"" ? ov “ ¢'(t) (7) 
F(t gu +0" qv +0 = oft) 


which will be used to derive Eq. (3). 
C ase / QO monotonically decre asing. 


Since ¢’ < O, therefore F’ = g'v < (so that 


F(0) = F(t) = o(tvo(t)? = m’r(t)’. (8) 


7 | , eee 
v(t) | < — F(O)'? = — [o(Ov(0)’ + v'(0)?]'” (9) 
m m 
which transforms to Eq. (3) upon introducing Eq. (4). 
Case II, @ monotonically increasing. 
Since ¢’ > 0. therefore ¢(t) > (0) = m* and F’ = ¢’v" = Oso that 


F(t) = F(O) = ¢(0)v(0)? = m’v(0)’. (10) 
Consequently Eq. (7) may be integrated to 
¢(t) F(t) . o(t)v(t)* 
o> In > n— Ss 11 
In 50) =o F(0) — ” 5(0)v(0) + v (0)° (11) 
Chis shows that if “ l and 
¢(0) > F(0) > 0 (12) 


then Eq. (11) again yields Eq. (9), and similarly, Eq. (3). Now F > 0 for ¢ > 0 since 
vy and v’ cannot be simultaneously zero, however ¢@ > F only for sufficiently small initial 
displacements from equilibrium so that v* < 1 and v” < (1 — v’)¢@X ¢. 

Consequently the bound provided for Eq. (1) by Eq. (3) is always valid when 
¢(= m* > 0) is monotonically decreasing, but is only valid for sufficiently small initial 
values of w and u’ so that Eq. (12) is satisfied when ¢ is monotonically increasing. This 
latter requirement is satisfied whenever 


o(0) > {o(O0)u(O)? + [u’(0)* + p(O)u(O)/2]}*} (13) 
in which case Eq. (3) could be replaced by 


lu(t) | < 1 ¢(0)'”* exp (- | p/2 at) (14) 
0 


m 
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or in the usual situation wherein ¢(0) = m” when ¢ is monotonically increasing, by 
at 
| u(t) | S exp (— | p/2 at). (15) 
0 


The limiting behavior of Eq. (5) as t — , see Bellman [1] or Kamke [2], shows 
that in both cases the solutions of Eq. (1) can be written for t > © as 


u— C, sin (at + C.) exp € [ p/2 at) , (16) 


whenever ¢(/) — a°(= m* > 0) ast > © and 


Application to dynamic stability analysis of a missile. As a typical application of 
the upper bound provided by Eq. (3) let us consider the case of a missile having a nearly 
vertical trajectory with no control forces or thrust misalignment. In this case the equa- 
tions presented in [3], which define the variation of the small incremental angle of attack 
a, during constant acceleration, may be written as 


a’’(t) + p(a’(t) + q(da(t) = 0. (17) 


Therefore as long as p and q are such that the ¢ defined by Eq. (2) is monotonic, the 
upper bound of a can be determined. This is generally true for constant atmospheric 
density since the forward velocity U is continually increasing while the rate of accelera- 
tion is generally small. However if the magnitude of the acceleration is sufficiently large, 
or the height above the earth’s surface is rapidly varying, then it is necessary to trans- 
form the independent variable from the time to the distance along the trajectory, as 
in [3], before Eq. (3) can be applied. 

A particularly interesting case occurs for the vertical entry of a missile into the 
earth’s atmosphere at hypersonic speeds. Then the ratio of the acceleration to the 
velocity squared is nearly constant so that in terms of x, the non-dimensional distance 
along the trajectory normal to the earth’s surface, the variation of the small incremental 
angle of attack may be approximated to the first order by 


a’’(x) + bp(x)a’(x) + cp(x)a(x) = 0, (18) 
where b and ¢ are constants and p is the mass density of the air. 
Consequently, Eq. (2) becomes 
o(x) = [ep — (bp/2)* — bp’/2] 
and it is easily shown that ¢(x) is always bounded and monotonically increasing for 
any physically possible variation in air density as long as ¢(0) > 0 for the constant 


mass missile at hypersonic speeds. 
For example, the usual exponential approximation to the density variation may be 
written in terms of p. , the sea-level mass density of the air, as 


p(x) = pao{l — exp [—B(x + xo)]} 
so that ¢ is monotonically increasing and bounded by 


0 < ¢(0) = d(x) < [cp. — (bp 2)? — bBp../2] 
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since 
po K 1 lb. sec.’ ft.~* 
o<2<c18h.", 
for any physically possible atmosphere. 
Therefore a(x) is oscillatory as long as ¢(0) > 0, and the upper bound is defined by 
Eq. (3) as long as 
d(0) > {¢(0)a(0)? + [a’(0)? + bp(O)a(O)/2)°}, 


that is, provided Eq. (13) is satisfied. 
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BOOK REVIEWS 
(Continued from p. 72) 
The mathematics of diffusion. By J. Crank. Oxford University Press, New York, 1956. 


v1 + 347 pp. $8.00. 


n the diffusion-heat conduction equation to appear since the publication of 
slaw and Jaeger. Almost all the situations treated involve the transport of 


This is the first book « 
the well known work by ( 
matter rather than of heat, 


and many show a number of features, such as concentration dependence of 
coefficient, which are normally unimportant in the corresponding heat conduction problems. 





the diffusior 
The diffusion processes are for the most part one dimensional, although quite a few cases with spherical 


and cylindrical symmetry are presented. 
The first six chapters constitute a necessary review of fairly standard diffusion problems, but from 
this point on most of the material has never been presented in book form before; indeed, much of it is 
from the author’s own publications. The chief topics discussed are moving boundary problems, diffusion 
with simultaneous chemi reaction, concentration-dependent diffusion coefficients, and numerical 
methods. In addition there is a chapter dealing with the different kinds of diffusion coefficients measured 
by various experimental procedures, plus one on the somewhat specialized problem of coupled diffusion 
heets ol wool fiber. 
the book is the inclusion of a very large number of plots of concentration- 
distance and sorption-time relationships; these greatly aid the reader’s intuition, and many of them 
should be useful in actual computations. Also mentioned should be the appendix, which contains tables 
of a variety of useful quantities, such as integrals of the error function, as well as sample computation 
» numerical methods presented in the text. 
way of criticism, it is not about the material which has been included, 


and heat conduction in s 


One ol t he best feature Ss ol 


sheets for some of the 

If there is anvthing to be aid by 
but rather about that which has been left out. Thus, coupling between the diffusion of different species 
nent systems is not mentioned, and no problems involving convection or external forces 


in multicor 
I sion of these and other topics would have made the book excessively long. 


are treate Perhaps the 
In any case, it is a valuable contribution as it stands, and no one engaged in or contemplating diffusion 
measurements should be unfamiliar with its contents. 


STEPHEN PRAGER 


Introd ction to finite mathematics. By John cx, Kemeny, " B Laurie Snell, and Gerald L. 
Thompson. Prentice-Hall, Inc., Englewood Cliffs, N. J., 1957. xi + 372 pp. $5.00. 


Here is 0k which n safe ly be described as unique among college mathematics texts 
In recent years there has been an increasing demand for “‘terminal’’ courses in mathematics, running 


for one or at most two semesters and intended to provide appropriate material for students who may 
not intend to go furthe1 thematics. The usual introductory analytic geometry and calculus course 


nel 


because one or two terms is hardly 


is from some points of vi ins itisfactory for this purpose, firstly, 
enough time to get start the subject, and secondly, the material lacks motivation without a con- 


ces. ‘‘Introduction to Finite Mathematics’’ over- 


siderable at II OI ap} tion to the phy sical scien 
lifficulties, first, by taking for its mathematical subject matter not analvsis but 


comes both of these difficulties, 
branches of algebra), and second, by choosing motivating 


“finite mathematics,’ tually certain 
material from the social 1 than the physical sciences, thus dealing with matters which are close to 
the students’ every day ¢ perience 
to introduce a student to some concepts in modern 


} ] ? | 


In the ¢ thors own word 
mathematics early in his college career. While primarily a mathematics course, it was to include appli- 


nd social sciences and thus provide a point of view, other than that given by 


1e book’s purpose is * 


cations to the biological 
physics, concerning the possible uses of mathematics. Our aim was to choose topics which were initially 
close to the student’s experience, which are important in modern day mathematics, and which have 


interesting and important applications. 
As to the somewhat mystifying title of the book, the authors explain 


7 


‘our purpose. . . was to develop 
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several topics from a central point of view. In order to accomplish this on an elementary level, we re- 
stricted ourselves to consideration of finite problems, that is, problems which do not involve infinite 
sets, limiting processes, continuity, etc.” 

The program of the book in very brief outline, is the following: The first chapter “Compound 
Sentences,”’ deals with the elementary propositional calculus, the central idea being that of truth tables. 
Also in this chapter the notion of the “logical possibilities” associated with a statement is introduced. 
This concept which appears to be a device invented by the authors (we return to this later) provides 
a link between the first and second chapter entitled “Sets and Subsets.”’ Its main concern is the Boolean 
Algebra of sets and its relation to the propositional calculus of the previous chapter. The set notion 
leads naturally to Chapter ITI, ‘“‘Partitions and Counting,” which is concerned with elementary combina- 
torics, permutations, etc. The material of all the preceding chapters is brought together in Chapter IV, 
“Probability Theory,”’ which, however, is quite a departure from the usual treatments of this subject 
in first year texts as indicated by the presence of sections entitled “Finite Stochastic Processes,” “The 
Law of Large Numbers,” “Markov Chains.” This last section leads to Chapter V on “Vectors and 
Matrices’ which, after introducing the standard terms and operations, applies them to the theory of 
Markov chains. These five chapters comprise the “basic core” of the book and represent material which 
could be conveniently covered in one college semester. The book concludes with two “optional” chapters, 
(so indicated by an asterisk). The first entitled “Linear Programming and the Theory of Games,’ 
contains an elementary treatment of some aspects of these two subjects. The final chapter ‘‘Applications 
to Behavioral Science Problems” treats five problems, ‘‘selected for their interest both to mathematicians 


and to behavioral scientists. One topic was chosen from each of five sciences; sociology, genetics, psy- 
cholog inthropology, and economics.” 
Chis book represents such a radical departure from traditional college texts that it would take 
lerable space to give it a really adequate review. Perhaps some notion of its unusual flavor can be 
conveyed by listing a few of the many (over 1000) varied and excellent problems which appear at the 
end f each section. 
In a starred (optional) section of Chapter I on “Statements Having Given Truth Tables’’ the 
ies of exercises to the proof that every logical connective can be expressed in 
terms of ‘‘neither ... nor.”’ Ina section on “Valid Arguments” the authors develop the theory of the 


r 
Z 


syllogism so the student has a systematic procedure for checking the validity of arguments like ‘Father 

praises me only if I ean be proud of myself. Either I do well in sports or I cannot be proud of myself. 

If I study hard then I cannot do well in sports. Therefore, if Father praises me then I do not study hard.” 

The final section (starred) of this chapter is on “Application to Switching Circuits’’. Problem: 

\ committee has five members. It takes a majority vote to carry a measure except that the chairman 

to. Design a circuit for the committee, so that each member votes for a measure by pressing 
ton,.and a light goes on if and only if the measure is carried.” 

In the chapter on sets the following problem is typical. “In a survey of 100 students the number 
studying various languages were found to be: Spanish, 28; German, 30; French, 42; German and French, 
5; all three languages, 3. (a) How many students were studying no language? (b) How many students 
had French as their only language?” ete. In the chapter on partitions we find, “On a transcontinental 

irliner there are 9 boys, 5 American children, 9 men, 7 foreign boys, 14 Americans, 6 American males 


les. What is the number of people on the plane?” 


and 7 foreign fema 

(s a final example we quote one of the problems on Markov chains. “A professor tries not to be late 
for class too often. If he is late one day, he is 90% sure to be on time the next. If he is on time, then the 
next day there is a 30% chance of his being late. In the long run how often is he late?”’. 

The foregoing examples show the way in which each new mathematical concept is illustrated by a 
non-trivial application. Even such seemingly mechanical matters as the rules for matrix multiplication 
are illustrated by realistic examples in which multiplying matrices is precisely the natural thing to do. 
Indeed, of the text’s many virtues, perhaps the outstanding one is the remarkable balance which the 
wuthors have preserved between abstract theory and “common sense’’ applications. This is “applied 
mathematics’’ in the best sense. 

\ second notable feature of the book is its emphasis on pedagogical devices of a graphical nature. 
Simple problems in set theory are solved by the use of Venn diagrams. Another device which occurs 
in the first chapter and is used especially effectively in the probability chapter is that of ‘‘three diagrams” 
to represent either partitions, or “logical possibilities” at various stages of some process. This use of 
“picture writing” is thoroughly exploited to the great advantage of the student. 
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Having pointed out some of the books strong points, we now turn to what, in the reviewer’s opinion, 
is a very serious flaw, The authors have chosen as their central notion the concept of a “‘statement.”’ 
Almost all other notions in the book are derived from or closely related to this one. In the chapter on 
probability the fundamental notion is not the usual idea of the probability of an event but the “prob- 
ability of a statement.’’ It is therefore somewhat of a shock to find that this notion is used in a completely 
inconsistent manner. Specifically, the book opens with the following sentences, ‘‘A statement is a verbal 
or written assertion. In the English language such assertions are made by means of declarative sentences.”’ 
But on page 21 we find, ‘Hence, also, the statement [the italics are the reviewer’s] ‘(Draw three white 
balls” is logically false.” (Problem for the authors: what is a declarative sentence?). Next, on page 2 
we find “The fundamental property of any statement is that it is either true or false (and that it cannot 
be both true and false),”” but on page 19 we read “Normally, for a given statement there will be many 
cases in which it is true and many in which it is false.”” Which sentence is the student to believe? Worse, 
what is the instructor to tell him? This example does not represent a momentary lapse on the authors’ 
part, but a deliberate ambiguity which is carried through the whole book. 

There is no need to dwell on the morass of inconsistency caused by the authors’ refusal to distinguish 
between sentences like ‘“‘z? = 4’’ and “‘p implies qg’”’ which are not statements (they are ‘“‘sentential func- 
tions’) and are true in some cases, false in others, and honest to goodness statements like ‘‘Socrates is 
aman” and “2 + 2 = §” which are either true or false,—no cases about it. The result of the ambiguity, 
however, is to put the instructor in an unhappy dilemma. In the authors’ phraseology, he will be faced 
with the follawing “logical (or pedagocial) possibilities’ (we refrain from drawing the tree diagram) 
he must either play along with the authors, hoping the students won’t detect the deception, or he must 
rewrite a substantial portion of the exposition himself. It is unfortunate that this quite unnecessary 
mix-up had to occur in what is in many ways a superlative text. Perhaps these matters can be remedied 





in an eventual second edition. 

The last two chapters of the book should really be considered as a separate unit since the level of 
difficulty, especially in the final chapter, is much higher than in the first five. The reviewer found the 
exposition here rather muddy in spots, especially the section on Este’s learning model. 

The examples. treated were, however, uniformly interesting and, I suspect, entirely unfamiliar to 
nine out of ten teachers of college mathematics. 

Finally, this book should be welcomed eagerly by college teachers in departments other than math- 
ematics, especially behavioral scientists, as it provides precisely the sort of elementary mathematical 


training which will be useful to students whose main interest lies in this field. 
Davip GALE 


Physic and mathematics. By R. A. Charpie, J. Horowitz, D. J. Hughes, and D. J. Littler. 
McGraw-Hill Book Co., Inc., New York, and Pergamon Press Ltd., London, 1956. 
x + 398 pp. $12.00. 


This book is labeled Volume One of Series One (of eight series) arising mainly as the result of the 
review articles presented at the Geneva Conference. This volume contains articles on neutron physics 
(especially cross sections) mainly in connection with the physics of fission. There are eleven chapters 
with titles and authors as follows (1) Summary of Data on Cross Sections and Neutron Yields of U*, 
U*, and Pu™°—Harvey and Sanders; (2) Resonance Structure of U**, U*, and Pu**—Egelstaff and 
Hughes; (3) Theoretical Analysis of Neutron Resonances in Fissile Materials—Bethe; (4) Techniques 
for Measuring Elastic and Non-Elastic Neutron Cross Sections—Cranberg, Day, Rosen, Taschek and 
Walt; (5) The Cross Section of the Fission Product Poison z, as a Function of Energy—Bernstein 
6) Resonance Capture Integrals—Macklin and Pomerance; (7) Delayed Neutrons— 


and Smith; 
Homogeneous Critical Assemblies—Callihan; (9) The Physics of Fast Reactors—Codd, 


Keepin; (8 
Sheperd and Tait; (10) Heterogeneous Methods for Calculating Reactors—Feinberg; (11) Highly En- 
riched Intermediate and Thermal Assemblies—Hurwitz and Ehrlich. 
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SURVEYS IN APPLIED MATHEMATICS 


During the next six months, John Wiley & Sons will publish a series of five survey 
volumes of interest to mathematicians, engineers, physicists, and mathematical astrono- 
mers. Prepared at the request of the Office of Naval Research, this series will provide re- 
search workers in various fields of applied mathematics with an up-to-date account of 
modern progress and an appraisal of future promising research directions. The volumes 
concentrate to some extent on Russian contributions to the literature. The authors, with- 
out exception, are internationally recognized authorities in the areas of applied mathe- 
matics covered by their surveys. 


The Series . 


Vol. I—ELASTICITY AND PLASTICITY 


By J. N. Goodier and P. G. Hodge, Jr. Approx. 160 pages. In press. 


Vol. II—DYNAMICS AND NONLINEAR MECHANICS 


By E. Leimanis and N. Minorsky. Approx. 278 pages. In press. 


Vol. III—MATHEMATICAL ASPECTS OF SUBSONIC AND 
TRANSONIC GAS DYNAMICS 


By L. Bers. Approx. 278 pages. In press. 


Vol. IV—SOME ASPECTS OF ANALYSIS AND PROBABILITY 


By I. Kaplansky, E. Hewitt, M. Hall, Jr., and R. Fortet. Approx. 320 
pages. In press. 


Vol. V—NUMERICAL ANALYSIS AND PARTIAL DIFFERENTIAL 
EQUATIONS 


sy G. E. Forsythe and P. C. Rosenbloom. Approx. 286 pages. In press. 


Write to reserve your examination copies 


JOHN WILEY & SONS, Inc. 440 Fourth Avenue, New York 16, N. Y. 


Se 


LOGICAL DESIGN of 
DIGITAL COMPUTERS 


By MONTGOMERY PHISTER, JR., Thomp 
son-Ramo-W ooldridge Products, Inc. An introduc 
tion to the art and science of designing digital sys 
tems, offering a thorough description of the logi 
cal-equation school of computer design. Describes 
and interprets the work of various men in the 
field, and applies their methods and approach to 
a wide variety of problems in the logical design 
of digital computers. 1958. 408 pages. Prob. $8.50. 


QUEUES, INVENTORIES and 
MAINTENANCE 


By PHILIP M. MORSE, Massachusetts Insti 
tute of Technology. Describes the application of 
queuing theory to maintenance and inventory 
problems, treats general queuing problems, and 
discusses the effect of changes of arrival and 
service distributions on queuing results. The first 
volume in the series: Publications in Operations 
Research. 1958. 202 pages. $6.50. 


GAMES and DECISIONS: 
Introduction and Critical Survey 


By R. DUNCAN LUCE, and HOWARD 
RAIFFA, both of Harvard University. While 
keeping complex detail to a minimum, the book 
describes the mathematical attempts to deal with 
such problems as value (or utility), subjective 
probability, conflict of interest (game theory), 
decision making under uncertainty, and social 
resolutions of value conflicts. 1957. 509 pages. 
$8.75. 
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DIGITAL COMPUTER PROGRAMMING 


By D. D. McCRACKEN, General Electric 
Company. The first full-length, general introduc- 
tion to programming, with emphasis on actual 
work with computers. The book is entirely self- 
contained and most of the material covered can 
be applied to any computer. One of a series writ- 
ten by General Electric authors for the advance- 
ment of engineering practice. 1957. 253 pages. 


$6.50. 


ELEMENTARY THEORY of 
ANGULAR MOMENTUM 


By M. E. ROSE, Oak Ridge National Labora- 
tory. A clear presentation of the ideas underly- 
ing the theory of angular momentum and the 
means of using this theory. The author is con- 
cerned with the properties of rotations due to 
their intimate connection with this theory. A 
publication in Wiley’s Structure of Matter Series, 
Maria Goeppert Mayer, Advisory Editor. 1957. 
248 pages. $10.00. 


VECTOR ANALYSIS 


By LOUIS BRAND, University of Houston. 
Supplies the basic tools of vector algebra and 
calculus and gives an introduction to their appli- 
cations. All theory is illustrated by examples and 
reinforced with numerous problems. The final 
chapter on vector spaces, with its postulational 
approach, shows the nature of some modern 
generalizations of the vector concept. 1957. 282 
pages. $6.00. 


MATHEMATICS and WAVE MECHANICS 


By R. H. ATKIN, Northern Polytechnic, Lon- 
don. Sorting out the relevant mathematics in- 
volved, the book deals with a large number of 


theoretical and practical problems of atomic 
physics, both physical and chemical. 1957. 348 
pages. $6.00. 


Send for examination copies today 


JOHN WILEY & SONS, Inc. 


440 Fourth Avenue, New York 16, N. Y. 





THIRD U. S. NATIONAL CONGRESS OF APPLIED MECHANICS 


Brown University, June 11-14, 1958 


The congress program, which is to be distributed in May, will feature over 
one hundred papers in Applied Mechanics, including mechanies of rigid bodies, 
mechanics of deformable solids, mechanics of fluids and gases, thermodynamics, 
and heat transfer. In addition, there will be four general lectures by leading 


authorities. 


To reduce the cost of attending the symposium to the utmost, the following 
all inclusive living arrangement is offered. This covers rooms in the new West 
Quadrangle of Brown University, and all meals from Breakfast on Wednesday, 
June 11, through Luncheon on Saturday, June 14, but NOT the Banquet on Thurs- 


day night. The rates will be $24.00 for each of two persons sharing a room, or 


$30.00 for a person occupying a single room. Rooms will be available without 


additional charge for Tuesday night and Saturday night. Please note that this 
arrangement was only made possible by excluding all deviating arrangements for 
those planning to live on the Brown campus. Those planning to stay at Hotels or 


Motels will be able to lunch in the Brown Refectory at the price of $1.50 per 
; j ] 


luncheon. The banquet ticket will be $5.00. 
Reaistration Fees 


Advance Registration on or before June 7th, $4.00. 


Registration on June 8th or later, $6.00 


Please make checks payable to: APPLIED MECHANICS CONGRESS, not to any 


person. 


Requests for copies of the program and letters concerning advance registration 
should be addressed to Miss E. M. Addison, Box F, Brown University, Provi- 
dence 12, Rhode Island. In registering, please indicate number of persons in party 
and family relationship, desire for single rooms or willingness to share rooms, 


estimated arrival and departure times, and need for parking space. 
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